Granular Electronic Systems 

I. S. Beloborodov, A. V. Lopatin, and V. M. Vinokur 

Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 
K B. Efetov 

Theoretische Physik III, Ruhr-Universitat Bochum, 44780 Bochum, Germany and 
L. D. Landau Institute for Theoretical Physics, 117940 Moscow, Russia. 

A granular metal is an array of metallic nano-particles imbedded into an insulating matrix. Tuning 
the intergranular coupling strength a granular system can be transformed into either a good metal 
or an insulator and, in case of superconducting particles, experience superconductor-insulator tran- 
sition. The ease of adjusting electronic properties of granular metals makes them most suitable 
for fundamental studies of disordered solids and assures them a fundamental role for nanotechno- 
logical applications. This Review discusses recent important theoretical advances in the study of 
granular metals, emphasizing on the interplay of disorder, quantum effects, fluctuations and ef- 
fects of confinement in formation of electronic transport and thermodynamic properties of granular 
materials. 
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A. Why are granular electronic systems Interesting? 

Granular conductors form a new class of artificial ma- 
terials with tunable electronic properties controlled at 
the nanoscale and composed of close-packed granules 
varying in size from a few to hundred nanometers (of- 
ten referred to as nanocrystals) . The granules are large 
enough to possess a distinct electronic structure, but 
sufficiently small to be mesoscopic in nature and ex- 
hibit effects of quantized electronic levels of confined elec- 
trons. Granular conductors combine the unique proper- 
ties of individual- and the collective properties of coupled 
nanocrystals opening a new route for potential novel elec- 
tronic, optical, and optoelectronic applications. Applica- 
tions range from light emitting devices to photovoltaic 
cells and bio-sensors. The intense interest is motivated 
not only by the important technological promise but by 
the appeal of dealing with the experimentally accessible 
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model system that is governed by tunable cooperative 
effects of dis order, electron co r relations, and quantum 
phen omena llGanonenk4 Il998t iMowbrav and SkolniclJ. 
l2005HMurrav et QiJ . [l99l . 

Among traditional methods of preparation of such ma- 
terials, the most common are the thermal evaporation 
and sputtering techniques. During those processes metal- 
lic and insulating components are simultaneously evapo- 
rated or sputtered onto a substrate. Diffusion of metal- 
lic component leads to the formation of small metallic 
grains, usually 3-50 nm, see Fig. ^ Variations of grain 
sizes within a sample produced by these methods can 
be ac hieved to as low as about ~ 10% llGerber et all 
Il997j) . Depending on the materials used for the prepara- 
tion, one can obtain magnetic systems, superconductors, 
insulators, etc. 

Recent years have seen a remarkable progress achieved 
in the design of granular conductors with the controllable 
structure parameters. Granules can be capped with or- 
ganic (ligands) or inorganic molecules which connect and 
regulate the coupling between them. Altering the size 
and shape of granules one can regulate quantum con- 
finement effects. In particular, tuning microscopic pa- 
rameters one can vary the granular materials from be- 
ing relatively good metals to pronounced insulators as a 
function of the strength of electron tunnelling couplings 
between conducting grains. This makes granular con- 
ductors a perfect exemplary system for studing metal- 
insulator transition and related phenomena. 

One emerging technique to create granular sys- 
tems is through self-assembling colloidal nanocrys- 
tals llCollier et all '1998; 'Mur rav et aH l2000|) . For in- 
stance ifLin et describe almost perfectly peri- 

odic two-dimensional arrays of monodisperse (with grain 
size variations within ~ 5%) gold nanoparticles covered 
by ligand molecules that play a role of an insulating 
layer. Such samples are produced via self-assembling 
of gold nanoparticles at the liquid-air interface during 
the e vaporation of a colloidal droplet (Naravanan et al.. 
l2004f) . By changing the experimental conditions, 
multilayers of na nocrystals c an al s o be created, see 
Fig. 121 ( Parthas arathv et all l2004 l200lt iTran et all 
|2005). Other i mportant exarnples include Laiigmuir films 
of colloidal Ag (Collier a/.l.ll997t iDu et aZ.l. l200^ and 
arrays of seni i conductor quantum dots llDu et all 2002 ; 



Murrav et all . Il993t IWehrenberg et all l2002t lYu etd. 
20041) . 



Of the range of techniques developed for the fabri- 
cation of semiconductor quantum dot arrays, the most 
successful appeared the self-assembled technique of epi- 
taxial growth of two semiconductors having signifi- 
cantly different lattice constants. For the prototype 
system of In As on GaAs, where the lattice mismatch 
is 7%, In As initially deposited on GaAs grows as a 
strained two-dimensional layer (referred to as the wet- 
ting layer) jMowbrav and Skolnick, 2005). These mate- 
rials along with In20^ : Sn (known as indium tin ox- 
ide, ITO) are the most widely studied systems. They 





FIG. 1 Scanning electron microscope photo graphs of indium 
evaporated onto 5*402 at room temperatures llYu et aTlllQQlT) . 




FIG. 2 Transmission electron micrographs showing periodic 
granular (a) bilayers, (b) trilayers, (c) tetralayers, and (d) 
thick films. The insets on the left sides of panels (a) and (c) 
are the zoomed-in images. The scale bars co rrespond to 200 
nm (a)-(c) and 40 nm (d) ijTran et all\200^ 



exhibit a high visible transparency and a good electri- 
cal conductance. They are used as electrodes in light 
emitting diodes llKim et all Il998t IZhu et all l2000l) so- 
lar cells ()Gordon| . |2000() , smart windows and flat panel 
displays. 

All these experimental achievements and technological 
prospects call for a comprehensive theory able to pro- 
vide quantitative description of transport and thermody- 
namic properties of granular conductors and can, there- 
fore, serve as a ground for a clever design of devices for 
emerging new generation nano-electronics. 

It has been realized for quite some time that granu- 
larity can bring new physics extending already wealthy 
list of remarkable effects exhibited by disordered sys- 
tems even further. One of the early observation 
was the stretched exponential temperature behavior of 
the conductivity in the strongly disordered films and 
arra ys of meta llic granules in the insulating regime 
(see l)Abeles et a l. 1975) for a review): 



a(r)~exp(-VVr 



(1.1) 



This 



with To being a material dependent constant, 
behavior resembled the Mott-Efros-Shklovskii vari- 
able range hopping conducti v ity in semiconduc- 
tors ijEfros and Shklovskiil Il975t IShklovskii and Efrol 
Il988|) and appeared to be generic for granular ar- 
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rays both metallic a nd sem iconductinK - either 
irregular llAbeles et all Il975^ a nd strictly peri- 
odic llRomero and nrndid \2mFk lTra,n p.t all 
IYu et all l2004j) . Several explanations have been ad- 
vanced for this striking behav ior, but real unders t anding 
was achieved only recent ly llBeloborodov et all ^Q^^ 
iFeiKel'man and loselevich", l2005t iTran et all l2005t 
[zhang a nd Shklovs kii. 200j). 

A most recent incitement to intense and deep study of 
physics of gr anular m aterials w as given by the i nfluential 
works ( Gerb er et "d^Ximi; Si mon et oil 119871) where a 
logarithmic dependence of the conductivity a (T) 

a{T) = a + b\nT, (1.2) 

with a and b being material dependent constants, 
was observed in the metallic conductivity domain. 
This logarithmic behavior, has been observed in 
both two and three dimensional samples thus rul- 
ing out the tem pting explanation in terms of the 
weak localization l)Abrahams et all Il979l: iGorkoy et al. , 

or interaction co r rectio n s llAltshulerTT'ainTgSOt 
Belitz and KirkpatrickL Il994t iLee and Ramakrishnan , 

which result in logarithmic behavior in two di- 
mensions only. The call for u nderstanding these r esults 
brought to life ne w models fBeloboro dov et all . l2003l 
|2Q04b; Efetov and T schcrsich. 2002, 200^ which evolved 
into a new direction of research that will be one of the 
major topics of our review. 

The feature that plays a fundamental role, especially 
at low temperatures, is the pronounced discreteness of 
the electronic levels due to electron confinement within 
a single grain. The mesoscopic scale of the grains brings 
about the levels statistics and all the wealth of the related 
effects erin, 1986.; Nagaev. _1992) . The mean level 

spacing 5 in a single grain is inversely proportional to the 
volume of the grain, 

S={iyV)~\ (1.3) 

where V is the volume of the grain and v is the density 
of states at the Fermi energy. For metal particles of size 
of several nanometers, the parameter S is typically of 
order of several Kelvin. For example, for an aluminum 
particle with radius R = 5 nm one has S ~ IK. However, 
we concentrate on the temperature range T > S where 
the quantum size effects are not important. In fact both 
types of conductivity behavior, Eqs. (|1.1|) and H1.2|l . were 
observed at temperatures T > S. 

Superconductivity brings yet another aspect to 
diversity of the inter es ting ef fects in granu la r ma- 
terials (iGerber all . Il997t iHadacek a/.l l2004 
iJaeeer et all Il989t IShaoira and Deutsched . Il983|) 
One of them is a counterintuitive suppression of the 
conductivity due to superconducting fluctuations. 

In this review we summarize the recent theoretical 
progress in understanding of the phenomena observed 
in granular metals and superconductors. These are ef- 
fects that do not demand ultra low temperatures (and 




FIG. 3 a) Sketch of a nanocrystal monolayer and in-plane 
electrodes, b) Highly ordered superlatti ce between electrodes 
visib le at the upper left and lower right. iParthasarathv et all 
1200 ID 

thus very sophisticated experimental setups) and are thus 
promising from the standpoint of possible applications. 
Physics of these phenomena is not particularly material 
specific and our consideration will be based on the cor- 
respondingly general models of the granular systems. 



B. Physical quantities characterizing granular materials 

This Review deals with the systems which we call a 
granular metals and which are well modelled by an array 
of the identical in size and shape but mesoscopically dif- 
ferent metallic particles, the intergranular electron cou- 
pling being described via the tunnelling matrix. The 
grain arrangements may be either periodic or irregular. 

The effect of irregularities in the grain positions and 
in the strengths of the tunneling coupling on the physi- 
cal properties of granular systems is different for metal- 
lic and insulating samples. If the coupling between the 
grains is sufficiently strong and the system is well con- 
ducting, the irregularities are not very important. On the 
contrary, irregularities become crucial in the limit of low 
coupling where the system is an insulator. As a matter of 
practice, recent advances in fabrication techniques now 
permit very regular arrays with the f luctuations in the 
granule size wi thin the 5% precision ijTran et all l2005t 
IYu et ad . l2004|) . As regular arrays promise many techno- 
logical applications the further progress in making perfect 
systems is expected. 

At the same time, disorder related to either internal 
defects inside and/or on the surface of the individual 
granules and to the charged impurities in the insulating 
substrate/matrix is unavoidable. Even in metallic grains 
with the mean free part exceeding their size, the electrons 
that move ballistically within the granules yet scatter at 
the grain boundaries irregularities; the eletfron motion 
becomes chaotic with the resulting effect equivalent to 
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the a ction of the intragranular disorder (see, e.g. l|Efetovl 
Il997j) l. This surface chaotization would have been absent 
in the atomically perfect spherical (or, say, cubical) gran- 
ules. In such ideal granules the additional degeneracies of 
the energy levels would have led to singularities in phys- 
ical quantities. However, the slightest deviations (even 
of the order of lA or of the order of the electron wave- 
length) of the shape of the grains from the ideal spheres 
or cubes would have lifted the accidental degeneracy of 
the levels (this would be equivalent to adding an internal 
disorder). Thus one can justly assume that the grains 
are always microscopically irregular. Moreover, the cou- 
pling between the grains is the source of the additional 
irregularities. We thus can adopt the model with the dif- 
fusive electron motion within the each grain without any 
loss of generality. In this model the use of expression of 
Eq. H1.3|l for the mean level spacing i5 is fully justified. 

The key parameter that determines most of the phys- 
ical properties of the granular array is the average tun- 
nelling conductance G between the neighboring grains. It 
is convenient to introduce the dimensionless conductance 
g (corresponding to one spin component) measured in the 
units of the quantum conductance e^/h : g = G/ (2e^/fi.). 
As we will see below, the samples with 5^1 exhibit 
metallic transport properties, while those with 5^1 
show an insulating behavior. 

One of the most important energy parameters of the 
granular system is the single grain Coulomb charging en- 
ergy Ec- This energy is equal to the change in the en- 
ergy of the grain when adding or removing one electron, 
and it plays a crucial role in the transport properties in 
the insulating regime when electrons are localized in the 
grains, such that the charge of each grain is quantized. 
The physics of the insulating state is closely related to 
the well known phenomenon of the Coulomb blockade of 
a single grain connected to a metallic reservoir. 

The behavior of a single grain in a contact with a 
reservoir has been discussed in m any articles and re- 
views (se e e.g. (Alcincr et^a l., 2002; Averin and Nazarovl 
[l992; Av erin and Likhar"evtll99ll) l. The main features of 
the Coulomb blockade can be summarized as follows: (i) 
If the grain is weakly coupled to the metallic contact, 
(7 <C 1, the charge on the grain is almost quantized; this 
is the regime of the so-called Coulomb blockade, (ii) In 
the opposite limit, g ^ I, the effects of the charge quan- 
tization are negligible and the electrons freely exchange 
between the granule and reservoir. 

Although the systems we consider are arrays of inter- 
connected granules rather than a single grain coupled to 
a bulk metal, a somewhat similar behavior is expected: 
In the regime of the strong coupling between the grains, 
(7 S> 1, electrons propagate easily through the granular 
sample and the Coulomb interaction is screened. In the 
opposite limit of the low coupling, g <ti 1, the charge on 
each grain gets quantized as in the standard Coulomb 
blockade behavior. In this case an electron has to over- 
come electrostatic barrier of the order Ec in order to hop 
onto the neighboring granule. This impedes the trans- 



port at temperatures T lower than Ec- 

Throughout this review we will be always assuming 
that the mean level spacing S from Eq. (|1.3() is the small- 
est energy scale. In particular, in all cases we take the 
condition Ec ^ S he satisfied. This is the most realis- 
tic condition when dealing with the metallic particles of 
a nanometer-size scale, given that the charging energy 
Ec is inversely proportional to the radius a of the grains, 
whereas the mean level spacing 6 is inversely proportional 
to the volume V and also taking into account the high 
density of energy states in metals. Note that this may 
not be the case in arrays of semiconductor dots, where 
Ec and S may appear of the same order, but this goes 
beyond the scope of our review. 

Another important thing to remember is that the in- 
tergranular (tunneling) conductance g is much smaller 
than the intragrain conductance go by the very meaning 
of the notion of "granular system" : 

9 < 50- (1-4) 

The intra-grain conductance go is brought by scattering 
on impurities or on the boundaries of the grains and the 
inequality (|1.4|l means that the grains are not very dirty. 
The case g ^ go can be viewed as a homogeneously dis- 
ordered system and we do dwell on this limit here. The 
single grain conductance go can be most easily defined as 
the physical conductance of a granule of a cubic geom- 
etry measured in the units of the quantum conductance 
/h. In mesoscopic physics it is customary, however, to 
relate the single grain conductance go to the single grain 
Thouless energy Exh as go — Exh/S, where the energy 
Exh is defined as 

ETh = Do/a^ (1.5) 

and Do — Vpr/d is the classical diffusion coefficient. 
Other parameters are vp, the Fermi velocity, t, the elas- 
tic scattering time within the grains, and d, the dimen- 
sionality of the grain. The length a in Eq. (|1.5|) is the 
linear size (radius) of the grain. If the grains are not very 
dirty, such that the electrons move inside the grains bal- 
listically, the mean free path I = vpT should be replaced 
by the size of the grains 2a. The Thouless energy Eth, 
Eq. (|1.5I) , is proportional to the inverse time that it takes 
for an electron to traverse the grain. The energy Eth, 
Eq. H1.5|l . exceeds the mean level spacing S, Eq. 1)1. 3() . and 
therefore the intragrain conductance go is always large, 
go S> 1 , whereas the inter-granular conductance g can be 
either larger or smaller than unity. 

The above parameters make a full set of variables de- 
scribing properties of the normal granular metals. If 
the constituent particles are made out of the supercon- 
ductor material, a wealth of new interesting phenomena 
arises. The behavior of such a system can be quantified 
by adding one more energy parameter, the superconduct- 
ing gap A of the material of a single granule. Now, if the 
inter-grain coupling is sufhciently strong, the system can 
turn a superconductor at sufficiently low temperatures. 
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Properties of such a superconductor are not very different 
from those of bulk superconductors. 

On the contrary, in the opposite Hmit of weak coupHng 
between the granules, an array of superconducting grains 
can transform into an insulator at T = 0. In this regime 
the Cooper pairs are locally formed in each grain but 
remain localized inside the grains due to the strong on- 
site Coulomb repulsion leading to the Coulomb blockade. 
Because of localization, the number of the Cooper pairs 
in each grain is fixed and, according to the uncertainty 
principle, this leads to strong phase fluctuations. Thus 
the system does not develop the global coherence, and 
the global macroscopic superconductivity is suppressed. 
One can describe this effect using the model of supercon- 
ductor grains coupled via the Josephson junctions. Such 
a system is characterized by three energy parameters: the 
superconductor gap, A, characterizing a single grain, the 
Josephson coupling J, and the grain charging Coulomb 
energy Ec- Strong Josephson coupling, J ^ Ec, sup- 
presses the phase fluctuations leading to the globally co- 
herent superconducting state at sufficiently low temper- 
atures. If J ^ Ec, the Coulomb blockade prevails, the 
Cooper pairs get localized at T ^ 0, and the system falls 
into an insulating state. 

Note that even in the insulator state the superconduct- 
ing gap A still exists in each grain and its value is close 
to the gap magnitude in the bulk provided A ^ i5. If 
the latter inequality is not fulfilled the superconducting 
gap in t he grains is suppressed or can even be fully de- 
stroyed l)AndersonLll959iri964() . There are interesting ef- 
fects in this regime but we do not consider them here be- 
cause the relevant region of the parameters corresponds 
to either too small grains or too weak superconductors. 

At low temperatures, the Josephson coupling J is ex- 
pressed via the tun nelling c onduc tance g as J = 7rgA/2 
ijAmbeeaokar and'BaratoflL Il963() and, at first glance, 
one could have concluded that the transition between 
the insulating and superconducting states should have 
happened at g ^ Ec/ A. However, this simple estimate 
holds only in the weak coupling regime, g < 1, that as- 
sumes Ec < A. At stronger coupling, the Coulomb en- 
ergy Ec is renormalized down to the value Ec ^ E ^ 
A/(7, ijChakravartv et glLllQSTtlLarkin and Ovchinnikovl 
Il983|) due to electron tunnelling between the neighbor- 
ing grains. Therefore, in the limit of the strong cou- 
pling, g ^ 1, the effective Coulomb energy E is always 
smaller than the Josephson coupling, E < J, imply- 
ing the superconducting ground state. Shown in Fig. ^ 
is the schematic phase diagram summarizing the above 
consideration (Chakravartv et ai, 1987). For many su- 
perconducting granular samples available experimentally 
the ratio Ec/ A is large. In this case, as one can see from 
phase diagram in Fig. ^ the transition between the su- 
perconducting and i nsulating at T — » oc curs at g 1 
l)Chakravartv adlT987tlOrr et a/lll986(l . 

To conclude our brief introduction to granular super- 
conductors we note that for the most of experimentally 
available samples, the grains are much smaller than the 




FIG. 4 Schematic phase diagram for granular superconduc- 
tors at temperatures T = 0. Symbols S and / stand for 
superconducting and insulating phases respectively. 

bulk superconducting coherence length of the granule 
material 

a « ^0- (1.6) 

This allows one to neglect variations of the superconduct- 
ing order parameter A inside the grains and treat a single 
grain as a zero dimensional object. As the condition H1.6|) 
is at this point the most common experimental situation, 
we will further concentrate mainly on this regime. 

As we now see, depending on the parameters of the sys- 
tem, many different physical situations appear and the 
phase diagram of a granular material is quite reach. In 
the next two sections, Sec. Inland Sec. IIIII wc will focus 
on the systems consisting of the normal and supercon- 
ducting grains respectively. 

II. NORMAL GRANULE ARRAYS 

A. Transport properties 

We start the discussion of the properties of the granu- 
lar systems consisting of the normal grains by presenting 
the main results and their qualitative explanations. This 
may help the reader to understand the basic physics of 
the system and to learn important formulae suitable for 
a direct comparison with experiments without going into 
the technical detail of theoretical models, which will be 
presented at the subsequent parts of each section. 

1. Classical conductivity 

The key parameter that determines the transport prop- 
erties of the granular materials is the tunneling conduc- 
tance g. In the strong coupling regime, g ^ 1, a granular 
array has metallic properties, while in the opposite case, 
g <C 1, the array is an insulator. The insulating state ap- 
pears as a result of the strong Coulomb correlations that 
block electron transport at low temperatures. Generally 
speaking, apart from the Coulomb interaction effects one 
has also to consider the effects of quantum interference 
that may play an important role in the low conducting 
samples. In homogenously disordered systems the in- 
terference effects play an important role leading to the 
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local ization of electro n states in the absence of interac- 
tion l)Andersonlll958() . 

In the metallic regime and at high enough tempera- 
tures both Coulomb correlation and interference effects 
are weak. In this case the global sample conductivity ctq 
is given by the classical Drude formula, in particular, for 
a periodic cubic granular array 



to r as 



(2.1) 



where a is the size of the grains and d is the dimension- 
ality of the sample. Formula H2.HI has a straightforward 
meaning: in order to obtain the physical conductance of 
the contact, one should multiply its dimensionless con- 
ductance, g, by 2e^ (the factor 2 is due to spin and h = 1) 
and, then, multiplying the result by a^~'^ one arrives at 
the conductivity per unit volume. 

Although Eq. H2.1|) is written for the periodic array, 
the electron system is not translationally invariant, oth- 
erwise the sample conductivity would be infinite. Equa- 
tion 1)2. l|l is valid for grains with an internal disorder, 
such that the electron motion inside the grains is chaotic. 
When hopping from grain to grain, the electron momen- 
tum is not conserved and this leads to the finite con- 
ductivity ao, Eq. H2.1|l . The conductance of the contact 
g depends on the microscopic properties of the contact, 
and we will consider it in most cases as a phenomenolog- 
ical dimensionless parameter controlling behavior of the 
system. Note that the model of a periodic array assumes 
no variation in tunneling conductances. 

Upon decreasing temperature the Coulomb interac- 
tions become relevant and formula Eq. H2.1|l does not hold 
any more. Below we will discuss Coulomb effects in more 
detail and find that their manifestation in granular sys- 
tems may differ noticeably from that in "homogeneously 
disordered" metals. 



2. Metallic regime 

In the metallic regime electrons tunnel easily from 
granule to granule. The time tq that the electron spends 
inside a grain plays an important role in the metallic 
regime: the corresponding characteristic energy F = Tq"^ 
is related to the tunneling conductance and the mean 
energy level spacing as 



F = gS. 



(2.2) 



The energy F can also be interpreted as the width of the 
smearing of the energy levels in the grains. In the limit of 
large conductances, g 3> 1, this width exceeds the energy 
spacing S and the discreteness of the levels within a single 
grain ceases to be relevant. 

Since the electron motion on the scales well exceeding 
the granule size is always diffusive (even in the case of the 
ballistic electron motion inside each grain), the electron 
motion on the time scales larger than F~^ can be de- 
scribed by the effective diffusion coefficient D^ff related 



D 



(2.3) 



Then the Einstein relation gives the conductivity of the 
granular sample as 



(2.4) 



For periodic arrays Eq. H2.4|) is equivalent to Eq. (|2.1|l . 
which can be seen from Eqs. 12.21 12.3|l . At the same 
time Eq. (|2.4() is more general than Eq. H2.1f) since with 
the properly defined diffusion constant Def / it applies to 
arrays with the arbitrary grain arrangement as well. 

The energy scale F plays a very important role; many 
physical quantities have qualitatively different behavior 
depending on whether they are dominated by the energies 
higher or lower than F. 

From the experience with the homogeneously disor- 
dered metals one can envision two major causes that 
may alter the classical conduct ivity gp in Eq. I|2.1[l: 
(i) electron-elec tron interactions ('Alt shuler and Aronovl 
119851 iLee and Rain akrishnan, 1985) and ( i i) quantum in- 
terfer ence effects lAbraham s et all Il979t iGorkov et al\ . 
Il979() . Accordingly, constructing the theory of granular 
conductors with the reference of highly advanced theory 
of disordered metals, one can expect two corresponding 
distinct corrections to ctq- Since in the metallic domain 
electrons effectively screen out the on-site Coulomb inter- 
actions the bare magnitudes of which are high because 
of small sizes of the grains, the notion of the interaction 
corrections to conductivity is well justified. 

To gain the qualitative understanding of the interac- 
tion effects, we introduce the characteristic interaction 
temporal- and the corresponding spatial scales tt ''^ h/T 
and Lt ~ associated with this time. One 

expects that the behavior of the interaction correction 
is different on the distances exceeding the granule size, 
Lt > a, and within the granule, Lt < a. Using Eq. 1)2. 3|l 
for the effective diffusion coefficient Deff, one immedi- 
ately sees that these conditions imply the existence of two 
distinct temperature regions T > F and T < T with re- 
spect to interaction contributions. Accordingly, the cor- 
rection to the conductivity due to Coulomb interaction 
can be written as a sum of contributions coming from the 
large, e > F, and the low, e < F, energies. This "contri- 
bution separation" naturally follows, as we will see below, 
from the diagrammatic approach where the two contri- 
butions in question are represented by the two distinct 
sets of the diagrams. Deferring the details for later, the 
result is as follows: denoting the corrections coming from 
the high and low with respect to F energies as dai and 
6a2, respectively, we write: 



o-Q + Sai + 5(72- 



with 
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(2.6) 
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FIG. 5 Resistance of 3d granular Al — Ge sample as a function 
of temperature on a log- log scale, as measured at (zero) (x) 
and 100 kOe field ( open circles). Samp le room temperature 
resistance is 500 Q (iGerber et a/.L ^1997^. 



l|Efetov and Tschersichl l2002ll2003() and 
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d = 
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d = 
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/3 IT 
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where a ~ 1.83 and /3 « 3.13 are numerical con- 
stants ()Beloborodov et a/J .b003). The high energy con- 
tribution Sai in Eq. 1)2. 5|l contains the dimensionality of 
the array d merely as a coefficient and is, in this sense, 
universal. On the contrary, the low energy contribution 
6a2 in Eq. (|2.5|) has a different functional form for differ- 
ent array dimensionalities (note that for the 3d correction 
we have kept the temperature-dependent part only). 

At high, T > r, temperatures the correction Sai, 
Eq. (|2.5|l . grows logarithmically with decreasing tempera- 
ture. Upon further lowering the temperature this correc- 
tion saturates at T « F and remains constant at T < F. 
Then the correction Sa2 in Eq. H2.7|) from the low energy 
scales, e < F, where the coherent electron motion on the 
scales larger than the grain size a dominates the physics, 
and which is si milar to that derived for homogeneous dis- 
ordered metals (lAltshuler and AronovL Il985l) comes into 
play. In the low temperature regime it is this term that 
entirely determines the temperature dependence of the 
conductivity. At the same time, the contribution Sai, al- 
though being temperature independent, still exists in this 
regime and, as a matter of fact, can even be larger in mag- 
nitude than (5(72 • Equation 1)2. 6|l can be written for T > F 
in the form Eq. H1.2I1 that has been observed in a num- 
its llFuiimori et all [l993 iGerber et all 



ber of experiments . . 

I1997I: iRotkina et ad 120051: ISimon et aZ.LIl987|) . (see, e.g 
Fig. 13 . 

Rewriting the low energy contribution 5(72, in 
terms of the effective diffusion coefficient -De// from 



Eq. (|2.3ll one reproduces the Al tshuler-Aronov correc- 
tions (Altshuler and Aronovl [1985 ). This reflects a uni- 
versal character of a large scale behavior of a disordered 
system and we will show it rigorously in Sec. III.DI that 
indeed the granular metal model can be reduced to an ef- 
fective disordered medium on distances much larger than 
the single grain size (i.e. coming from the low energy, 
T < F excitations). The contribution 6ai, which is dom- 
inated by the energies e > F, is the consequence of and 
specific to the granularity - and does not exist in the 
homogeneously disordered metals. 

Now we turn to quantum interference (weak local- 
ization) effects that exist also in the systems without 
any electron-electron interaction. In the metallic regime 
where the perturbation theory with respect to inverse 
tunneling conductance, 1/g, the interaction and weak lo- 
calization corrections can be considered independently in 
the leading order. 

The weak localization correction is of a purely quan- 
tum origin: it stems from the quantum interference of 
electrons moving along the self-intersecting trajectories 
and is proportional to the return probability of an elec- 
tron diffusing in a disordered medium. In one- or two- 
dimensional conductors the probability of infinite times 
returns is unity and returning trajectories can be in- 
finitely long. Thus a fully coherent electron propagation 
would lead to a divergent weak localization correction. 
Finite phase relaxation (dephasing) time disables long 
trajectories and limits the correction. Adapting here the 
concept of the effective di sordered medium we can as- 
sume that the results of ijAltshuler and Aronovl Il985j) 
for homogeneously disordered metals apply to granular 
metals with the proper renormalization of the diffusion 
constant, we introduce the effective dephasing length 
Lcf, = ^yOeffT^ ~ a-y/Fr^, where is the dephasing 
time, which dtermines the scale for the interference ef- 
fects. At low temperatures the dephasing time is large 
and the decoherence length can exceed the size a of a 
single grain, > a. In this regime the relevant electron 
trajectories pass through many granules and quantum in- 
terference effects are similar to those in homogeneously 
disordered metals and contribute essentially to the con- 
ductivity. With the increasing temperature the decoher- 
ence length decreases and as soon as it drops to the 
grain size, ~ a, the trajectories contributing to weak 
localizations traverse to only one neighbor and the weak 
localization correction is suppressed. Using Eq. H2.3|l we 
write the last condition separating the domains of "rel- 
evance" and "irrelevance" of weak localization effects as 
Tff) ~ F^^. For r0 > F^^ (or > a) the quantum in- 
terference effects are important, while for < F~^ (or 
< a) the corresponding correction drops rapidly with 
temperature. 

The final result for the quan tum interferenc e cor- 
rection s reads ((Beloborodov et al~ ■2004c Biagini et al\ . 
I2005a|) 
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for granular films, and 



,1/2 



for granular wires. In both equations is the dephas- 
ing time. Within this time the wave function retains its 
coherence. The most common mechanism of dephasing, 
the el ectron-electron interac tions gives for the dephasing 
time ijAltshuler et a/.lll982fl : 



g 



(^) 



1/3 



d = 2, 
d=l. 



In this case the condition F « ""^ defines (in 2d case) yet 
another characteristic energy scale, T* = g^5, that marks 
the interval of relevance of weak localization effects. 

The quantum interference correction Sa^L are sup- 
pressed by applying even the relatively weak magnetic 
field; the dependence upon the magnetic field can serve 
as a test for identifying the weak localization effects. At 
sufhcicnt fields thus the main temperature dependence 
of the conductivity will come from the electron-electron 
interaction effects, Eqs. (|2.6|l and (|2.7|) . 

Both, the electron-electron interactions and quantum 
interference effects decrease conductivity of a granular 
system at low temperatures similarly to the same effect 
in homogeneously disordered systems. The novel impor- 
tant feature is that the granularity restrains screening 
thus enhancing the role of Coulomb interaction: this is 
reflected by the contribution 6ai to conductivity which 
is specific to granular conductors but absent in the ho- 
mogeneously disordered metals. As a result, in 3d [and, 
to some extent, in 2d, see below], granular systems the 
Coulomb interactions can become a main driving force of 
a metal-insulator transition. 

This question worth more detailed discussion. The 
"granular" contribution, Sai, comes from the short dis- 
tances and is actually the renormalization of the tunnel- 
ing conductance between the grains: 



1 



9^ 9^9 



2TTd 



In 



9Ec 



max (T, r) 



Using renormalization group methods one can show that 
Eq. (|2.11() represnts in fact the solution of the renormal- 
ization group equation for the effective conductance g 
rather than a merely perturbative correction. As such, 
it holds not only till the second term in the r. h. s. of 
Eq. H2.11|l is much smaller than the first one but in a 
broader temperature region - as long as the renormalized 
conductance is large, ^ > 1. It is important that the loga- 
rithm in Eq. (|2.11|) saturates at temperatures of the order 
of r. Then one sees from Eq. H2.11|l that the renormal- 
ized conductance g may remain large in the limit T — > 
only provided the original (bare) conductance g is larger 
than its critical value 



gc = {l/2Trd) HEJS). 



It g < gc the effective conductance g renormalizes down 
to zero at finite temperature. Of course, as soon as g be- 
(2.9) comes of the order of unity, Eq. (|2.11|l is no longer vahd, 
but it is generally accepted - in the spirit of the renor- 
malization group approach - that the conductance fiow 
to low values signals the electron localization (a recent 
exact solution for a model equivalent to a single grain 
connected t o a metallic contact lends confidence to this 
conclusion ijLukvanov and Zamolodchikovl [2004(1 '). The 
result H2.12|l can be also obtained by the analysis of the 
stability of the insulating state that we discuss below. 

(2.10) The physical meaning of the critical conductance gc is 
most transparent for a 3d system: in this case, there are 
no infrared divergencies and both the localization and 
Altshuler-Aronov corrections are not important even in 
the limit T — > 0. This implies that so long as g > gc the 
granular conductor remains metallic. 

At g < gc the conductance renormalizes to low values 
at temperatures exceeding F signaling thus the develop- 
ment of the Coulomb blockade. Upon the further tem- 
perature decrease the system resistance begins to grow 
exponentially at certain characteristic temperature Tch 
indicating the onset of the insulating behavior. The tem- 
perature Tch tends to zero in the limit g gc- Thus, in 
3d the value gc marks the boundary between the insulat- 
ing and metallic states at T ^ 0. 

The interpretation of the critical value gc is less 
straightforward in the case of granular films, since at 
g > gc the low temperature conductivity corrections due 
to interaction and localization effects diverge logarithmi- 
cally. This means apparently, that the system turns an 
insulator without a sharp transition. Yet the notion of 
gc still makes sense as a mark distinguishing between the 
systems that are strong Coulomb insulators at low tem- 
peratures {g < gc) and those that are weak insulators 
{g > 9c)- 

All the above results have been obtained within the 
model of a periodic granular array (cubic lattice) neglect- 
ing the dispersion of the tunneling conductances. In re- 
ality, the typical granular samples are disordered. It is 
then important to understand what effect the irregular- 

(2.11) ities in the granule's arrangement may have. It is plau- 
sible and intuitive that the universal regime of low ener- 
gies is hardly affected by the irregularities since all the 
physical characteristics in this regime can be expressed 
through the effective diffusion coefficient of the medium. 
At the same time, the physical results that are controlled 
by the local physics, in particular, the critical value gc 
are sensitive to the grain arrangements. In some cases 
the irregularities can be incorporated at almost no ex- 
pense. For example, breaking some finite fractions of the 
junctions between the granules merely changes the aver- 
age coordination number z, which appears as the l/2d 
factor in the expressions (|2.11l I2.12|l for the critical con- 
ductance. This means that the dispersion in tunneling 
conductances is not expected to change noticeably, for 
example, the position of the metal-insulator transition 

(2.12) (which may acquire the percolation character) and can 
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be accounted for via replacing 2d in Eqs. (|2.11l I2.12|l by 
the effective coordination number ^e// < z. 

In general the proper treatment of disorder in the grain 
arrangement may appear more tedious. Yet we do not 
expect that it can change the physics of the metallic state 
qualitatively. Irregularities of the grain displacements 
and of other quantities characterizing the system play 
much more important role in the insulating regime which 
we briefly review in the following subsection. 

3. Insulating regime 

We begin with a consideration of a periodic granular 
array which, in the regime of a weak coupling between the 
grains, is an exemplary Mott insulator at low tempera- 
tures. The electron transport is mediated by the electron 
hopping from grain to grain. However, leaving a neutral 
grain and entering its neighbor costs a considerable elec- 
trostatic energy, and the electron transport is blocked 
at low temperatures by the Coulomb gap in the elec- 
tron excitation spectrum Am- At very small tunneling 
conductances this gap is simply the Coulomb charging 
energy of the grain. Am — Ec- Virtual electron tun- 
neling to neighboring grains leads to a reduction of the 
Mott gap Am and, in the limit of noticeable tunneling, 
it decays exponentially in g (Bcloborodo v et ai, 2005c;3 
until it reaches the inverse escape rate from a single grain 
r. At Am ~ r the system falls into a regime of weak 
Coulomb correlations and the insulator-metal transition 
occurs in 3d. Using the estimate Am T for estimating 
the transition point one arrives, within the logarithmic 
accuracy, at the same result for the critical conductance 
gc as that derived from a renormalization group consid- 
erations (see for details Section FlI.Fp . 

The presence of the hard gap in the excitation spec- 
trum leads to the activation dependence (Arrhenius law) 
of the conductivity on temperature: 

a ~ e"'^"/'^, T < Am. (2.13) 

Indeed, the finite temperature conductivity is due to the 
electrons and holes that are present in the system as real 
excitations. Their density is given by the Gibbs distri- 
bution that results in the exponential dependence of the 
conductivity, Eq. H2.13I) . 

However, the activation behavior is usually not ob- 
served in real granular samples at low temperatures. In- 
stead, the experimentally observed resistivity follows the 
law, Eq. Hl.l|) . that resembles the Efros-Shklovskii law 
derived for doped semiconductors. The fact that the 
observed conductivity behavior cannot be explained in 
terms of the periodic model suggests that disorder plays 
the crucial role in formation the low temperature con- 
ductivity in the insulator state. 

The stretched exponential Shklovskii-Efros-like con- 
ductivity behavior in granular conductors remained a 
challenging puzzle for a long time. Several explana- 
tions had been advanced, in particular, it was pro- 




,-1 * . Thick film 



0.08 0.12 0.16 0.20 

FIG. 6 Conductance go vs inverse temperature T~^^^ 
for per iodic gra nular multilayer and thick-film shown in 
Fig. |5| llTran et a i. 2005). Inset: for the high-temperature 
range go has been replotted as a function of , indicating 
Arrhenius behavior from 100 — 160 K. 



posed l|Abeles et all Il975(l . that the capacitance ran- 
dom variations resulting from grain size dispersion 
could provide the dependence, Eq. (11.111 . However, 
as was pointed out in Refs . ijPollak and Adkins, 199^ 
IZhang and ShklovskiiL |2004|) . the capacitance disorder 
can never lift the Coulomb blockade in a single grain com- 
pletely and therefore cannot give rise to the finite density 
of states at the Fermi level. Furthermore, the stretched 
exponential dependence, Eq. Hl.l|l . was r ecently observed 
in th e periodic arrays of quantum dots ijYakimov et al\ . 
l2003j) and artificially manufactured metallic periodic 
granular systems (Tran eJ. a /,, 2005), see Fig. El where 
the size of granules and the periodicity in the dots ar- 
rangement were controlled within a few percent accuracy. 
Either of those systems does not posses a noticeable ca- 
pacitance disorder at all, yet the dependence Hl.l|l has 
been observed. These indicates that it should be electro- 
static disorder unrelated to the grain size variations but 
caused most probably by charged defects in the insulating 
matrix/substrate that is responsible for lifting Coulomb 
blockade and formation the finite density of states near 
the Fermi level resulting in the dependence Hl.l(l . 

There is however another ingredient necessary for the 
variable-range-hopping type (VRH) conductivity to oc- 
cur (apart the finite density of states at energies close 
to Fermi level): the finite, although exponentially decay- 
ing with the distance, probability, for tunneling to the 
spatially remote - and not only to the adjacent - states 
close to the Fermi level. Accordingly, the problem of the 
hopping transport in the granular conductors is two-fold 
and should contain: (i) explanation the origin of the fi- 
nite density of states near the Fermi-level and the role of 
the Coulomb correlations in forming this density of states 
and (ii) constructing and quantitative description of the 
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mechanism of tunneling over long distances through the 
dense array of metallic grains. 

The behavior of the density of states, as we have 
mentioned above, can form due to the onsite random 
potential, which in its turn is induced by the carrier 
traps in the insulating matrix in the granular conduc- 
tors. The traps with energies lower than the Fermi level 
are charged and induce the potential of the order of / nr 
on the closest granule, where k is the dielectric constant 
of the insulator and r is the distance from the granule 
to the trap. This compares with the Coulomb block- 
ade energies due to charging metallic granules during 
the transport process. Such a mech anism was consid- 
ered by [Zhang and Sh klovskiil l)2004j) . Speaking about 
the 2d granular arrays and/or arrays of quantum dots, 
one can expect that the induced random potential origi- 
nates also from imperfections and the charged defects in 
the substrate. 

We model the electrostatic disorder via the random 
potential Vi , where i is the grain index. Such a potential 
gives rise to a flat bare density of states at the Fermi 
level. In a complete analogy with semiconductors, the 
bare density of states shoul d be suppressed by th e long- 
rang Coulomb interaction llEfros and Shklovskiil 119751 
ihklovskii and Efros, 1988). 

Next we have to consider the electron hopping over 
the distances well exceeding the average granule size in 
a dense granular array. This processes can be realized 
as tunneling via the virtual electron levels in a sequence 
of grains. The virtual electron tunneling through a sin- 
gle granule or a quant um dot (the so called co-tun neling) 
was first considered bv lAverin and Nazarovl l)l990() where 
two different mechanisms for a charge transport through 
a single quantum dot in the Coulomb blockade regime 
were identified. Namely, there are elastic and inelastic 
co-tunneling mechanisms. In the course of elastic co- 
tunneling the charge is transferred via the tunneling of 
an electron through an intermediate virtual state in the 
dot such that the electron leaves the dot with the same 
energy as it came in. In the latter mechanism (inelas- 
tic co-tunneling) an electron that comes out of the dot 
has a different energy from that of the incoming one. 
After inelastic co-tunneling the electrons leaves behind 
in the granule the electron-hole excitation absorbing the 
in- and out- energy differences. Note that both these 
processes are realized via classically inaccessible inter- 
mediate states, i.e. both mechanisms occur in the form 
of the charge transfer via a virtual state. The inelas- 
tic cotunneling dominates at temperatures larger than 
Ti ^/Ec6 (Avcrin and Nazarov, 1990). 

These two co-tunneling mechanisms can be general- 
ized to the case of the multiple co-tunneling through 
several grains. The tunneling probability should 
fall off exponentially with the dist ance (or the num - 
ber o f granules left behind N) llBeloborodov et al. 
| 2005d lFeigcrma ~and loselevichl boOSt iTran et al 
|2003) and this is equivalent to the exponentially 
decaying probability of the tunneling between the 



states near the Fermi surface in the th eory of 

Mott-Efros-Shklovsk i i llEfros and Shklovskij Il975t 
IShklovskii and Efro5lll98^ 

Thus the hopping processes in the amorphous semi- 
conductors and granular materials are alike - up to the 
specific expressions for the localization lengths - and min- 
imizing the hopping probability in the same manner as in 
the classic Efros-Shklovskii work, one ends up with the 
hopping conductivity cr in a form 



exp 



(2.14) 



where Tq is a characteristic temperature depending on 
the particular microscopic characteristics. Explicit ex- 
pressions for this temperature in different regimes are 
given in Section lil.GI Eq. 12.141) explains the experimen- 
tally observed conductivity behavior in poorly conduct- 
ing granular materials (see Eq. (|l.l|l 'l. 

Our discussion of the variable range hopping via virtual 
electron tunneling through many grains is based on the 
assumption that the hopping length r* exceeds the size 
of a single grain a. This length decays when the temper- 
ature increase and reaches the grain size at some charac- 
teristic temperature T. Then the VRH picture does not 
apply any longer, the hops occur between the adjacent 
granules only. Once the probability of a single jump is 
defined, the quantum effects can be neglected and one 
can use a classical approach. In particular, at tempera- 
tures T > T one expects the conductivity to follow the 
simple Arrhenius law. 

The classical approach to the transport in granular 
metals was developed in Refs. (iJha and Middletonl . l2005t 
iMiddleton and Wingreenl Il993() . One of the results of 
this study is the presence of the threshold voltage below 
which the conductivity is exactly zero at T = 0. The ex- 
istence of such a threshold voltage is a consequence of the 
classical approach where multiple cotunneling processes 
are not taken into account. 

In order to match the results of the classical and 
hopping t heories one has t o gen eralize the approach 
of Middlct on and Wingreenl |)1993|) to include the multi- 
ple cotunneling processes. Development of such a theory, 
in our opinion, is an interesting and important task. 



B. Model and main theoretical tools 

From the theoretical point of view a granular con- 
ductor is an appealing exemplary system whose behav- 
ior is governed by a non-trivial interplay of electron- 
electron interactions, disorder and quantum fluctuations. 
The powerful approaches that allowed for recent break- 
throughs in our understanding the physics of granular 
media are based on the effective field theories; in par- 
ticular, the phase action technique appeared especially 
suitable. For the situations that fall outside its range of 
applicability, the appropriate diagrammatic techniques, 
generalizing those for the homogeneously disordered sys- 
tems, can be developed. In this subsection we present 
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in detail a model for description of granular metals and 
introduce both complimentary methods that serve as a 
foundation for quantitative description of granular con- 
ductors. 



Another important block for the diagrammatic technique 
is the diffusion propagator (diffusion) which is just the 
impurity averaged particle-hole propagator. In bulk dis- 
ordered metals it is given by 



1. Hamiltonian 

We model the granular system as an array of metallic 
particles connected via tunneling contacts. The Hamil- 
tonian H describing a granular conductor has the form 



H — ^ Ho,i + i?/ + lit 



(2.15) 



where i/o.i stands for the Hamiltonian of non-interacting 
electrons in a grain represents the interactions, and 
Ht describes the electron tunneling between the grains. 
We now discuss each term in Eq. H2.15|l . 

The i?o,i term describes the free electrons within the 
each grain in the presence of impurities 



0,t 



(r) ( + udv) (r) dr, (2.16) 



where V'i^(r), ?/;i(r) are the electron creation and anni- 
hilation operators, \x is the chemical potential and u (vi) 
represents disorder responsible for the electron scattering 
inside the i-th grain. We adapt the Gaussian distribution 
for u (r^) with pair correlations 



(M,(r)Mj(r')) = 



1 



-5(r-r')<5^, 



(2.17) 



imp 



Throughout our review we assume that all the grains are 
in the diffusive limit, i.e. the electron mean free path I 
within each grain is smaller than the grain size a. This as- 
sumption simplifies our calculations because it allows to 
avoid considering the electron scattering from the grain 
boundaries that becomes the main "disorder" mechanism 
in the case of ballistic grains. At the same time, most of 
the results obtained in the diffusive limit are expected to 
hold also for the ballistic grains with an irregular surface 
as long as the single grain diffusion coefficient Dq does 
not enter the final result. This happens in the normal 
grains provided all relevant energies are smaller than the 
Thouless energy Etk of a single grain, Eq. (|1.5|) . Actu- 
ally, for typical grain sizes of the order 100^ the mean 
free path / is comparable with the granule size a. 

A single grain may be considered within the standard 
diagram matic appro a ch de veloped for disordered met- 
als (,Abrikosov et alY Il965|) . The main building block 
for the diagrams is a single electron Green function aver- 
aged over disorder and can be derived by the self consis- 
tent Born approximation (SCBA). The diagram shown 
in Fig. [7^ represents the relevant contribution to the self 
energy: 



Goe (p) = - C( P 



.sgn (e) 



2Ti 



imp 



(2.18) 



i^(^,q) 



(2.19) 



where Dq is the classical diffusion coefficient and uj is the 
Matsubara frequency. 

In a grain, however, the term has to be replaced by 
the Laplace operator with the proper boundary condi- 
tions ( Aleiner et al 2002 : Efetov, 1983, 1997). This pro- 
cedure leads to the quantization of the diffusion modes 
and to the appearance of the lowest excitation energy of 
the order of Dq/o? with a being the grain size. This is 
just the Thouless energy Exh that has been defined in 
Eq. fT3| . 

The Thouless energy of a nanoscale grain is large 
and this allows us to simplify calculations neglecting all 
non-zero space harmonics in the diffusion propagator, 
Eq. (|2.19|) . Throughout our review we will assume that 
Exh is the largest energy scale associated with our sys- 
tem, and will use the zero dimensional approximation for 
the diffusion propagator 



D{u)^1/\lo\. 



(2.20) 



Now we turn to the description of the electron-electron 
interaction, the second term Hi in the Hamiltonian, 
Eq. (|2.15|) . We begin with the consideration of the in- 
teraction effects in an isolated grain and then we will 
turn to the discussion of the intergranular terms that 
are important because of the long range character of the 
Coulomb interaction. 

The most general form of the electron-electron inter- 
action Hi in an isolated grain is 



p,q,r,s 



(2.21) 



where the subscripts p, q, r, s stand for the states in the 
grains and a, /3 label electron spins. 

The matrix Tipqrs is a complicated object, but in the 
low energy region of the parameters, not all of the matrix 
elements are of the same order. In the case of disordered 
grain that we consider, o nly the elements wit h the equal 
in pairs indices survive l)Aleiner et all l2002|^ , all others 
are small in the parameter 1/go, where go is the single 
grain conductance. Thus in the leading order in 1/qn the 
in-grain interaction term is l)Aleiner et all |2002|) 



H 



(0) 



(2.22) 



where n = J (r) V'i (r) dr — No is the number of exces- 
sive (with respect to the electron number in the charge 

neutral state Nq) electrons in the grain, S is the total 
spin of the grain and f t,!" are the Cooper pair creation 
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and annihilation operators: T = V'p.tV'p.i- The inter- 
action strengths in these three channels are controlled by 
coupling constants Ec,Js,^ respectively. The coupling 
constants in the spin and BCS channels, J5 and A, can- 
not be large and are, at most, of the order of 5. At the 
same time the charging energy Ec is usually much larger 
than 5. For this reason, in the absence of the super- 
conductivity, the most interesting and noticeable effects 
come from the Coulomb correlations. 

The above considerations of the interactions in a single 
grain can be easily generalized to the case of a granular 
array. It is clear that the bare interactions in the spin and 
Cooper pair channels are short-range and thus they have 
to be diagonal in the granular indices. At the same time 
the Coulomb interaction is long-range and it's of-diagonal 
components cannot be neglected. Thus, we arrive at the 
following Hamiltonian that describes the Coulomb corre- 
lations 

Hc=^^Y.^.C-^n,, (2.23) 
y 

where C~j^ is the capacitance matrix which can be found 
by solving a classical electrostatic problem of metallic 
particles embedded into the insulating matrix. Note that 
since metallic grains have infinite dielectric constant, the 
effective dielectric constant of the whole sample can be 
considerably larger than the dielectric constant of its 
insulating component. Thus the effective single grain 
charging energy can be much less than the electrostatic 
energy of a single grain in a vacuum. 

The Hamiltonian He describing the long range part 
of the Coulomb interaction, 12.23|l . has been derived 
for g ranular superconductors in the early work ( Efctov|. 
ll98Cl() . As far as the normal grains are concerned, the 
proper derivation is given in the review ijAleiner et al\. 
|2002() . The term H^, Eq. ifT^ . commutes with the free 
part -ffo.i) Eq. H2.16|l and therefore does not describe dy- 
namics in a single insulated granule. For the macroscopic 
transport to occur the electron trasfer from grain to grain 
has to be switched on, and we turn now to the third 
term Ht in Eq. (|2.15() describing the electron tunneling 
between the neighboring grains. We write this term in 
the form 

Ht= (2.24) 

where the summation is performed over the states p, 
q o f each grai n and o ver t he neighboring grains i and 
7. llAbrikosoyl Il988l: lAmbeeaokar and BaratoflL Il963t 
ICohen et oJl I1962D . It is assumed that the linear size 
of the inter-grain contact area well exceeds atomic dis- 
tances, i.e. there is a large number of conducting channels 
between the grains. At the same time, the magnitude of 
the potential barrier between the grains can be large, this 
can easily be achieved experimentally by the process of 
oxidation. Therefore, the tunneling conductance g, which 
is roughly proportional to the area of the contact and to 



the square of the tunneling matrix element tij , can be set 
both larger and smaller than unity. 

While the dependence of the tunneling elements Uj-^pq 
on the state numbers p, q is not important for transport 
through th e conventional point contacts and is usually 
neglected llAbrikosovl Il988t lAmbeeaokar and Baratofl 
Il963t ICohen et oJl I1962D . it may become relevant for 
granular arrays. Namely, the importance of this depen- 
dence is related to whether the electron motion in the 
grains is chaotic or can be viewed as an integrable bil- 
liard. In order to clarify this point we write the matrix 
elements tij-pq as 

tij\pq ~ j ^ij i^ij) 4'p (^ij) 4'q (^ij) dSij , (2.25) 

where is the coordinate on the junction between the 
grains i and j and (f>p (ji) are the wave functions. The 
integration over Sij in Eq. p.25|l extends over the area 
of the junction. 

If the grains do not contain any disorder and have 
a perfect shape (e.g. cubes or spheres), the array is a 
completely periodic system and according to the Bloch 
theorem the total resistivity is zero. The presence of 
the internal disorder or irregularities of the shape of the 
grains change the situation giving rise to a finite resis- 
tivity. In the latter c ase one may employ the random 
matrix theory (RMT) l|Alhassidl '2000: B eenakkeJ. Il997t 
IMehtal Il99l|) for an isolated single grai n or, which i s 
equivalent, the zer o dimensional a- model l|EfetovllT997(l . 
It is well known ("Mehtal, Il99l[l that this theory, when 
applied to the eigenfunctions leads to a Gaussian distri- 
bution W^{|(/)p (r)|} of their amphtudes 

l^{|(/.p(r)|} = exp(-|^p(r)|V) , (2.26) 

where V is the volume of the grain. Further, the Gaus- 
sian distribution of the amplitudes of the wave functions 
(j)p (r^) and the fact that their correlations decay rapidly 
with the increasing both spatial distance and the en- 
ergy levels spacing (the characteristic scales are the wave- 
length and the mean level spacing, respectively) leads to 
a Gaussian distribution of the matrix elements tij-pq. 
A simple transformation brings these correlators to: 

(^pigi ip2 92) ~ T; " (^PlP2 '^9l<32 + '^P1?2 '^P29l )i 
ZTT 

(tpq) = 0. (2.27) 

In Eq. H2.27|) the matrix elements tpq are taken for the 
same contact between the grains i and j, otherwise the 
correlations vanish. As we will see, the constant gij is 
nothing but the tunneling conductance of the contact. 

Equation H2.27|) is written for time reversal invari- 
ant systems (orthogonal ensemble). This, in particular, 
means that there is no magnetic field and/or and no mag- 
netic impurities present and therefore all eigenfunctions 
of an insulated grain can be chosen real. In the limit of 
comparatively strong magnetic fields, one arrives at the 
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a) b) 

FIG. 7 Self energy of the election Green function averaged 
over the impurity potential inside the grains and over tun- 
nelling elements between the grains. Averaging over the im- 
purity potential is represented by the dotted line (a) while 
tunnelling elements are represented by crossed circles (b). 



unitary ensemble and the first term in the correlation 
function of the matrix elements in Eq. H2.27|l vanishes. 
However, for the metallic grains of the size of the order 
of b-lOnm the characteristic magnetic field that require 
the description in terms of the unitary ensemble is of the 
order of several Tesla. Thus in the subsequent discussion 
we will be using Eq. (|2.27|l . assuming that the applied 
magnetic fields are not that high. 

To conclude this part, we have reformulated the ini- 
tial theory with impurities and regular tunneling matrix 
elements in terms of the model with random tunneling el- 
ements. This equivalence holds for the zero dimensional 
grains, or in other words, for the situations when all the 
characteristic energies are smaller than the Thouless en- 
ergy Exh, Eq- (|1.5I) . Working with the random tunneling 
elements turns out the more convenient approach. 

The model introduced in this subsection, Eqs. H2.15I 
ITTCl [TTTl VTM [T77|) describes a strongly corre- 

lated electronic disordered system which cannot be solved 
exactly. Below we introduce two complimentary ap- 
proaches to explore this problem in different regimes. 
One of these approaches - the diagrammatic technique - 
is especially useful provided the perturbative expansion 
with respect to one of the relevant parameters is possi- 
ble. The second method is based on the gauge transfor- 
mation which allows to get rid of the explicit Coulomb 
interaction term in the Hamiltonian at the expense of 
an appearing phase field. As we will show below, this 
approach offers a very powerful tool in its domain of ap- 
plicability. 



2. Diagrammatic technique for granular metals 

There are two routes to construct diagrammatic tech- 
nique for granular metallic systems. First, one can work 
in the basis of the exact single grain eigenfunctions and 
use the distribution of the tunneling elements H2.27|l for 
the description of the scattering between the different 
states of the neighboring grains. In this case the impuri- 
ties inside each grain are treated "exactly," and, by def- 
inition, the diagrams that represent impurity scattering 
within each grain do not appear in this representation. 
The alternative approach is, in essence, equivalent to the 
conventional cross technique: one begins with the mo- 
mentum representation and then carries out the standard 
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FIG. 8 Diagrams represent the Dyson equation for a sin- 
gle grain diffusion propagator Eq. 12.201 (a) and for the whole 
granular system Ea. 12.291 (b). Dotted lines represent the im- 
purity scattering while crossed circles stay for intergranular 
tunnelling elements. 

averaging over impurities within the each grain. The tun- 
neling between the grains can be then viewed as an inter- 
granular scattering. The corresponding matrix elements, 
as in the first method can, be viewed as Gaussian ran- 
dom variables. We prefer to follow the latter approach 
since it is straightforwardly related to the standard dia- 
grammatic technique for homogeneously disordered met- 
als making it easier to make comparisons to well known 
situations in disordered metals when possible. 

Following the guidelines of Abrikosov et al. (jlOBfj) we 
construct the action expansions with respect to both 
types of disorder: potential disorder within the each grain 
u (r) and the intergranular scattering matrix elements 
tpq. Shown in Fig. are the lowest order diagrams repre- 
senting both contributions to the self energy of the single 
electron Green function in the granular metals obtained 
by averaging over u (r) and tpq . The diagram a) describes 
the potential scattering within a single grain, while the 
diagram b) is due to the intergranular scattering. Both 
processes result in a similar contribution proportional to 
sgncj to the electron self-energy. This shows, that on the 
level of the single electron Green function; the intergran- 
ular scattering results merely in the renormalization of 
the relaxation time r 

=T-^ + 2rd, (2.28) 
where tq is the electron mean free time in a single grain. 

The next step is to consider the diffusion motion of 
electron through the granular metal. The diffusion mo- 
tion inside a single grain is given by the usual ladder 
diagram that results in the diffusion propagator D(uj), 
Eq. H2.2()|l . The tunneling between the grains does not 
change the selection rules for the diagrams. Typical di- 
agrams are shown in Fig. El b). The diagrams are gener- 
ated by connection the tunneling vertices, and only the 
diagrams without the intersection are to be kept. This is 
similar to what one has in the standard impurities tech- 



14 



nique ijAbrikosov et all Il965j) . In order to derive the 
expression for the total diffusion propagator one should 
sum up the ladder diagrams shown in Fig.|Hlb). For the 
periodic array of the grains we arrive at the following 
formula for the diffusion propagator D{uj,q) 



where 



DiLO,q)=T-'{\Lu\+TX^)-\ 



2 \^(1 — cos qa), 



(2.29) 



(2.30) 



a is the (super) lattice vectors, and q is the quasi- 
momcntum. In Eq. 12.29|l only the zero space harmonics 
of the diffusion motion in the grain is taken into account. 
This approximation is valid in the limit 



Th 



(2.31) 



where the energies F and Exh are given by Eqs. (|2.2|) and 
(|1.5|) . respectively. In the limit of small quasimomenta, 
q ^ a~^, we have a^q^, such that the propaga- 

tor (|2.29|) describes the diffusion motion on the scales 
much larger than the size of a single grain a with an 
effective diffusion coefficient 



(2.32) 



The other necessary blocks of the diagrammatic tech- 
nique can be constructed analogously to the diffusion 
propagator. 



3. Coulomb interaction and gauge transformation 

Taking an alternative route, one can eliminate the 
explicit Coulomb term in the model for the granular 
metals, Eqs. [TTEl ITT71 TTM [TTTji . by in- 

troducing a new phase field via a gauge transforma- 
tion. This approach proves to be especially useful in the 
regime of strong Coulomb correlations, where the dia- 
grammatic technique developed in the previous subsec- 
tion may break down. 

To introduce this approach it is convenient to adapt 
the formalism of the functional integration that al- 
lows replacing the calculations with the Hamiltonian H , 
Eq. (|2.15|l . by computation of a functional integral over 
classical fermion fields V'i (^) and their complex conju- 
gate tp* (X), where X — (r,T). These fields must satisfy 
the fermionic antipcriodicity condition 



V'* {r) = (r + /3) , 



(2.33) 



where (3 — \/T \s the inverse temperature. We define the 
field V'i {X) such that it is different from zero only in the 
ith grain. 

The Lagrangian L entering the functional integral 
for the partition function Z 



cxp 



(2.34) 



takes the form 

L[^P]=Y,Lo^['4'] + Lt[1p]+L4^P] , (2.35) 

i 

where 

(2.36) 

Ltm= E t^^r,P<lrp^{r)i^,J{T), (2.37) 



and 



Lcm^^T.''^ir)C-,'n,iT) , (2.38) 



where ^^(t) = J ip* {X)ipi{X)dr is the density field of the 
grain i. 

The Coulomb interaction Lc [t/j] , Eq. (|2.38() , can be de- 
coupled usi ng a Gaussian int egration over the auxiliary 
field V, (r) ijEfetov and Tschersich. .2003,) 



y exp f nT)i,+ iX)MX)dX 



:Z^}cxp [-S{V)]dV„ 



(2.39) 



where Vi{T) is the bosonic field obeying the boundary 
condition Vi (r) — Vi {t + (3) and Zy is the partition 
function, 

Zy - J exp [-S (V)] d%. 
The action S (y) has the following form 



We see from Eqs. (|2.32l - I2.39|l that the Lagrangian 
becomes quadratic in fields ip after the decoupling, 
Eq. H2.39|l . Instead of dealing with the Coulomb inter- 
action, £c[V'], Eq- (|2.38|l . one should consider now an 



effective Lagrangian 

[^P^y] for the gram i: 

[V', V] = Lo. m -I f% (r) ^+ (r) V-, (r) dr. 
Jo 

(2.41) 

The effective action L^f^ ^ip, V~\ is now expressed in terms 
of electron motion in a granular matter in the presence 
of the fluctuating potential Vi (r) of the grains. 
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We remove the field T^i(r) from the Lagrangian 
L^(^ Eq. (|2.41l) . using a gauge transformation of 

the fermionic fields 

(r, r) ^ e-^-^'^^V* (r, r) , ^, (r) = V, (r) , (2.42) 

where the phases ipi (r) depend on the imaginary time 
r but not on the coordinates inside the grains. This is 
a consequence of the form of the Coulomb interaction, 
Eqs. H2.23I I2.38|) . Since the action of an isolated grain 
is gauge invariant, the phases (fiir) enter the whole La- 
grangian of the system only through the tunneling matrix 
elements 



Vi{T), and the effective Lagrangian assumes the form 



(2.43) 



where ipij (r) — ipi (r) — ipj (r) is the phase difference of 
the i— th and j— th grains. 

At the first glance, the transformation, Eq. (|2.42(l . have 
removed completely the potentials V (r) from the effec- 
tive Lagrangian Ll{^ [i'-.V], Eq. In fact, this 
is not the case since the gauge transformation as de- 
fined in Eq. (|2.42|) , violates the antiperiodicity condition, 
Eq. (|2.33|) and the resulting effective Lagrangian for the 
single grain changes. In order to preserve the bound- 
ary condition, Eq. H2.33|l . the certain constraints should 
be imposed on the phases (r) and potentials Vi (r). 
Namely, the phases should obey the condition 



tpi (r) ^ ipi{T + (3) + 27rfci 



(2.44) 



which leads to the following constraint for the potential 
V^{r) 



V^{T)dT = 2TTh, 



(2.45) 



where fci = 0,±1,±2,±3.... 

The constraint on the phases (fi (r), Eq. I|2.44|l . can 
be reformulated via introducing a function (pi (r) assum- 
ing arbitrary real values from — oo to oo and obeying a 
periodicity condition 



(2.46) 



With the aid of this function the phase (r) can be 
written in the form: 



(pi (r) = 4).i (r) + 27rTkiT, 



(2.47) 



which satisfies Eq. H2.44|l . The potential Vi (r), in its 
turn, is taken in a form 

ViT)=p, + V,{T), (2.48) 

where the static variable pi varies in the interval 

- ttT < < ttT, (2.49) 

and the dynamic variable 1^(t) satisfies the constraint 
(|2.45|l . Note that the static part of the potential pi can- 
not be gauged out, contrary to dynamic contribution. 



^of^ [V', p]- In the limit of not very low temperatures. 



(2.50) 



the static potential pi drops out from the fermionic Green 
functions and does not influence the system behavior. 
This can be understood by noticing that at moderately 
high temperatures l|2.50|l the discreteness of the levels in 
the grains is not important. Then, in the Green func- 
tions, one may replace the variable £p — /i — s- ^, where Sp 
are eigenenergies and ^ is a continuous variable varying 
from — oo to oo. Integrating over ^ one may shift the con- 
tour of the integration into the complex plane and remove 
Pi provided it obe ys th e inequality (12.491) . More details 
can be found iu: Efetov and Tschersich ( 2003). Note that 
the variable pi cannot be neglected in the Lagrangian 
L^{^ [ip, p] in the limit of very low temperatures. 

An integer ki in Eqs. p. 441 12.451 12.471) represents an 
extra degree of freedom related to the charge quantiza- 
tion and is usually called the "winding number" . The 
physical meaning of the winding numbers becomes espe- 
cially clear in the insulating regime where they simply 
represent static classical electron charges. 

Thus, at not very low temperatures, H2.50|l . one can 
write the partition function Z in the form 



Z = j exp 
where 



- / (io M + Li [iP, Lp] + L2 [ip])dT 



D'ipD(p, 
(2.51) 



with Loi [-0] from Eq 
Li [ip, if] is 



Lo[V'] =^io.M, (2.52) 
The tunneling term 



LibP, 'P]^y2 d.TUj.pgip*p{T)tljjq{T) exp [iip.j (r)] , 
Jq 

(2.53) 

and the term L2 [(p] describing the charging effects reads 



Cij difiir) dipjir) 



2e2 dr 



dT 



(2.54) 



One should integrate in Eq. 1)2. 5 1() over the anticom- 
muting variables ip with the antiperiodicity condition 
Eq. H2.33|l . The integration over (f includes the inte- 
gration over the variable (pi (see Eq. (I2.47() ^ and summa- 
tion over ki. Equations (|T5]T - 1^511 ITI^ITTt)! completely 
specify the model that will be studied in the subsequent 
sections. 



4. Ambegaokar-Eckern-Schon functional 

The partition function Z, Eqs. (|2.51I - [T^ . can be fur- 
ther simplified by integration over the fermion fields -0. 
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This can hardly be done exactly but one can use a cumu- 
lant expansion in the tunneling term Li[ip, ip], Eq. H2.53|l . 
Of course, even in the absence of the tunneling the ran- 
dom potential in Lo Eqs. (jTCTIT^ . can make the 
problem highly non-trivial. Yet in the limit of not very 
low temperatures, Eq. (|2.5Q(I . the problem can be simpli- 
fied because of the possibility of neglecting interference 
effects and considering the disorder within the SCBA tak- 
ing Green functions from Eq. H2.18|l . We perform the 
cumulant expansion in the tunneling term 5*1 ["0, ip] up to 
the lowest non- vanishing second order, integrating over 
fermionic degrees of freedom and averaging over the tun- 
neling matrix elements with the help of Eq. (|2.27|l . 

The resulting action, originally derived 
in lAmbegaokar et al\ (1982) for a system of two 
weakly coupling superconductors, contains only the 
phases ipi (r) but no fermionic degrees of freedom left. 
Then for the granular metals the partition function Z 
can be written as 



Z 



exp (—5) Dip, S = Sc + St 



(2.55) 



where Sc describes the charging energy 



and St stands for tunneling between the grains 



St — Trg'S^ / drdr'a (r — t') sii 



(2.57) 



The function a (r — r') in Eq. H2.57|l has the form 

a (r - r') = T^Resin"^ j^yj-^ - t' + iij)] . (2.58) 

where rj +0 and one should take the real part in 
Eq. 

Despite the fact that the above functional was obtained 
via an expansion in the tunneling term Si[ip,ip], its va- 
lidity is not limited by the insulating regime, 5 <C 1, only. 
The functional can be used in the metallic regime at tem- 
peratures T ^ r, where F is given by Eq. (|2.2() . 

Of course, the phase action looses some information 
about the original model, Eqs. H2.34I - IT^ . and its ap- 
plicability in each particular case has to be carefully ana- 
lyzed. In general, the functional S, Eqs. (|2.55l - [T55|l does 
not apply in cases where the coherent diffusive motion of 
an electron on the scale of many grains is important. For 
example, the weak localization correction cannot be ob- 
tained within this approach. 

Moreover, one has to be careful in using this approach 
even in the weak coupling regime, 5 <C 1. For example, 
the hopping conductivity in the low temperature elastic 
regime is also beyond the accuracy of the phase action, 
Eqs. H2.55I - since it requires consideration of the 

elastic multiple co-tunneling processes which the phase 
action misses. 



Yet the phase functional approach is extremely pow- 
erful for many applications. For example, it enables ob- 
taining non-perturbative results for the conductivity in 
the metallic regime at temperatures T 3> F. 



C. Metallic properties of granular arrays at not very 
temperatures 



low 



In this section, we will derive the logarithmic correc- 
tion to the conductivity, Eq. (|2.6|l . in the temperature 
regime T 3> F. While this correction has the same origin 
as the Althsuler-Aronov (AA) re sult derived for the ho- 
mogeneously disordered metals |Altshuler and Aronovl 
Il985t II ee and Ramakrishnanl Il985|) . its form is specific 
to granular metals because the AA correction. In con- 
trast to Eq. (|2.6|l . AA corrections are sensitive to the 
sample dimensionally and exhibit the logarithmic behav- 
ior only in two dimensions. 

The logarithmic correction, Eq. (|2.6|l . as we will show 
below, can be obtained as a result of the renormalization 
of the tunneling coupling g due to the Coulomb corre- 
lations. It comes from the short distances, and this ex- 
plains its insensitivity to the sample dimensionally. At 
temperatures smaller than F the coherent electron mo- 
tion on the scales of many grains becomes important, and 
the conductivity correction acquires the form similar to 
that of the homogeneously disordered metals. 

For the sake of simplicity we will consider a periodic 
system assuming that the granular array forms a cubic 
lattice. This implies that all the grain sizes and inter- 
granular conductances are equal to each other. At the 
end of this section we will discuss the applicability of the 
result (|2.6|l to systems containing irregularities. We will 
start with perturbation theory in 1/g, then generalize our 
derivation of Eq. H2.6|l using the renormalization group 
approach. 



1. Perturbation theory 

In the limit of large tunneling conductances 5^1 the 
tunneling term St, Eq. (|2.57|l . suppresses large fluctua- 
tions of the phase </3(r). In this regime the tunneling 
term St can be expanded in phases fijir) because they 
are small. The charge quantization effects in this regime 
are not pronounced allowing to neglect all non-zero wind- 
ing numbers fc^. At the same time, the phase fluctuations 
can change considerably the classical result, Eq. (|2.1|) . 

The quadratic in (pi^r) part of the action S in Eq. (|2.55l 
- 12.5711 will serve as the bare action for the perturbation 
theory. Keeping terms of the second order in (pi^r) in 
Eqs. I|2.55l - 1^.571) and performing the Fourier transfor- 
mation in both the coordinates of the grains and the 
imaginary time we reduce the action S in Eq. (|2.55l) to 
the form 



S'o = T ^ ipq^n 
q,n 



(2.59) 
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FIG. 9 This diagram represents the conductivity in the lead- 
ing order in 1/g. The crossed circles represent the tunnelling 
matrix elements e"^'^ where phase factors appear from 
the gauge transformation. 



where 



G„^^u;l/4E{q) + 2g\uJn\X^. 



(2.60) 



where Aq is given by Eq. (|2.3Q(I . Here q are the 
quasi-momenta of a periodic array, £'(q) is the Fourier 
transform of the charging energy related to the Fourier 
transform of the capacitance matrix C (q) as E{q) — 
e^/2C(q). The conductivity cr (w) of the array of the 
grains can be obtained from the current-current corre- 
lation function via the standard Kubo formula. The con- 
ductivity in the leading order in 1 /qis shown in Fig. El Its 
analytical expression is given by ijEfetov and Tschersichi 

l2noa 



a{uj) = 



2'Ke gT ia 



2~d 



dr 



1 



■ exp 



(-Ga(r)) 



to Jq sin^(7rTT) 

(2.61) 

where the analytic continuation to the real frequency is 
assumed as £7 ^ —itu and the function Ga (t) is 



Ga (r) = ATd 



UJ„>0 



d q „ . 2 qa . 9 lli„ t 

jGan sm — sm" 

{27rf ^ 2 



^2 



(2.62) 

where d is the dimensionality of the array. The factor 
exp ^— Ga (t)^ in Eq. (|2.61(l is responsible for the inter- 
action effects and appears as a result of the averaging 
of the phase exponents e"^'^'^^ emerging after the gauge 
transformation, Eq. H2.42|l . One can see from Eqs. I|2.6()l 
12.621 ) that the function Ga (t) behaves as a logarithm 
\n{gEcT) at the values of r of the order 1/T that are 
essential for calculation. With the logarithmic accuracy 
one can neglect the term in G~^, Eq. I|2.6U|1 . and 
reduce Eq. to 



rp •^c 

do — ' cj. 



1 — cos (WjiT) 



(2.63) 



In Eq. H2.63|l one should sum over positive Matsubara fre- 
quencies up to the cutoff ujc ^ gEc- For larger frequencies 
the first term in Eq. 1)2. 6U|) is no longer small and there is 
no logarithmic contribution from these large frequencies. 
Equation (|2.63() shows a remarkable independence of the 
result on the structure of the lattice. What is also impor- 
tant, there are no "infrared" divergencies in the integral 
over q in any dimensionality including 2c? and Id. 



With the logarithmic accuracy, one can replace r by 
1/T in the function Ga(T) and calculate the remaining 
integral over r in Eq. ()2.61(l ignoring the dependence of 
the function Ga on r. Taking the limit a; — > we obtain 
the result (|2.6() in the temperature interval T 3> F : 



cr = (To 1 - 



1 



2Tidg 



In 



T 



(2.64) 



It was shown in lEfetov and TschersichI l)2003j) that the 
terms of the order (l/5)^ln^ (gEc/T) are cancelled out 
in the expansion of the conductivity correction, which 
means that the accuracy of Eq. (|2.64|l exceeds the ac- 
curacy of the first order correction. Furthermore, this 
cancellation is not accidental and the result H2.64I) turns 
out to be applicable even at temperatures at which the 
conductivity correction becomes of the same order as ctq 
itself. This fact is explained in the next subsection where 
the temperature dependence of the conductivity is con- 
sidered within the RG scheme. 



2. Renormalization Group 

In order to sum up all the logarithmic corrections 
to the conductivity we use RG scheme suggested for a 
one-dimensional XF-model long ago fiKosterlitz, ll97d) 
and used later in a number of works llBulgadaevL Il984t 
iFalci et oi.LllQQ.'iHGuinea and Sch5iiLll986^ . As the start- 
ing functional we take the tunneling action St 



St = Trg 



drdr'a (r — t' 



(P2j (r) - ip,j (r')" 



(2.65) 

The charging part Sc is not important for the renormal- 
ization group because it determines only the upper cutoff 
of integrations over frequencies. In the limit T the 
function a{T — r') is proportional to (r — r') ^ and the 
action is dimensionless. We neglect in this subsection 
the non-zero winding numbers ki and replace the phases 
ipi (t) by the variables (j>i (r), Eq. H2.47|l . 

Following standard RG arguments we want to find how 
the form of the action St changes under changing the cut- 
off. Generally speaking, it is not guaranteed that after 
integrating over the phases 4> in an interval of the fre- 
quencies one comes to the same function sin^ in the 
action. The form of the functional may change, which 
would lead to a functional renormalization group. 

In the present case appearance of terms sin^ 20, 
sin^ etc. is not excluded and, indeed, they are gen- 
erated in many loop approximations of the RG. Fortu- 
nately, the one loop approximation is simpler and the 
renormalization in this order results in a change of the 
effective coupling constant g only. 

To derive the RG equation we represent the phase </> 
in the form 



(2.66) 
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where the function 



is the slow variable and it is 



not equal to zero in an interval of the frequencies < 



uj < Xujr, while the function 



is finite in the interval 



Xujc < Lu < ujc, where A is in the interval < A < 1. 
Substituting Eq. lfTHH|l into Eq. l(TH5|l we expand the 
action St up to terms quadratic in Integrating in 

the expression for the partition function 



Z = Jcxp i-StmD(^ 



over the fast variable (jj^^j^ with the logarithmic accuracy 
we come to a new renormalized effective action 5*4 



2TTg 



P rl3 



drdr'a (r — t') 



(t) - (r') 



2iTgd 



,(2.67) 



where <^ = — In A. It follows from Eq. (|2.67|) that the form 
of the action is reproduced for any dimensionality d of the 
lattice of the grains. This allows us to write immediately 
the following renormalization group equation 



1 

2nd' 



(2.68) 



The solution of Eq. (|2.68|l is simple. Neglecting the 
Coulomb interaction in the action St, Eq. (|2.65|l . is justi- 
fied only for energies smaller than gEc and this gives the 
upper cutoff. Then, the renormalized conductance g (T) 
takes the form 



9{T)^g - — In—. 

27ra 1 



(2.69) 



and we come using Eq. lEH to Eq. for the conduc- 

tivity. Both the quantities depend on the temperature 
logarithmically. 

Equation H2.69|l is obtained in the one loop approxi- 
mation and should be valid so long as the effective con- 
ductance g (T) remains much larger than unity. This 
assumes that the result of the perturbation theory, 
Eq. (|2.64ll can be extended to a wider temperate interval 
and, in fact, it is valid as long as 



g-{2nd) hn{gE,/T) » 1. 



(2.70) 



At high temperatures, when the inequality (|2.7l)|) is ful- 
filled, one can speak of metallic behavior of the system. 
At lower temperatures T < Tc, where 



Tc = gEc exp {-2Trgd) , 



(2.71) 



the system is expected to show insulating properties. 

The result for the conductivity correction H2.64|l was 
obtained within a cubic lattice model. It is important 
to understand how it will change in a more realistic case 
of an irregular array. First, we note that Eq. (|2.64|) can 



easily be generalized to the case of an arbitrary periodic 
lattice with the result 



(T = do 1 In 

TTZg 



gEc 
T 



(2.72) 



where z is the coordination number of the arbitrary pe- 
riodic lattice. 

The role of t he dispersion of the tunneling cond uctance 
was studied in iFeigel'man and loselevichi l)2005(l within 
the 2d model that assumed regular periodic positions and 
equal sizes of grains but random tunneling conductances. 
It was shown that the dependence H2.72|l holds in a wide 
temperature range in the case of a moderately strong 
conductance dispersion. However, the distribution of the 
conductances broadens and this effect becomes especially 
important close to the metal-insulator transition where 
it was suggested to describe the transition in terms of 
percolation, (FciRcl'man and losclcvich, 2005). One may 
expect that in 3d samples the effect of the conductance 
dispersion on the macroscopic transport will be less im- 
portant than in 2d ones, while, on the contrary, it will be 
more dramatic in Id samples. 

Finally, we would like to comment on the validity 
of Eq. H2.72|l in the case of an arbitrary irregular ar- 
ray. First, we note that in the RG approach the ca- 
pacitance disorder leads to a local renormalization of 
the tunneling conductances. Then, the general irregu- 
lar system can be viewed as an irregular array of equal 
size grains with random tunneling conductances. The 
ma in difference of such a system fr om that considered 
in ijFeigel'man and loselevichi l2005(l is that the coordi- 
nation number of each grain varies from grain to grain. 
In such a case one expects that the coefficient in front 
of the second term in Eq. H2.72|l will be determined by 
an effective coordination number in a system that is ex- 
pected to be close to the average number of neighbors of 
a grain. 

Following the RG scheme we did not take into account 
non-zero winding numbers fc^, Eqs. 1)2.441 - [2.47|l . This 
is a natural approximation because the contribution of 
such configurations should be exponentially small in g 
(~ exp(— eg), where c is of the order of unity). So long 
as the effective conductance g is large in the process of the 
renormalization, the contribution of the non-zero winding 
numbers ki can be neglected. They become important 
when g becomes of the order of unity. This is the region 
of the transition into the insulating state and apparently 
the non-zero winding numbers play the crucial role in 
forming t his state. (However, the authors of recent pub- 
lications l)Altland et all l2004 l2005t iMever et al\ l2004j) 
argued that the non-zero winding numbers might become 
important in Id samples at temperatures T* parametri- 
cally exceeding T^, Eq. (|T7T|) '). 



(a) 



(b) 






FIG. 10 Diagrams representing the vertex correction (a) and 
the renormahzed Coulomb interaction (b). 



D. Metallic properties of granular arrays at low 
temperatures 

The phase functional technique may not be used for ob- 
taining the conductivity corrections at temperatures T < 
r since in this regime the coherent electron motion over a 
large number of grains , which is mi ssed in the functional 
approach (Ambcgaokar et all . ll982t) . becomes important. 
In this section we will derive the conductivity correction 
within the diagrammatic perturbation theory described 
in the Sec. lII.CTl Unlike the phase functional approach, 
this technique does not allow us to obtain the non pertur- 
bative results in an easy way but it has an advantage of 
being applicable at arbitrary temperatures. In particular, 
we will see below that it reproduces at T ^ F the pertur- 
bative result, Eq. H2.64|l . As we will show, the conductiv- 
ity correction agrees in the low temperature regime with 
the one obtained in Ref s. Ja tshuler and Aronov*. 'l985l 
feclitz and KirkDatrick, Il994 : lEee and R amakrishnan. 
[l985i|l for homogeneously disordered samples 

In Section FlFCll we have already described the main 
building blocks of the diagrammatic technique includ- 
ing, in particular, the diffusion propagator defined by 
the ladder diagram in Fig. |S1 The same ladder diagram 
describes the dressing of the interaction vertex as shown 
in Fig.llOb. The dressed vertex can be used to obtain the 
polarization operator that defines the effective dynami- 
cally screened Coulomb interaction fFig. IlUb'): 



y(f],q) 



C(q) 



2A„ 



AqF 



(2.73) 



The conductivity of the granular metals is given by the 
analytical continuation of the Matsubara current-current 
correlation function. In the absence of the electron- 
electron interaction the conductivity is represented by 
the diagram (a) in Fig, ill! that results in high tempera- 
ture (Drude) conductivity ctq defined in Eq. 1)2. First 
order interaction corrections to the conductivity are given 
by the diagrams (b-e) in Fig. Th ese diagrams are 
analogous to those considered in lAltshiiler and Aronovl 
(l985.1 for the correction to the conductivity of homoge- 
neous metals. 

We consider the contributions from diagrams (b,c) and 
(d,e) separately: The sum of the diagrams (b,c) results 



(a) 



(b) 



(c) 





FIG. 11 Diagrams describing the conductivity of granular 
metals: the diagram (a) corresponds to ctq in Eq. 12.511 and 
it is the analog of Drude conductivity. Diagrams (b)-(e) de- 
scribing first order correction to the conductivity of granular 
metals due to electron-electron interaction. The solid lines 
denote the propagator of electrons and the dashed lines de- 
scribe effective screened electron-electron propagator. The 
sum of the diagrams (b) and (c) results in the conductivity 
correction Sai in Eq. 12.511 . The other two diagrams, (d) and 
(e) result in the correction Sa2. 



in the following correction to the conductivity 

^Im^ J d^7HAq V{io,q), (2.74) 



5(Tl 



where 7(0;) = ^wcoth^, and the potential V{uj,c[) is 
the analytic continuation of the screened Coulomb poten- 
tial V{^1, q) with dressed interaction vertices including 
those attached at the both ends 



V{cu,q) 



ix^r-iu) (4Aq£;c(q) 



(2.75) 



The expression for the screened Coulomb interaction 
V{uj,q), Eq. (|2.75|l . is written in a simplified form us- 
ing the fact that the charging energy £'c(q) — e^/2C(q), 
expressed in terms of the Fourier transform of the capac- 
itance matrix C(q), is much larger than the escape rate 
F. 

Performing the integration over the frequency and 
summing over the quasimomentum q in Eq. I|2.74|l with 
the logarithmic accuracy we obtain the correction to the 
conductivity, Eq. H2.6|l . One can see from Eq. (|2.74() that 
the contribution Sai in Eq. (|2.6|) comes from the large 
energy scales, e > F. At low temperatures T < F, the 
logarithm is cut off by the energy F and is no longer 
temperature dependent. 

To obtain the total correction to the conductivity of 
a granular sample the two other diagrams, (d) and (e) 
in Fig. ^2 should be taken into account. These dia- 
grams res ult in the following cont ribution to the con- 
ductivity ijBeloborodov et ajl\200^ 



5(72 
CTO 



2g5_ 



E 

q 



duj 7(w) Im^ 



y(cc>,q)EaSin"(qa) 
AqF — iuj 



(2.76) 

In contrast to the contribution 6ai, Eq. H2.74|l . the main 
contribution to the sum over the quasimomentum q in 
Eq. (|2.7t)|) comes from the low momenta, q <C 1/a. In 
this regime the capacitance matrix C(q) in Eqs. (|2.75|) 
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and H2.76|l has the asymptotic form 



Hl/qa) 






■n/q 


d = 2, 


(2.77) 


2^/g2 


d = 3. 





Using Eqs. 12.75I - IT77I) . we obtain the result for the cor- 
rection (5(72 ill Eq. 1)2. 7|l . This correction has a physical 
mean ing similar to that of the Alt shuler-Aronov correc- 
tion ijAltshuler and Aronovl Il985j) derived for homoge- 
neously disordered metals. 

Comparing results in Eqs. (|2.5|l - (|2.7|l with those ob- 
tained in the previous subsection using the AES func- 
tional one can see that the correction to the conductivity 
obtained in Sec. III. d is equivalent to the correction Sai 
in Eq. (|2.5|) . which corresponds in the diagrammatic ap- 
proach to the sum of diagrams (b) and (c) in Fig. 
The correction Sa2 in Eq. H2.5fl becomes important only 
at low temperatures, T <T where AES functional is not 
applicable. 

Using the final results of the calculations, Eqs. (|2.6|) 
and H2.7|l , we can make a very important statement about 
the existence of a metal-insulator transition in three di- 
mensions. It follows from Eq. H2.7I) that for a 3d granular 
array there are no essential corrections to the conductiv- 
ity at low temperatures T <C T coming from the low 
energies uj <T, since the correction Sa2 is always small. 
This means that the result for the renormalized conduc- 
tance, Eq. H2.69|l . for 3d samples can be written within 
the logarithmic accuracy in the form 



~9iT) 



67r 



■In 



gEc 



max {gS, T) 



(2.78) 



such that it is valid for all temperatures as long as the 
renormalized conductance is large, g{T) 3> 1. 

One can see from Eq. H2.78|l that for a large bare 
conductance, g 3> (l/67r) ln((7i?c/<5), the renormalized 
conductance g is always large and the system remains 
metallic down to zero temperatures. In the opposite 
limit g < {l/6TT)lii{gEc/6), the system flows when de- 
creasing the temperature to the strong coupling regime, 
^ ~ 1, which indicates the onset of the insulating phase. 
One can see that with the logarithmic accuracy the criti- 
cal value of the conductance gc is given by Eq. H2.12|l . 
and t his value separates the metallic and insulating 
states l|Beloborodov et allEoO^ . 

One can see from Eq. H2.7|l that in Id and 2d sam- 
ples the correction to the conductivity, being negative, 
grows with decreasing the temperature and diverges in 
the limit T 0. Such a behavior is usually attributed 
to a localization in the limit T — s- 0. 

However, the situation is more interesting. In a re- 
cent paper fBas ko et a/., 2 005) , it was demonstrated that 
there must be a metal-insulator transition at finite tem- 
perature Tk in systems with a weak repulsion provided 
all one-particle states are localized. This means that 
the conductivity is strictly zero at T < Tk and be- 
comes finite at T > Tk- The Id and 2d granular systems 



with the Coulomb interaction discussed in the present 
review should belong to t he class of the models consid- 
ered in ijBasko et al\ . \200^ (all one-particle states should 
be localized in Id and 2d systems for any disorder) and 
one can expect such a transition. Clearly, the results of 
Eq. (ETjl should hold for T > Tfc. 

At the same time, the 3d samples do not f all into the 
class of the models studied by Basko et all |2p05) be- 
cause there exist both localized and extended one-particle 
states. Therefore, for 3d granular systems, one can speak 
of the metal-insulator transition at T = only. The tran- 
sition point can be varied by changing, e.g. the tunneling 
conductance between the grains or their size. 

Finally, we would like to note that the low temperature 
conductivity correction da2 is less sensitive to a partic- 
ular model of a granular sample than the high temper- 
ature result, Eq. H2.64|l . As we will discus in the next 
section, the low temperature conductivity correction is 
determined by the scales larger than the single grain size 
and can be expressed in terms of the effective diffusion 
coefficient D^f / . This assumes that the underlying struc- 
ture of the model is not important in this regime. 



E. Universal description of granular materials 

One can see from the previous subsection that at low 
temperatures, T <^ T, the dependence of the correction 
(5ct2, Eq. (|2.7|) . to the conductivity of granular metals 
coincides exactly with the corresponding result for the 
conductivity of the homogeneously disordered samples. 
A question immediately arises: is it an accidental coin- 
cidence that the two different physical systems exhibit 
in the main approximation the identical low tempera- 
ture transport behavior, or there is an underlying deep 
connection between the two? Can on describe both the 
region T < T and T > F in an unified manner? The 
main goal of this subsection is to answer these important 
questions. 

A very convenient description of the low tem- 
perature behavior of homogeneously disordered met- 
als is based on the so-called (T-model. There are 
several formulatio ns of t he cr-model based on the 
replica trick (lEfetove^/l 
persvmmet rv llEfetovl 11983. 



19801: IWegneii [Hzi, 



su- 



1997ft . and K eldysh Green 
functions (Kamen ev and Andreevl Il999j) . The ap- 
pr oach of ..Efetov et al\ lll980D has been generalized 
by iFinkelsteinTi^gol) to inclu de the interaction effects. 
It is clear that the cr-model of lFinkelsteinI l)l990f) cannot 
be used for the granular materials because it does not 

contain the scale F. 

Fortunately, the approach of IFinkelsteinI |)1990|) can 
rather easily extended to include the granularity. This 
can be done using Eqs. I|2.15ll - 112.2311 and f ollow- 
ing the usual way llAndreev and Beloborodovl l2004t 
iBeloborodov et a/.ll20o'l|) : 1) we write a generating func- 
tional in terms of functional integrals over anticommut- 
ing Grassmann variables, 2) average over the disorder 
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using a replica trick, 3) we decouple a -0^ effective in- 
teraction appearing after the averaging by a Gaussian 
integration over a Q- matrix field and do the same for 
the Coulomb interacti on term in Eq. (12.2311 using axil- 
lary fields V (Efetov eFoH . Il980t iFinkclstein. .199Q) , 4) 
we integrate thus obtained exponential of a quadratic 
form over the "0 fields and finally, 5) find a saddle point 
in the free energy functional containing the fields Q and 
V only and expand around this saddle point. The final 
expression for the effective low-energy action reads: 



(2.79) 



Here the sums are performed over the grain indices, the 
symbol (...) means summation over the nearest neigh- 
bors, Cij is the capacitance matrix, e = idr and the 
symbol "Tr" means the trace over spin and replica in- 
dices and integration over r. The field V in Eq. (|2.79|) 
is a time dependent vector in the replica space and the 
corresponding "scalar product" is implied: {ViVj) — 

lo^^ (■'') ^ja {t) dr , whcrc a is the replica in- 
dex (0 < a < TV). The Q- matrix in Eq. l(T7^ is 
the matrix in the time (it depends on two times r and 
r'), spin, and replica spaces subject to the constraints 

— 1, TrQ = 0. The variable r enters Q in the 
same way as the replica and spin ones, which means that 

iQlQ2)r,r' =lQl (r, r") Q2 (t", t') dr" . 

For energies smaller than the Thouless energy Erh of 
one grain, Eq. ( 11.51) . the Q-matrices in Eq. H2.79|l are 
coordinate independent within each grain. They can be 
written in the form 



where the function A equals 



Ar,r' = 



iT 



sinTrT (r — r') ' 



(2.80) 



(2.81) 



and Ui is an unitary matrix. The action >5'[g,y] in 
Eq. H2.79|l describes the entire region of both low, T < F, 
and not very low, T > F, temperatures discussed in the 
previous subsections. For T 3> F, an essential contribu- 
tion comes from U (r, r') in Eq. (|2.80|l diagonal in time, 
spin and replica spaces. This means that we write this 
function as 



(r, t') = Sa/iS (r - r') exp [iipat (r)] , 



(2.82) 



where a, (3 stand for both replica and spin indices. 

Substituting Eqs. (|2.8UI2.82|I into Eq. l(?7^ . calculat- 
ing the integral 



J exp {-S[Q,V])DV 



and taking the limit N = 0, which is trivial in this limit, 
we come to the AES action, Eqs. I|2.55l I2.5S|| . 

The (7- model, Eq. H2.79|l , has a well-defined continuum 
limit. It can be obtained for slow spatial variations of the 
Q-niatrix by expanding the second term in Eq. H2.79|l to 
the second order in gradients of Q, and replacing the 
summation over i by the integration over r. The tunnel- 
ing term in Eq. (|2.79|l describes in the continuum limit 
the diffusion and we reduce the action, Eq. (|2.79() . to the 
form 



Tr 



{e + V)Q-^{\7Q)' 



dr 



drdr' , 

n2d 



-Tr 



2e2 



(2.83) 



where ly is the density of states. The coefficient D in the 
second term of Eq. (|2.83() is in this approximation just 
the classical diffusion coefficient determined by Eq. H2.3|l . 

However, we can derive Eq. (|2.83(l more accurately, not 
just neglecting the contribution of energies exceeding F 
but taking them into account. Calculating the contribu- 
tion of the energies exceeding F we can use Eq. (|2.82|l for 
the excitations with such energies. As a result, we come 
again to Eq. (|2.83|) but with a renormalized coefficient D 
that can be written in the form 



where geff equals (c.f. with Eq. 



9efS ^ 9 



— In^ 



(2.84) 



(2.85) 



Since the effective model H2.83|l operates with matrix 
fields Q that have only long range degrees of freedom, it 
applies after an appropriate high energy renormalization 
to any disordered metal including a homogeneously dis- 
ordered one. Integrati ng over the potent ials V one can 
come to the cr-model of lFinkelsteinI l)l99Cf) (the limit of a 
long range interaction is implied). 

Thus, all the information about the granularity of the 
sample is hidden in the low temperature limit in the coef- 
ficients of the effective model ()2.83(l . The conductivity of 
the sample is related to the effective diffusion coefficient 
through the usual Einstein relation 



2e^D{a'^5)- 



(2.86) 



The effective model, Eq. (|2.83|l . together with equa- 
tion H2.69|l for the renormalized conductance naturally 
explains the results for the low temperature, T ^ F, 
conductivity described in section III.DI The interaction 
correction to conductivity consists of two terms: 1) The 
first contribution is temperature independent and comes 
from the high energy renormalization of the diffusion co- 
efficient D, Eqs. (jTSBl EHUl- It can be written in the 
form (c.f. with Eq. 



5ai = -CTo 



1 

2'Kdg 



In 
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5 



(2.87) 
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where ctq is the classical Drude conductivity defined in 
Eq. (|2.1|l . This contribution is specific for the granular 
materials. 

2) The second contribution to conductivity, Eq. I|2.7|l . 
is temperature dependent and comes from the low energy 
renormalization of the diffusion coefficient D in the ef- 
fective model (|2.83|> . It coincides with the corresponding 
correction to cond uctivity obtained for homog eneously 
disordered metals l|Altshuler and Aronovl llQSSji . 

Making use of the effective description of the granu- 
lar metals in terms of effective the cr- model, Eq. H2.83|l . 
applicable at low temperatures we conclude: all the phe- 
nomena that are described in terms this model including 
localization effects, magnetoresistance, Hall conductivity, 
etc., are universal. 

Using the cr-model, Eq. H2.88|l . we also describe quan- 
tum interfere nce (weak localizatio n) corrections to the 
conductivity, l)Gorkov et al\ . Il979|) . In the leading order 
in the inverse tunneling conductance, 1/g, the interaction 
and weak localization corrections can be considered in- 
dependently. The quantum interference corrections may 
be obtained from the effective fi-model, Eq. (|2.83|l . us- 
ing directly the cor responding results f or homogeneously 
disordered metals l)Efetov et all \l9H(t^ with the proper 
effective diffusion coefficient, D = ga^S. For 2D and 
ID samples it is important to take into account de- 
phasing effects since the weak localization correction di- 
verges. The dephasing time can also be extracted 
directly from the results for homogeneously disordered 
metals (^Altshuler et all Il982(l with the proper effective 
diffusion coefficient, D. The final result for the weak 
localization corrections is given by Eqs. and (jTi^ . 



F. Insulating Properties of Granular Metals: Periodic model 

In the previous section we discussed the transport 
properties of metallic granular arrays with strong inter- 
granular tunneling coupling. Now we turn to the oppo- 
site limit of weakly coupled grains, g <C 1. We start our 
consideration with a detailed description of a periodic 
granular array model that is presented in this section. 
To remind, the periodic model assumes periodic arrange- 
ments of equal size grains as well as the absence of the 
disorder in conductances and grain potentials. At the 
same time, the electron motion inside the grain is chaotic 
and this brings a disorder into the model. 

As we discussed in the Introduction, the periodic gran- 
ular array model predicts the insulating conductivity be- 
havior with a hard gap in the electron excitation spec- 
trum and, for this reason, it cannot describe the trans- 
port properties of realistic granular systems where the 
conductivity governed by the variable range hopping 
mechanism. Nevertheless, it is instructive to consider 
this model for two reasons: first, it is important to be see 
explicitly that such a model may indeed produce only the 
activation conductivity dependence; second, this model 
illustrates a general effect of the charge renormalization 



of the Coulomb energy due to the intergranular coupling. 

1. Activation conductivity behavior 

First, we will consider a periodic granular array with a 
small intergranular coupling, g <^ 1. The conductivity in 
this case can be found via the gauge transformation tech- 
nique described in Sec. III. B. 31 The phase action in the 
lowest order in tunneling conductance g is simply given 
by the Coulomb term Sc, Eq. With the help of 

the Kubo formula the conductivity in the first n onvanish- 
ing o rder in g may be written as ([Efetov and TschersichI 

= y / dr . 2, %. , n(T). (2.88) 

w Jo sm (ttJ t) 

Here the analytic continuation to the real frequency is 
assumed as i7 — s- —iuj and the function II(t) represents 
the phase correlation function 

n(ri - rs) - (exp (i {ip (n) - ^ (ts)))), (2.89) 

where the averaging is performed with the Coulomb ac- 
tion Sc in Eq. 

At first glance, calculation of the correlation func- 
tion n (ri — T2) with the action Sc reduces to compu- 
tation of a Gaussian integral. However, it is not so 
because at finite temperatures it is necessary to take 
into account all winding numbers, Eq. (|2.47|l . in order 
to obtain the correlation function H(r) correctly. At 
lo w temperature the straightfor ward calculation results 
in ^Efetov and TschersichI l2nn3r 

cr = 2cto exp i-Ec/T) , (2.90) 

where ctq is the Drude high temperature conductivity and 
Ec is the single grain charging energy. This result has a 
clear physical meaning: The conductivity of the insu- 
lating granular array is mediated by electrons and holes 
that are present in the system as real excitations. Their 
concentration is given by the Gibbs law and the factor 2 
appears due to the presence of both electrons and holes. 
From the result H2.9()|l we conclude that the activation 
exponent is determined by the Mott gap - that is the 
energy one pays to add or remove an electron from the 
system. 

One can extend the calculation of the conductivity to 
higher orders in g expanding in the tunneling term St, 
Eq. lE33l- It was shown bv iLoh et all l)2005ji that the 
activation behavior persists at least in the next order in g. 
Apparently, the activation behavior, Eq. H2.13|l is quite 
universal for the periodic arrays and what one should 
clarify is the dependence of the Mott gap Am on the 
conductance g. Now we turn to this question. 

2. IVIott gap at small tunneling conductances 

At finite intergranular coupling the Mott gap turns out 
to be reduced due to the processes of virtual tunneling 
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of electrons to the neighboring grains. This effect can be 
most easily studied in the limit of the low coupling be- 
tween the grains within the straightforward perturbation 
theory in the tunneling conductance. The gap Am can 
be formally defined as 



E 



N=l 



E 



(2.91) 



where ji is the chemical potential, E^^q is the ground 
state energy of the charge neutral array and En=i is 
the minimal energy of the system with an extra electron 
added to the neutral state. The correction to the Mott 
gap due the to the electron tunneling to the neighboring 
grains can be found in the second order of the perturba- 
tion theory 



ae^-j:- 



Ek ~ Eq' 



(2.92) 



where the matrix elements of a perturbation V are taken 
between the ground state and excited states k. 

The correction to the energy due to the finite inter- 
granular coupling should be included in both the terms 
En=i and En=q in Eq. (|2.91(l . We consider in the zero 
approximation isolated grains and the tunneling Hamil- 
tonian Ht, Eq. (|2.24|) . is our small perturbation. Using 
Eq. H2.92|l we have to calculate matrix elements of the 
tunneling between neigboring grains. As all grains are 
equivalent, we consider the tunneling between the grains 
1 and 2 assuming that the charge N on the grain 1 can be 
or 1, while the grain 2 is initially neutral. Contribution 
of the hops between the grain 1 and all its other neigh- 
bors leads merely to the factor z (coordination number) 
in the final result. Using Eq. H2.92|l we should calculate 
the corrections to the energies En^q and E^^i. 

First, let us consider the correction to the ground state 
energy En^q. In this case the matrix elements Vkfi corre- 
spond to electron tunneling between the initially neutral 
neighboring grains. The excitation energy of this process 
equals e^^ -I- e^-g -I- E^h, where e^^ is the bare (with no 
Coulomb energy included) energy of an electron excita- 
tion in grain 2, Sk-^ is the bare energy of a hole excita- 
tion in grain 1, and E^h is the electrostatic energy the 
electron-hole excitation that should be determined from 
Eq. 1223 ■ It is equal to E^h = Ef^ + E^^ - 2-5^2 > where 
Eij can be obtained from Eq. H2.23I) putting rii = 1 and 
Uj — —1. For the periodic array of the grains under 
consideration, the energy E^h reduces to 

Eeh = 2Ec-2E^^. (2.93) 

The energy correction corresponding to such a process is 



ki ,k2 



£ki + £fc2 + Eeh ' 



(2.94) 



where the factor 2 takes into account the equivalent pro- 
cess of the electron hopping from grain 2 to grain 1 . 



Analogously, we find the correction to the energy 
En^i. The excitation energy of the process of the elec- 
tron tunneling from the grain 2 to the grain 1 equals 
£fei + £k2 + 3£;fi + E^^ - 4£;f2 = + £/c2 + ^E^h, while 
the excitation energy of the process of electron tunneling 
from the grain 1 to the grain 2 is ek^ + Efea + Efi — E22 — 
£ki + Efea- The corresponding correction to the energy 
i?jv=i reads 



AE 



N=l 



E 

ki ,^2 



\tkit 



£ki + £k2 



2E, 



£ki + £k2 



(2.95) 



One can see that corrections to the energy levels, 
Eqs. H2.94I I2.95|l . are ultraviolet divergent. However, 
their difference is finite and it gives the correction to the 
Mott gap. Am- Subtracting Eq. ifTMjl from Eq. 
replacing the summation over the states by integrals over 
a continuous variable e, such that Sk^(k2) ~^ ^1(2)1 we can 
easily caclulate the Mott gap. Am- Taking into account 
the spin degenera cy we obtain the following expression 
for the Mott gap l)Beloborodov et a^J . 12005^ 



2z 

Am = Ec gEeh In 2, 

TT 



(2.96) 



where z is the coordination number of the array of the 
grains and Eeh is the energy necessary to create an 
electron-hole excitation in the system by removing an 
electron from a given grain and putting it on a neighbor- 
ing one. In the case of a diagonal Coulomb interaction 
Eij — Ec 6ij the energy of an electron-hole excitation E^h 
is simply twice the Coulomb charging energy E^.- 

Derivation of Eq. (|2.96|l was carried out neglecting the 
fact that an extra electron added to the neutral system in 
the presence of a finite intergranular coupling can move 
diffusively over the sample. Contributions corresponding 
to such processes are suppressed by an extra small factor 
6/Ec ^ 1 and, thus, can be neglect ed. More details of 
the ca lculations can be found in Ref. <|Beloborodov et all . 
[20053). 



3. Mott gap at large tunneling conductances 

From Eq. H2.96|l one can see that the Mott gap Aj\/ is 
significantly reduced at values ~ 1, where the pertur- 
bation theory becomes inapplicable. Suppression of the 
gap Aj\/ at large tunneling conductances g can be found 
with the help of the phase action, Eqs. H2.55I - [!i.57|l . 
In fact, the renormalization of the charging energy at 
gz ^ I can be obtained with the renormalization group 
approach that we have described in subsection III. C . 2l 

Indeed, the RG Eq. (|2.68|l is written assuming that 
g is a function of the independent variable ^. We can 
invert this equation and assume that ^ is a function of g. 
The running variable ^ determines an effective charging 
energy E^-^-f as ^ = InE^-^^ and we obtain solving the 
latter equation 



Eff ~ Aexpi 



-ngz) 



(2.97) 
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where A is a constant. 

In order to determine the Mott gap Am we notice that 
at conductances ~ 1 the charging energy is of the 
order of the Mott gap Au, as one can see from Eq. H2.96|l . 
Matching Eq. with Eq. ifT^ we conclude that the 

Mott gap Am at large gz is reduced exponentially and 
can be written as 

Am ^ Ec exp {-ngz) . (2.98) 

The RG approach used to obtain Eq. (|2.98(l is applicable 
only at energy scales larger than the inverse escape rate 
from a single grains F, which means that Eq. H2.98|l is 
valid as long as the Mott gap Am is larger than T. In 
this region the system should be the insulator. In the 
opposite limit one can expect that the granular material 
becomes a metal. 

Taking Am ^ E and resolving Eq. (|2.98() with respect 
to the conductance g we obtain the critical conductance 
gc, Eq. I|2.12|l . that marks the boundary between the insu- 
lating and metallic phases at low temperature. Thus, we 
come to the conclusion that the Mott gap Am remains 
finite as long as g < g^ which assumes the activation 
behavior, Eq. H2.13|l . The Arrhenius law, Eq. H2.13|l . is 
typical for crystalline insulators. 



in both the cases, and one may expect that the meth- 
ods developed for the study of the Mott transition in the 
system of strongly interacting electrons on the lattice. 

In the case of Id and 2c? granular samples the ques- 
tion of the Mott transition at T ^ cannot be for- 
mulated since, even in the metallic phase, the conduc- 
tivity corrections are divergent at low temperatures, see 
Eqs. (|2.6() - (|2.9|l . However, the critical conductance gc 
still has a meaning of the boundary between the hard 
gap insulator behavior at g < gc and a weak insulating 
behavior g > gc that essentially coincide with insulating 
properties of the ordinary disordered metals. We note 
that the fate of the ordinary disordered metals with in- 
teraction at low temperature is not quite clear due to 
difhculties r elated to the dive rgence of the conductivity 
corrections l)Finkelsteiri Il990() . In particular, the metal 
to ins ulator transition recently observed in disordered 
films Hb rahams et al\ . l200ll: iKra vchenk o and SarachiH 
l2004j) has not been yet theoretically understood, al- 
though some pla usible scenarios have been sugges ted re- 
cently (see, e.g. ijPunnoose and Finkelsteinll2005|) '). 

It is clear that understanding the metal-insulator tran- 
sition in the granular materials is one of the most difficult 
theoretical problems and considerable efforts are neces- 
sary to solve it. 



4. Metal - Insulator transition in periodic granular arrays 

In the previous subsection we have shown that in the 
case of three dimensional arrays there is a critical conduc- 
tance gc, such that the samples with g < gc are insulators 
at T — !■ 0, while those with g > gc are metals. The analy- 
sis was performed from both metallic and insulator sides 
and both the approaches agree as concerns the value of 
the critical conductance gc, Eq. H2.12|l . Unfortunately, 
they are not applicable in the vicinity of the phase tran- 
sition and do not allow us to find the critical behavior of 
the Mott gap Am near gc- 

The order of the Mott transition in 3d periodic gran- 
ular array is also not known, since one cannot exclude 
a weakly first order phase transition. For, example the 
Mott transition in the Hubbard model in dimensionality 
d > 3 is believed to be of the first order, and since physics 
of the transition is similar, the same can be expected for 
the periodic granular array. We note however that in 
spite of the similarities between the periodic granular ar- 
ray model and the Hubbard model, there are essential 
differences between the two. Namely, in granular metals, 
even in the case of the periodic samples, the electron mo- 
tion within a grain as well as on the scales exceeding the 
intergranular distance is diffusive, while in the Hubbard 
model electrons move through a periodic lattice that al- 
lows to label their states by quasimomenta. 

The model of the granular metals has also an addi- 
tional physical parameter 6- the mean energy level spac- 
ing in a single grain, that has no analog in the Hubbard 
model. Yet, the physics of the transition looks similar 



G. Hopping conductivity in Granular Materials 

In the previous section we have shown that a strictly 
periodic granular model predicts and explains the ac- 
tivation conductivity behavior only. Experimentally 
observed temperature behavior of the conductivity of 
granular metals and arrays of quantum dots at low 
temperatures, however, is usually not of the acti- 
vation typ e but resembles the Efros-Shklovski i law, 
Eg. itm. l)Liao et all I200I iRomero and Drndil 120051: 
iTran et ad l2005HYakimov et a/.ll2003HYu et a^.ll2004^ 
As we have already discussed in the introduction, it is 
crucial for understanding the hopping transport in gran- 
ular metals to take into account an electrostatic disorder 
that seems to be present in any realistic system. This 
disorder can be viewed as a random potential Vi applied 
to each grain that lifts the Coulomb blockade on a part 
of grains in the sample such that the conductivity is me- 
diated by the electron hopping between the sites where 
the Coulomb blockade is almost removed. Below, we 
consider in more details the two essential ingredients of 
the hopping conductivity - the density of states and the 
mechanisms of the electron tunneling through the dense 
granular system in different regimes. 



1. Density of states 

The electrostatic disorder causes fluctuations in the 
electrostatic energy of granules and can thus lift the 
Coulomb blockade at some sites of the granular sam- 
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pie. This results, in its turn, in the finite density of 
states at the Fermi level and makes the variable range 
hopping the dominating mechanism of conductivity. The 
bare density of states induced by the random potential 
can be substantially suppressed due to the presence of 
the long-rang Coulomb interaction in the same way as it 
happens in semiconductors where the density of states 
(DOS ) Vn(s) is given by t he E f ros-Shklovskii expres- 
sion ijEfros and Shklovskil Il975t IShklovskii and Efrol 
Il988r 

:.,(£) ^ {k/eY \ef-\ (2.99) 

with e being the electron charge and k being the dielec- 
tric constant. This result has recently been confirmed 
bv iMu Ucr and loffc (2004); Pankov and Dobrosavliicvic 
1)20051) analytically using a locator approximation. 

As we have explained in Sec. III.BI the Coulomb inter- 
action in the granular matter may be considered writing 
classical electrostatic formulae for the electron-electron 
interaction, Eq. (|2.23l) . This approach can also be used 
in the limit of a weak coupling between the grains, g -^l. 
So, we describe the Coulomb interaction in the presence 
of the electrostatic disorder writing the following expres- 
sion for the Hamiltonian H 

H^J2 + II ^<>J ' (2-100) 

i ij 

where Vi is the random external potential, Ui is the num- 
ber of excess electrons on the grain i, and E![j is the 
Coulomb interaction between the grains i and j. Taking 
into account the asymptotic behavior of the Coulomb 
interaction matrix Ef^ ~ e^/2ry K, where rij is the dis- 
tance between the grains i and j and k is the effective 
dielectric constant of the granular sample we see that 
the classical model (|2.100|) for the granular materials is 
essentially equivalent to the one studied by Efros and 
Shklovskii. Therefore their results have to be applicable 
to the model of the granular array, Eq. p.lOOfl . as well. 

At the same time, one expects that the DOS in the 
array of the metallic granules is larger than that in a 
semiconductor since each metallic grain has a dense elec- 
tron spectrum. Indeed, one has to remember that there 
are many electron states in a grain that correspond to 
the same grain charge, while in the model of impurity 
levels in the semiconductors the charge is uniquely (up 
to the spin) identified with the electron state. 

The energy of an unoccupied state Si in the model spec- 
ified by Eq. H2.100|l is by definition the energy of an elec- 
tron placed on this state. In the granular metal, any 
state with the energy larger than Si is also available for 
the electron. Thus, in order to translate the ES result, 
Eq. H2.99|l . to the density of electron spectrum in gran- 
ular metals one has to integrate the dependence (|2.99|) 
over the energy s and multiply it by the bare DOS in a 
single grain. As a result, we obtain 

iy{e)^iyoi\e\k/eY, (2.101) 



where vq is the average DOS in a single grain (defined 
as the number of states per energy). However, the above 
expression cannot be used in the Mott argument for the 
hopping conductivity where one needs to estimate the 
distance to the first available state r within the energy 
shell £ via the relation 

r"^ f de'vg{e') ~ 1. (2.102) 

The problem with using the expression for DOS, 
Eq. H2.101|) . in Eq. (|2.102|l is that the expression (|2.101|) 
takes into account the fact that if there is a state avail- 
able for placement of an electron with a given energy 
then typically on the same grain, there will be plenty of 
other states available for electron placement. However, 
for application to the hopping conductivity we should not 
count different electron states that belong to the same 
grain since it is enough to find at least one state to in- 
sure the transport. 

Thus, when finding DOS relevant for the hopping 
transport, the lowest energy states within each grain are 
to be counted only. Then, we arrive at the conclusion 
that, although the electron DOS in granular metals is 
modified according to Eq. (|2.101|l . one has to use the ES 
expression for DOS in its form (|2.99|l even in granular 
metals in order to find the distance to the first available 
state within a given energy shell via Eq. (|2.102|l . 

Similar considerat i ons were presented 

in IZhang and Shklovskiil l)2004j) . Following this ref- 
erence we call the DOS that counts only the lowest 
excited states in each grain and that is relevant for the 
hopping conductivity "the density of ground states" . In 
order to distinguish this quantity from the conventional 
density of states from Eq. (|2.101|l we ascribe to it the 
subscript g as in Eq. (|2.99|) . 

The presented arguments demonstrate that one can 
obtain a finite density of states in a granular material 
in the insulating regime. Although this is the necessary 
ingredient for establishing the mechanism of the hopping 
conductivity, it alone is not sufficient for this type of 
transport. The problem is that en electron has to hop 
over several grains, which, at first glance, does not look 
probable in the case of closely packed granular array. For 
quite a long time, this fact did not allow to apply the 
Efros-Shklovskii theory for explanation of the behavior, 
Eq. (|l.l|l , in th e granu l ar ma terials and a rather artificial 
model of .Abeles et al ] lll975D was used. 

Only recently a resolution of this puzzle has been sug- 
gested independently in publications fBel qborodov et all . 
l2005d iFeigel'man and loselevick ■2005fK The 
authors of these works demonstrated that a 
well known phenomeno n of co -tu nnelin g (vir- 
tual hops) llAve rin a nd Nazarovl Il992l Il990t 
lAverin and Likharevl Il99lir through a grain might 
be responsible for the long range hops of the electrons. 
One has to distinguish between elastic and inelastic 
co-tunneling processes. The elastic co-tunneling is the 
dominant mechanism for the hopping conductivity at low 
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temperatures T < Tc, 



while at larger temperatures 



T > Tcross the electron transport goes via inelastic 
co-tunneling processes. The characteristic temperature 
Tcross of the crossover from the elastic to the inelastic 
tunneling is given by the following formula 



T 



(2.103) 



where c « 0.1 is the numerical coefficient. In the next 
subsections we discuss the elastic and inelastic contribu- 
tions to the variable range hopping. 



following combination of the electron and hole excitation 
energies 



(2.107) 



The presence of the delta function in Eq. H2.1()5|l reflects 
the fact that the tunneling process is elastic. 

The form of Eq. (|2.105|l for the tunneling probability 
Pel corresponds to an independent sequential tunneling 
from grain to grain. This equation means that on average 
the probability falls off exponentially with the distance s 
along the path: 



2. Hopping via elastic co-tunneling 

Considering the probability for an electron to tunnel 
from the site iq to the site in+i it is convenient to put 
the Coulomb interaction energy at these sites to zero and 
count the initial, and final, ^at+i, electron energies 
from the Fermi level. The presence of the electrostatic 
disorder on the grains is modelled by the random poten- 
tial Vi. The energy of the electron (hole) excitation on 
the site i is 



E^ 



Ef±fi,, 



(2.104) 



where = Vi + 2 El^ nj is the local potential that 
along with the bare potential Vi includes the potential in- 
duced by the neighboring grains. The probability P^i of 
a tunneling process via the elastic cotunncling can most 
easily be found for the case of the diagonal Coulomb in- 



teraction Ef, = E- 



Sij. We leave calculational details for 



the Appendix A and use here only the final result of the 
derivation. 

The probability of the elastic cotunncling through N 
grains can be written as 



Pel ~ exp(-2s/^e;) , 



(2.108) 



where the dimensionless localization length S,ei can be 
written as 



ln{ET:/cgS)' 



(2.109) 



Both the distance s and the length £^ei m Eqs. (|2.1()8I 
I2.109|l are measured in the units of the grain size a. The 
numerical constant c in Eq. (|2.109|l is equal to unity for 
the model of the diagonal Coulomb interaction, Efj — 
Ec Sij . We show in Appendix A that the inclusion of the 
off-diagonal part of the Coulomb interaction results in 
a renormalization of the constant c to a certain value 
0.5 < c < 1. 

Applying the conventional Mott-Efros-Shklovskii argu- 
ments, i.e. optimizing the full hopping probability that 
is proportional to 



exp[-(2s/^e;) - {e'/kTas)], 



(2.110) 



one obtains Eq. H2.14|l for the hopping conductivity with 
the characteristic temperature 



r 



Pel ^wgo [^] S{^N+i 



N 



(2.105) 



where the factor w — n{S,o)[l — ?^(^Af+i)] takes into ac- 
count the occupations n(^o) and n(^jv+i) of the initial 
and final states respectively (the initial state is filled and 
the final one is empty). The bar in Eq. H2.105|l denotes 



the geometrical average of the physical quantity along 
the tunneling path. For example, the average energy F 
is defined as 



InF = 



1 



iV+ 1 



N 

E 

k=Q 



InFfe, 



(2.106) 



where the summation extends over the tunneling path, 
T/c = gk^k and gk is the tunneling conductance between 
the fc-th and (fc-l- l)-st grains. [We note that the meaning 
of Ti as an escape rate from a grain does not hold in the 
insulating regime under consideration.] Yet, we use this 
notation even in the insulating sate. The energy E is 
the geometrical average \nE = X^feLi ln_Efc, of the 



Tn 



(2.111) 



where k is the effective dielectric constant of a granular 
sample, a is the average grain size, and ^e/ is given by 
Eq. (|2.109l) (when deriving Eq. (|2.111l) we considered the 
tunneling path as nearly straight). 

In the presence of a stron g electric field £ , a direct 
application of the results of ShklovskTil l)l973|) gives in 
the limit of low temperatures 



T ^ VE^ 



eS,eia 



the following expression for the nonlinear current depen- 
dence 



J ^ Jo exp 



1/2 



(2.112) 



where the characteristic electric field Eq is given by the 
expression 



£o To/eaS,ei. 



(2.113) 
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In a full a nalogy with the hopping cond uctivity in semi- 
conductors ijShklovskii and Efrosl Il988|) . also in granu- 
lar metals inelastic processes are required to allow an 
electron to tunnel to a state with a higher energy. In 
the granular metals such processes occur due to interac- 
tion with phonons as well as due to inelastic collisions 
with other electrons. It is clear that these processes 
were not considered when deriving Eq. H2.1(J5|I . There- 
fore, the results presented in this subsection are valid 
as long as the contribution of the inelastic cotunneling 
to the hopping conductivity can be neglected. This cor- 
responds to the limit of low temperatures and electric 
fields T,£ea <C Tcross- Below, we present the results for 
the opposite limit when the inelastic processes give the 
main contribution to the hopping conductivity. 



3. Hopping via inelastic co-tunneling 

The inelastic cotunneling is the process when the 
charge is transferred by different electrons on each ele- 
mentary hop. During such a process the energy trans- 
ferred to a grain by the incoming electron differs from 
the energy taken away by the outgoing one. 

As in the case of the elastic cotunneling we assume that 
the electron tunnels from the grain iq with the energy 
to the grain iN+i with the energy ^at+i. For the proba- 
bility of such a tunneling process through a chain of the 
grains we find the following exp ression [see Appendix B 
and Ref. (Bclobor odov et a/.L 12005c')] 



AttT 



E 



2N 



|r(iv + ^)| 



r(2A^) 



A 

2T 
(2.114) 

where r(a;) is the Gamma function and A = S,n+i — Co is 
the energy difference between the final and initial states. 
The appearance of the factor exp (— A/2T) is consistent 
with the detailed balance principle. Indeed, at finite tem- 
peratures an electron can tunnel in two ways: either with 
an increase or with a decrease of its energy. Accord- 
ing to the detailed balance principle the ratio of such 
probabilities must be exp(A/T) , that is indeed the case 
for the function Pin, (|2.114|l . since apart from the factor 
exp(— A/2r) the rest of the equation is even in A. 

In order to obtain the expression for the hopping con- 
ductivity one has to optimize Eq. 1)2.11411 with respect to 
the hopping distance N under the constraint 



1, 



(2.115) 



that follows from the ES expression for the density of the 
ground states (|2.99|) . Optimization of Eq. H2.114|l is a 
bit more involved procedure than the standard deriva- 
tion of Mott-Efros-Shklovskii law based on the Gibbs 
energy distribution function (see Appendix B). Never- 
theless, it leads to essentially the same result, i.e. the 
ES law, Eq. (|2.14() , where the characteristic temperature 
To (T) takes the form 



To(T) ~ e^/ak^^n{T), 



(2.116) 



with a dimensionless weakly temperature dependent lo- 
calization length ^in(T), 



ln[E^/16cTTT'^gy 



(2.117) 



where the coefficient c equals unity for a model with a 
diagonal Coulomb interaction matrix. 

Comparing the characteristic lengths ^e; (T) and 
£.in (T) for the elastic and inelastic processes, respec- 
tively, we see easily that £,ei (T) > £^in (T) at tempera- 
tures T < Tcross, where the temperature Tcross is given 
by Eq. H2.103|) . This means that the mechanism of the 
inelastic scattering is more efficient at not very low tem- 
peratures exceeding Tcross- 

As in the case of the elastic cotunneling, the long range 
electrostatic interaction results in the reduction of the 
coefficient to some 1 /4 < c < 1 . 

At zero temperature. Equation (|2.114|l for the inelastic 
probability Pin can be simplified for A < to the form 



P,n{T = 0) 



w2^^Tr |A 



2N-1 



{2N-l)l E^^ 



N+l 



(2.118) 



while for A > 0, Eq. 12.1141) gives zero since the electron 
tunneling process with an increase of the electron energy 
is prohibited at zero temperature. 

In order to obtain the hopping conductivity in the 
regime of a stro ng electric field £ and low temperature, 
we use, following IShklovskiill|l97.'^ . the condition (|TTTH)) 
that defines the distance to the first available electron site 
within the energy shell A. We use also the equation 



e£r 



(2.119) 



that relates A to the electric field £, and the fact that 
the distance between the initial and final tunneling sites 
is r ^ Na. Equations (|2.115|) and (|2.119|) allow us to de- 
termine the quantities A ^ £e^ jk and N ^ ^ elk£o? . 
Inserting their values into Eq. H2.118() we come to the 
following expression for the current 



J ~ Jo exp 



-(.£^l£fl^ 



(2.120) 



where the characteristic electric field £q is a weak func- 
tion of the applied field £ 



fo(f)^^ ln'[^VeWg]. 



(2.121) 



The results obtained in the insulating regime show that 
the hopping conductivity is the main mechanism of trans- 
port in the low temperature regime. The temperature 
dependence of the conductivity is determined by the 
Efros-Shklovskii law, Eq. H2.14|l . and is similar to the one 
in amorphous semiconductors. The exponential depen- 
dence of the current on a strong electric field is given by 
Eqs. (|2.112() and H2.120() . Again, this is the same depen- 
dence as in the amorphous semiconductors. Of course, 
the characteristic temperature Tg and electric field £o 
are model dependent. 
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The theory for the variable range hopping (VRH) via 
virtual electron tunneling through many grains is based 
on the assumption that the hopping length r* exceeds the 
size of a single grain a. This hopping length r* should 
be found minimizing the exponent in Eq. (|2.110|l . From 
this condition we see that the characteristic length r* (T) 
equals 

.*(T) = /^, (2.122) 

where the dimensionless localization length ^ is given by 
either Eq. H2.109|l or Eq. H2.117|) . and that it decreases 
with increasing the temperature. At temperatures T ~ 
e^^/2aK the hopping length becomes of the order of the 
grain size, r*{T) ^ a. At such distances the VRH pic- 
ture does not work anymore because the electrons hope 
between the neighboring grains only. This means that 
at comparatively high temperatures T <^ T <^ Am one 
should use the Arrhenius law, Eq. H2.13|l instead of the 
Efros-Shklovskii law, Eq. H2.14II . Such a crossover has 
been observed experimentally ( Tra n et al\ . [2005) . 

Thus, we have obtained a remarkable result: Even 
though the array of the grains may look very regular, 
i.e. have a periodic structure and equal size and shape 
of the grains, it is almost inevitable that at low tempera- 
tures the main mechanism of the conduction is the vari- 
able range electron hopping and the conductivity is de- 
termined by Eq. H2.14|l as if the system were amorphous. 
This concerns also other formulae like Eqs. ()2. 112(1 and 

mm . 

III. ARRAYS COMPOSED OF SUPERCONDUCTING 
GRAINS 

In this Section we consider properties of granular ma- 
terials consisting of superconducting particles. As in the 
Section ^ devoted to study of normal metals we start 
our presentation in the next subsection with a qualita- 
tive discussion of the most interesting effects that occur 
in these systems. A quantitative analysis of these effects 
will be given in the subsequent subsections. 

A. General properties of granular superconductors 

1. Single grain 

A granular superconductor, similarly to a granular 
metal, can be viewed as an array of superconducting 
granules coupled via electron tunneling. Superconduct- 
ing properties of an array in many ways are determined 
by the properties of granules which it is formed of. 
For this reason we begin our discussion with reviewing 
properties of superconductivity in a single small isolated 
grain. 

Th is question wa s first addressed by Ander- 
son ijAndersorJ . Il959(l who realized that the s-wave 



superconductivity is almost insensitive to disorder of 
general kind as long as it does not break the time reversal 
invariance. This includes, in particular, the diffusive 
scattering by the grain boundaries that is in many ways 
similar to the scattering by potential impurities in the 
bulk. As long as the time reversal invariance is not 
broken, one can choose the basis of exact eigenfunctions 
and then apply the Bardeen-Cooper-Schrieffer (BCS) 
theory in the usual way. The critical temperature can 
still be expressed via the BCS effective coupling constant 
and the density of states in the vicinity of the chemical 
potential. The BCS coupling constant does not change 
considerably when changing the basis and the average 
density of states is not sen sitive to disorder. Thus, one 
concludes ()Andersorlll959fl that the critical temperature 
of the single grain has to be close to the bulk value. 

Anderson has also pointed out that these arguments 
hold even for small grains provided the average distance 
between the energy levels 5, Eq. (|1.3() is s till smaller 
than the superconducting gap A, in the bulk l|Andersonl 
I1959D . As soon as the energy level spacing 5 reaches the 
superconducting gap, the conventional BCS theory is no 
longer applicable. 

This fact can easily be understood from the BCS equa- 
tion for the order parameter that can be written as 

A^^oTVy ^ 2, (3.1) 

where Aq is the standard dimensionless constant de- 
scribing the attraction between the electrons, lo = 
2ttT {n + 1/2) are Matsubara frequencies and 5 is the av- 
erage mean energy level spacing in a grain determined by 
Eq. 1(1. 3|l . The discrete variable ^fc equals 

Cfc = <^k - Ai, (3-2) 

where are energies of the electron states in the grains 
and /i is the chemical potential. 

If the grain size is very large, such that the mean level 
spacing 5 is much smaller than the superconducting crit- 
ical temperature TcOi the sum over the states in Eq. 1(3. 1() 
can be replaced by the integral 




and we come to the conventional BCS equation. 

At the same time, it is clear from Eq. (13.1(1 that a fi- 
nite difference between the energy levels ^j, plays a simi- 
lar role as a finite temperature, namely, it reduces the 
superconducting gap |A|. As soon as the mean level 
spacing 5 exceeds Tco (or A (T = 0) for the bulk su- 
perconductors), the superconducting gap |A| vanishes. 
To be more precise, it remains finite due to fluctua- 
tions but its value is sufficiently reduced and in this 
regime the grain cannot be considered as a "true" super- 
conductor. Nevertheless, the superconducting pairing is 



29 



still present manifest i ng:, for example, in the parity ef- 
fect, llvon Delft et all Il996t iMatveev and Larkhl Il997t 
ISmith anc^^nbegaokaii Il99(j) . 

The Anderson theory agrees well with experiments. 
In particular, electron spectrum of ult rasmall su- 
percon d ucting Al grains was st u died b y iBlack et al. 



1 19961) : iDavidovic and Tinkhainl l)l999|) : iRalph et all Il98(]|) leading to the result at T = 



1 19951) where a well defined superconducting gap was 
shown to survive down to grain sizes lOnm, while 
smaller grains < 5nm did not show any superconduct- 
ing features in their electron spectrum. 

It is worth noting that, while in the diffusive/chaotic 
grains the density of states in the vicinity of the chemical 
potential indeed well agrees with the corresponding bulk 
values, this is not necessarily the case for the ballistic 
grains possessing certain geometric symmetries. For ex- 
ample, in a spherical grain one naturally expects that the 
energy levels are strongly degenerate. This may change 
the density of states on the Fermi level significantly re- 
sulting in a considerable change of the transition tem- 
perature. The possibility to utilize the enhanced density 
of states that may appear from the orbital degeneracy 
in order to increase the single grain critical temperature 
was rece ntly suggested in Refs. l|Ovchinnikov and Kresinl 
1200533). 

Although all these effects are very interesting, in this 
review we follow our line of studying models with pa- 
rameters that can easily be achieved in existing granular 
materials. So, we assume everywhere that granules are 
sufficiently large to maintain the local superconductiv- 
ity within each grain, i.e. that the Anderson criterion is 
satisfied 



5< A. 



(3.4) 



We assume also that the grains are diffusive or do not 
have an ideal shape like spheres or cubes, which excludes 
the possibility of an accidental degeneracy of the electron 
spectrum. 

As a final remark in our discussion of superconductiv- 
ity in a single grain, we would like to point out that the 
thermodynamics of a single grain is not affected by the 
presence of the Coulomb interaction as long as the latter 
can be expressed in a standard way through the opera- 
tor of the total electron number. The reason is that the 
number of particles in an isolated grain is conserved and 
thus, the Coulomb term in the Hamiltonian commutes 
with the rest of the Hamiltonian and does not influence 
the thermodynamic quantities. 

At the same time, phase fluctuations of the order pa- 
rameter A (r) 



A(t) = |Ao|exp [2iip{T) 



(3.5) 



are very sensitive to the presence of the charging term, 
Eq. (1223 ■ These fluctuations can be conveniently char- 
acterized by the correlation function Hg (t) 



where the averaging should be performed over the states 
of the Hamiltonian including the Coulomb interaction, 
Eq. 1223. 

Actually, the average over the states of the Hamilto- 
nian in Eq. (|3.6|l reduces to an average with the ac tion Sr. , 
Eg. I|2.56|) . and can rather easily be calculated ijEfetovl 



(3.7) 



H, (t) = exp i-AErT) 



where Ec — En is the Coulomb energy of a single grain 
and a factor 4 appears because of the doubled Cooper 
pair charge. Equation (|3.7I) is written for the imaginary 
time T. Changing to the real time, r — > it, one can 
conclude that the correlation function H^ (t) oscillates 
in time. Although the exponential form of the function 
Hs (r) does not lead to changes in the thermodynamics, 
it is crucial for the problem of the coherence in an array 
of many grains. 



2. Macroscopic superconductivity 

Having listed the superconducting properties of a sin- 
gle grain we turn now to discussing the superconductivity 
in the whole sample. In order to simplify the presenta- 
tion we first restrict ourselves to the zero temperature 
case. 

It is clear that an array with sufficiently strongly cou- 
pled grains should be able to maintain the supercon- 
ducting coherence in a whole sample because the cou- 
pling reduces the phase fiuctuations. On the contrary, 
in the opposite limit of a weak coupling, one expects 
that the strong Coulomb interaction should lead to the 
Coulomb blockade of the Cooper pairs in analogy with 
the Coulomb blockade of electrons in granular metals in 
the low coupling regime. 

In o rder t o quantify this intuitive statement 

Abeles jAbelesl [ 1 97 7') . following the earlier idea of An- 
derson IjAndersonl . Il964) . suggested to compare the en- 
ergy of the Josephson coupling of the neighboring grains 
with the Coulomb energy. Indeed, the Josephson cou- 
pling tends to lock the phases of the neighboring grains 
and to delocalizc the Cooper pairs, while the Coulomb 
interaction tends to localize the Cooper pairs and thus 
enhance the quantum phase fluctuations. 

Comparing the Josephson energy J, 



J = 7r.gA/2, 



(3.8) 



H, (r) = (exp (2* (r) - ^ (0)])) , 



(3.6) 



with the Coulomb energy one would come then to the 
conclusion that samples with g > ~ Ec/A. should 
be superconductors, while those with g < gs - insu- 
lators. This result was later derived within a model 
for a superconducting granular array including both 
Coulomb and Josephso n interactions in a number of the- 
oretical works jEfetovl lTgSO: McLean and Stephenl . ll979t 
ISimaned Il979(l using different methods. We will review 
the main results of these studies in the Sec. IIII.CI 
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FIG. 12 Resistance as a function of tempera t ure of one of the 
granular samples of IShapira and Deutsche'il il983tl . In spite 
of the fact that the sample exhibits strong insulating behavior 
in the intermediate temperature range it finally turns to the 
superconducting state. 



The Anderson- Abeles criterion predicts that the su- 
perconductivity may survive even in such samples that 
would be insulators if they were made from normal grains 
under equivalent conditions. Indeed, for samples with a 
small ratio E^j ^ <C 1, there is a parametrically large 
interval of conductances E^l ^ <C 5 <C 1 where super- 
conductivity should exist in spite of the fact that the 
corresponding normal array would be an insulator. 

Qualitativel y, this prediction was s upported by exper- 
iments ( Shapira and Deutscheil ll98S) where some su- 
perconducting samples were found to show clear insu- 
lating behavior above the critical temperature. (We 
note, however, that the ratio Ed ^ in the samples used 
in this experiment was certainly large, this fact, as we 
will show below, may be explained via renormaliza- 
tion of the Coulomb energy). The corresponding de- 
pe ndence of the resistance on temperature according 
to lShapira and Deutsche J l)l983|) is shown in Fig.lT^ 

However, the Anderson-Abeles criterion was found 
later to be in conflict with man y other experiments. 
After the pioneering works (,Jaeger et alV [1989,. ,1986; 
lOrr et al 1 I198,'t[ 119861) it became clear that the bound- 
ary between the superconducting and insulating states at 
T = is determined by the normal state conductance of 
the film rather then by the ratio of the Josephson and 
Coulomb energies. 

This fact is illustrated by the experimental resistiv- 
ity temperature dependencies of samples with differ- 
ent thicknesses shown in Fig. Ga and Pb ultrathin 
films. One can see that there exists a critical value 
of the normal state resistance Rq close to the value of 
quantum resistance /i/(2e)^ « 6.4k such that sam- 
ples with resistance R < Rq eventually become su- 
perconducting under temperature decrease, while those 
with R > Rq show tendency to the insulating behav- 
ior. Such a behavior turned out to be generic and was 
observed in a number of granu lar (as well a s on homo- 
geneously disordered, (Frvdm an et all l2002t iLiu et gj] . 
Il993t iMarkovic et all Il999f) l films made from different 
materials. 




FIG. 13 Temperature dependence of resistivit y of granular 
(a) G a and (b) Pb ultra thin films according to I Jaeger et al\ 
ijlflSST) . Different curves correspond to the films with different 
normal state resistances that is tuned by the sample thickness. 
Low temperature state is controlled by the film normal state 
resistance. 



In this respect, the superconductor-insulator tran- 
sition in thin granular films resembles the analogous 
trans ition in a resist i vely sh u nted Josephson junc- 
tion l|Bulgadaevi Il984l: ISchmidl . Il98.'^ where the state 
of the system is known to be controlled by the value 
of the shunt resistance only, no matter what the ratio 
of t he Coulomb and Josephson energies is (for review 
see ISchon and ZaikinI §990)). The resistive shunt in 
this ap proach is included via the C aldeira-Leggett ap- 
proach l)Caldeira and Leggetti . ll98ll [198381 that allows to 
take into account the dissipative processes on the quan- 
tum mechanical level. Similar approach based on the 
inclusion of the phenomenological dissipative term that 
models the resistive coupling of the grains was applied 
to the granular system in order to resolve the disagree- 
ment with experiment llChakravartv et a^J . 1198 6: FisheJ 
119861: ISimanek and BrownL Il986(l . ' It was found that in 
the limit Ec/J > 1, in agreement with experiment, the 
boundary between the superconducting and insulating 
states is controlled by the resistance of the shunts only 
rather than by the ratio of Ec / J- 

However, the presence of the phenomenological resis- 
tive coupling between the grains is rather difficult to jus- 
tify, especially in the limit of low temperatures where 
all quasiparticles are frozen out. Th i s pro blem was fi- 
nally resolved by IChakravartv et~ai ] II198/D on the ba- 
sis of the model that took into account the direct elec- 
tron tunneling processes between the grains in addition 
to the Josephson couplings. The phase diagram obtained 
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FIG. 14 Illustration of the renormalization of the Coulomb 
energy of a normal grain placed in a contact with a bulk 
superconductor: Due to the virtual electron tunnelling pro- 
cesses the Coulomb energy of a normal grain Ec is renormal- 
ized down to the values Ec — A/2(7. The conductance of the 
contact is supposed to be large g ^ 1 and the initial Coulomb 
energy is assumed to be larger than the superconducting gap 
A. 



in this approach turned out to be similar to the one ob- 
tained within the " dissipative" model. At the same time, 
the considered model did not rely on any phenomeno- 
logical assumptions and the virtual intergrain tunneling 
processes were well described in terms of the original tun- 
neling Hamiltonian. 

Acc ording to the approach of IChakravartv et al\ 
the originally strong Coulomb interaction is re- 
duced by the electron tunneling to other grains. This 
renormalization of the Coulomb energy can be well un- 
derstood on example of a simplified model of a normal 
grain placed in a contact with the bulk superconductor 
shown in Fig. ^| It is clear that the Coulomb charging 
energy is reduced due to the possibility for an electron 
to be present virtually in the superconductor. Actually, 
this renormalization is nothing but the screening of the 
Coulomb interaction by free charges, which is usual for 
metals. 

The value of the screened charging energy in the case 
of a strong coupling between the grains, g 3> 1, and 
strong Coulomb interaction, Ec ^ A turns out to be 
~ A/g; and what is remarkable, is independent on the 
original energy Ec- Similar result holds for the su- 
perconducting grain as well. We see that in the limit 
g ^ I the Josephson energy J ~ gA is always larger 
than the effective Coulomb energy A/g meaning that 
a sample with g ^ 1 should always b e a supercon- 
ducto r at sufficiently low temperatu re l)Eckern et all 
Il984t iLarkin and Ovchinnikovl I1983D . Another, im- 
portant conclusion is that the superconductor-insulator 
transition occurs at g ~ 1, which is in agreement 
wi th the experiments. Th e phase diagram obtained 
bv IChakravartv et all l)l987l) is presented in Fig. El 

The fact that an array consisting of strongly coupled 
superconducting grains (5 3> 1) has to be in a super- 
conducting state at sufficiently low temperatures can be 
understood from a somewhat more general consideration. 
In the strong coupling regime, the granular system resem- 
bles the homogeneously disordered system. According to 
Anderson (Anderson, 1959), within the accuracy of the 
mean field BCS theory, the thermodynamics of such a 
system is essentially the same as that of a pure bulk elec- 




FIG. 15 Zero temperature phase diagram of the granular su - 
perconducting array according to IChakravartv et al\ |l983). 
Dimensionless parameter a is related to the conductance g as 
a = Tvg. 



tron system. This means that the granular sample turns 
to the superconducting state immediately after the su- 
perconducting gap appears locally in each grain. 

Thus, the insulating state where the superconductiv- 
ity exists locally, while being absent in the whole sam- 
ple, may be obtained only beyond the BCS approxima- 
tion. Corrections to the BCS theory may come from the 
Coulomb interaction and superconducting fluctuations 
but they are expected to be small in the limit g :s> I. 
The fact that the corrections due to the Coulomb inter- 
action are small in normal metals in the limit g 3> 1 has 
been demonstrated in Sec.^ This result allows to come 
to the conclusion that at g 3> 1 the BCS theory is a good 
starting point for description of the granular system. 

This in particular means that the critical temperature 
of a good conducting granular sample (5 ^ 1) is approx- 
imately given by the single grain BCS critical temper- 
ature and that the corrections to Tc can be studied by 
means of the perturbation theory in 1/g. In the case of 
homogeneously disordered films such a perturbative ap- 
proach to calculation of Tc is well kno\yn jFi nkelstcinL 
ll987HMaekawa et a;J.[T98alOvchinnikovlll973lj) . Analo- 
gously, one can consider the corrections to the critical 
temperature in granular systems using a perturbation 
theory in the inverse conductance. We discuss this ap- 
proach and results obtained in Sec. IIII.DI 

In the opposite limit, 5 <C 1, the critical temperature 
can be found from an effective model that includes the 
Josephson and Coulomb interaction only. As the prob- 
ability of the tunneling from grain to grain is small, the 
effect of the screening due to the electron tunneling can 
be neglected. In particular, in the simplest case of a 
very small charging energy, one can estimate the criti- 
cal temperature Tc by comparing the Josephson energy 
J ~ Ag with the temperature, which leads to the es- 
timate Tc ~ Ag. We will consider this question within 
a mean field approximation in more detail in the next 
subsection. 
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3. Granular superconductor in a magnetic field 

Now we present new features in the behavior of a gran- 
ular array that manifest in an applied magnetic field. It 
turns out that magnetic properties of granular metals are 
very different from those of the corresponding bulk sys- 
tems even in the limit of large tunneling conductances, 
g ^ 1. In particular, there appears a new characteristic 
conductance value 
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0.16(a(5)"VAo£', 



(3.9) 



where Aq is the superconducting order parameter for a 
single grain at zero magnetic field and zero temperature 
and D is the diffusion coefficient in the grain. 

The critical magnetic field He of samples with g < g* 
is close to the critical field -fff of a single grain, while 
for samples with g > g*, it is close to the one in a ho- 
mogeneously disordered metals with the diffusion coef- 
ficient Deff = gSa^ (c.f. Eqs. (I221E3- The latter is 
the effective diffusion coefficient of the granular array. 
We note that the conductance g* corresponds to the well 
conducting samples. Indeed, expressing g* in terms of 
the Thouless energy Eth ~ D/a^ we obtain 
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1/2 



> 1. 



(3.10) 



The dependence of the critical magnetic field on the 
sample conductance and temperature will be considered 
in detail in Sec. IIII.(]I Here we will only show how the 
new conductance scale H3.9|l appears: The critical mag- 
netic field H^^ destroying the superconductivity in a sin- 
gle grain of a spherical for m due to the orbital mechanism 
is given by ijLarkinL Il965^ (the derivation is presented in 

sec nrmi 



5Ao 



(3.11) 



where the grain radius R is related to to the period of a 
cubic granular array as R = a/2. At sufficiently strong 
coupling the critical field of the sample, on the contrary, 
has to be given by the bulk value expressed via the effec- 
tive diffusion coefficient: 



^bulk " ^' 



$0 Ao 



2TT D,ff ' 



(3.12) 



where $o — hc/2e is the fiux quantum. Comparing 
Eqs. (|3.11l) and (|3.12|) one can see that the two criti- 
cal magnetic fields iJf and i/^""^ become of the same 
order at tunneling conductances g g* . 

The orbital mechanism of the destruction of the su- 
perconductivity leading to Eqs. (|3.11|) and H3.12|l is dom- 
inant provided the grains are not very small. However, 
in the limit of very small sizes of the grain, the Zee- 
man mecha nism of the destruction c an become more 
important l)Beloborodov et all l2000t iKee et all Il998t 



iTinkhaiiil . Il996fl (everywhere in this review we consider 
the s-wave singlet pairing). 

If the Zeeman mechanism dominates the orbital one 
the critical magnetic field is given by the Clogston 
value 



HI = Ao/\/2/iB 



(3.13) 



wher e ^lb = e'h/2mc is the Bohr magneton, ijCloestonL 
Il96nh . We note that the phase transition between the 
superconducting and normal states is in this case of the 
first order. 

In the opposite limit, when the superconductivity is 
destroyed mainly by the orbital mechanism, the Zee- 
man splitting has no dramatic effect on the phase dia- 
gram. Comparing the orbital and Zeeman critical mag- 
netic fields, Eqs. H3.11I IXTS|l one can see that the latter 
can be neglected as long as the grain is not too small 
a > flc with the critical grain size Oc given by 
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(3.14) 



Throughout our review we assume that the grains are 
not too small such that the phase transition within a 
singe grain remains of the second order. Of course, we 
speak about the phase transition within the mean field 
approximation. Fluctuations of the order parameter A 
smear the transition in the isolated grain. 



4. Transport properties of granular superconductors 

Experimentally observed dependencies of the resistiv- 
ity of the granular superconductors on temperature and 
applied magnetic field are very reach and all details are 
currently far from being well theoretically understood. 
In our review we discuss only briefiy the general trans- 
port properties touching in detail only the latest develop- 
ments related to the low temperature magnetoresistance 
in granular superconductors. A broader discussion of the- 
oretical advances in understanding the transport proper- 
ties of granular superconductor s (networks of Josephson 
junct ions) can be found in Ref. l)Fazio and van der Zanll 
I2001D . 

In the absence of the magnetic field, the dependence of 
the resistivity of a granular superconductor on tempera- 
ture is typically non-monotonic. One of the examples is 
shown in Fig. E| Decreasing the temperature first leads 
to a growth of the resistivity, which is due to the enhance- 
ment of the Coulomb correlations that suppress the cur- 
rent. However, the resistivity starts decreasing closer to 
the superconducting transition temperature. In this re- 
gion the conductivity is enhance d due to the fluctuation 
contr ibution of the Cooper pairs ijAslamazov and LarkinL 
I1968D . 

Granular films with the normal state resistance close 
to the critical value Rq ~ 6.4fcil exhibit often a more 
complicated reentrant behavior as it is seen in Fig. 1131 
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FIG. 16 Low temperature resistance of 3d Al granular sam- 
ples as a function of magneti c field at low temperatures, 
T < Tc, ijCerber et ali Il997ft . The two curves are shown 
for two samples with different normal temperature resistance. 
The insert shows much less pronounce peak in the resistance 
at higher temperatures, T — Tc -C Tc. 



The resistance of such samples increases again after a 
drop associated with the fluctuating conductivity. This 
eff'ect can be a consequence of the competition between 
the conductance increase due to the fluctuating Cooper 
pairs and the freezing of the excitations due to the open- 
ing of the superconducting gap in the density of states. 

The systematic theoretical description of this phe- 
nomenon is not available yet. We can only note that 
the reentrant behavior of the conductivity seems to be 
a consequence of the granularity. Homogeneous films 
made from the same material having no (or less pro- 
nounced) granular structure do not show such a behav- 
ior (Liu et aL, 1993). 

Magnetic field applied to a granular system brings an 
extra control parameter allowing to change the resistivity 
and to tune the state of the granular system. In particu- 
lar, sufficiently strong magnetic field can always destroy 
the superconductivity even at low temperatures. 

What is interesting, the resistance of a granular su- 
perconductor has a non-monotonic dependence on the 
applied magnetic field at low temperatures. As an 
example we p resent results of a measurement on Al 
granular from Gerber et al\ (|l993). The samples used 
by ICcrbcr et al. (199jJ were three dimensional and 
monodisperse with a typical diameter of the grains « 
120A. A surprising feature of the experimental curves is 
that there is a region with a negative magnetoresistance 
and a pronounced peak in the resistivity at the magnetic 
field of several Tesla is seen. Only at extremely strong 
fields, H > 6T, the resistivity is almost independent of 
the field. 

The non-monotonic behavior of resistivity 
can be understood in terms of the corrections 
to the conductivity due to the supcrconduct- 
ing fluctuations LBc loborodov and Efetov^ 1999; 
iBeloborodov et a/1 l200d) . In the vicinity of the phase 
transition the resistivity decreases due to the opening of 
the additional transport channel via fluctuating Cooper 
pairs. This is what is called the Aslmazov-Lark in 
conductivity correction ijAslamazov and LarkinL [1968(1 . 



We note that this correction gives always a positive 
contribution to the conductivity and thus cannot explain 
the negative magnetoresistance at higher fields. How- 
ever, at stronger fields another conductivity correction 
that appears as a result of the reduction of the density 
of states due to fluctuati ng Cooper pairs begins to domi - 
nate electron transport ijBeloborodov and Efetovlll999fl . 
As a result, in this regime the conductivity becomes lower 
than that in the normal metal approaching the latter only 
in the limit H ^ He, where all the superconducting fluc- 
tuations are completely suppressed by the magnetic field. 
We discuss the magnetoresistance of granular supercon- 
ductors in more detail in Sec. IIII.EI 

Changing the magnetic field in the limit T ^ one 
can cross the boundary between the superconducting and 
metallic/insulating states and investigate the transport 
around the zero temperature (quantum) phase transition. 
The corresponding transition point at T = is referred 
as qu antum critical point ( Sachdev. ..200 k ISondhi et all 

Quantum phase transitions are analogous to the clas- 
sical phase transitions with the differences that at T — > 
the system has been described not only via the char- 
acteristic spatial (coherence length) ^ but also via the 
characteristic time scale r. Both r and ^, as in the case 
of classical transitions, are assumed to scale as powers of 
the parameter 5 that controls the closeness of the system 
to the quantum critical point. 

Scaling ideas were applied to the problem of the su- 
perconductor to insulato r transition in Refs. l|Fisheil 
1199(1 iFisher et all Il99(l that gave a substantial im- 
pact on the field. Scaling expressions for the resis- 
tivity seem to work well for many samples. How- 
ever, one of the central prediction of Refs. l|Fisheii 
1199(1 IFisher et all Il99(tl - the universal conductance at 
the superconductor to ins ulator transition is not sup- 
ported by all experiments (|C hcrvcnak_and^J^.^L_Valle£ 
200C : [Goldman and Markovid 119981: Ivan der Zant et al. 



19921) . 

Of course, the scaling approach to the resistiv- 
ity behavior is not restricted to the granular sys- 
tems and, on the contrary, homogeneously disordered 
films are much more suitable candidates for exper- 
imental study of such phenomena. For this rea- 
son we will not furthe r touch this question refer- 
ring the reader to Refs. ("Fazio and van der Zant', l200lt 
■Goldman and Markovic. .1998; .Sondhi et aL . .1997,1 . 



B. Phase diagram of granular superconductors 

In this section we discuss the phase diagram of gran- 
ular superconductors in the absence of the magnetic 
field. We will mainly concentrate on the mean field ap- 
proach and will not touch the critical behavior in the 
very vicinity of the phase transition and we will not dis- 
cuss in detail the superconductor- insulator transition. 
The latter phenomenon relates to a wide class of dis- 
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ordered superconductors and represents the field of re- 
search on its own. For reviews on this subject we re- 
fer the reader to Refs. llFazio and van der Zanti l200lt 
IGoldman and Markoviclfl99aHLarkinL ll999'). 

The arrays of the superconducting grains can be con- 
veniently described using a phase functional analogous 
to the functional S [ip] derived in Chapter ITU for the nor- 
mal granular metals. In the next subsection we introduce 
such a functional. 



1. Phase functional for the granular superconductors 

Description of a granular superconductor in terms of 
the phase action is very similar to that used for the nor- 
mal metals and the derivation can rather simply be ex- 
tended to include the superconducting order parameter. 
In order to simplify the discussion we consider the limit 
of low temperatures T <^ T^q, where T^o is the critical 
temperature in the bulk in the BCS approximation. We 
assume that the superconducting gap is still smaller than 
the Thouless energy, Eq. (|1.5|l . 

A^Eth, (3.15) 

which allows us to consider a single grain as zero- 
dimensional. Another way to write this inequality is 

Ce//»a, (3.16) 

where ^e// — -y/^o^, Co = vp/Tco, and I is the mean free 
path. 

The derivation of the action S [tp] can be carried out in 
the same way as in Chapter ^ decoupling the Coulomb 
interaction by an auxiliary field Vi and the making the 
gauge transformation. The phase ip is determined by 
Eqs. H2.44l - ITT7|l . Expanding the action in the tunneling 
amplitude Uj up to the second order we obtain the action 
S [ip] in the form of a sum of three terms 

S[ip]^S,[ip] + Sts[^]+Sj[ip] (3.17) 

The first term Sc [^p] stands for the charging energy and is 
given by Eq. H2.56(l . This means that the charging energy 
of the superconducting grain is the same as that of the 
normal one provided the same charge is placed onto the 
grains. 

The second term Sts M can be written in the form 
of Eq. H2.57|l but with another function as (r) instead 
of a (r) , Eq. (|2.58|l . In the low temperatu re limit this 
function takes the form l)Eckern et alL\l984\ 

aAr) = ^A'KUA\r\), (3.18) 

where Ki (x) is the modified Bessel function. This differ- 
ence between the functions a (r) , Eq. (|2.58l) , and as (r) 
is very important. The function as (r) (|3.18() decays ex- 
ponentially when increasing r, while the functions ^(t) 
decays only algebraically at large r. The good conver- 
gence of the kernel (|3.18|l allows us to expand the term 



Sst, Eqs. H2.57I I3.18() in powers of t — r' and to write 
the term Sts in a simplified form 

= ^ E (3-19) 

which should be correct for r ^ A~^. Only the third 
term Sj in Eq. H3.17|l is new as compared with normal 
metals; it describes the Josephson coupling between the 
grains and in the limit r ^ A^^ can be written as 




The sum in Eq. (|3.20|l should be performed over the near- 
est neighbors i , j and J — Jij for these values of i,j is 
determined by Eq. (|3.8|) . 

Equations (jXTTI OBI ITFtI IXTHl IX^ EHJ describe 
completely the granular superconductors in the low tem- 
perature limit. They can easily be extended to higher 
temperatures by writing proper expressions for the func- 
tions Us (t) and Jy . Trying to present the physical pic- 
ture in the simplest way we do not consider here these 
complications. 

The action S [tp] contains only the phases tp, whereas 
fluctuations of the modulus of the order parameter are 
neglected. This is justified at low temperatures provided 
the inequality (|3.4(l is fulfilled. A microscopic deriva- 
tion of the action Sc+^j was or iginally done using an- 
other approach in lEfetovl lll98Cil. I t was also written 
phenomenologically in ISimanekl (|l 979l) assuming however 
that the variable (p varied from — oo to oo. The exis- 
tence of the term Sts was realized later l)Eckern et al\ . 
I1984D and used for a system of Josephson junctions 
bvlChakravartv et nl\ ^987). 

Assuming that the most important are large times r 3> 
A^^ we use in the subsequent calculations Eq. H3.19|l for 
the term Sts- Then, denoting 

^0 M = Sc [ip] + Sts [^] , (3.21) 

we write the quadratic form 5*0 [ip] as 

^0 M - J E / i^c') ^ . ^^ (r) ip, (r) dT, (3.22) 
where 

Ec (q) = ^ [C (q) /e^ + (37r.g/16A)Aq] (3.23) 

and Aq is to be taken from Eq. I|2.3Q() . The function 
E (q) is the Fourier transform of the function Eij and q 
is quasimomentum. When integrating over the phases if, 
one should remember about the winding numbers. Look- 
ing at Eq. (|3.23|) one can understand that the term Sts, 
Eqs. (|2.57l I3.18|) . leads to a screening of the Coulomb 
interactions. 
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2. Mean field approximation 

In spite of many simplifications, the action 

S [^] = So [v] + Sj [^] (3.24) 

is still complicated and it is not possible to describe the 
system without making approximations. A common way 
to understand, at least qualitatively, the properties of 
a model is to develop a mean field theory. Following 
this approach one should replace the initial Hamiltonian 
(action) by a simplified one and determine parameters of 
the effective model self-consistently. 

Now we will consider the action S [ip] , Eg. I|3.24ll in 
the mean field approximation. Following lEfetovl Ill98(il) 
we make in the action Sj, Eq. (|3.2Uf) . the following re- 
placement 

Jij cos[2 {(fii — ifj)] — * Jz(cos 2(p) MF cos 2</?j. 

(3.25) 

The average {cos ip)MF in Eq. (|3.25() should be calcu- 
lated with the action Sef / [^] that can be obtained from 
S [ip], Eq. H3.24|l by making the replacement, Eq. H3.25|l . 
in Sj [ip], Eq. H3.20|l . Actually, this average is the order 
parameter for the macroscopic superconductivity. 

Assuming that the transition between the supercon- 
ductor and insulator is of the second order we can find 
the critical point from the condition that at this point 
the mean field {cos2ip) mf vanishes. The corresponding 
equation determining the boundary of the superconduct- 
ing state takes the form 

/3 

l = ^lu,{r)dT, (3.26) 



where z is the coordination number of the lattice and 
Us (r) is determined by Eq. (|3.6|l assuming however that 
the average (. . .) in that equation should be calculated 
with the action So [^p], Eqs. and ifT^ . 

In order to simplify calculations further we will neglect 
the off-diagonal terms in the matrix Eij , Eq. (|3.22(l , keep- 
ing the value Eo {g) = only. This value can be written 
as 

Eo{9)-a^ j EM)^,- (3.27) 

where Ec (q) is determined by Eq. (|3.23l) . Then for lis (t) 
instead of Eq. (|3.7|) we obtain the following expression 

n, (r) = exp {-AEot) . (3.28) 

Solving Eq. H3.26|l at T — > one obtains the relation be- 
tween the Josephson and Coulomb energies at the phase 
transition 

zJ = 8Eo (g) . (3.29) 



Equation (|3.29|) agrees at small 5 <C 1 (u p to the coefH - 
cient) with the Anderson- Abeles criterion llAbelesl . ll977i) . 
In this limit we obtain using Eq. (j3.8(l the critical value 
for the conductance in the form 

g* = 16£:o (0) /ttzA. (3.30) 

Samples with the tunneling conductances g < g* are in- 
sulators, while those with g > g* are superconductors. 
We note that this result suggests that even samples with 
weakly coupled granules g* < g 1 are still supercon- 
ductors in spite of the fact that the corresponding array of 
normal grains with the same Coulomb interaction would 
be an insulator in this regime. 

It is clear from Eq. H3.26|l that the zero temperature 
result (|3.30l) holds also at finite temperatures so long as 
T <^ Ec- In the opposite limit T :s> E^ (but stiU T < 
Tco), one obtains lis (t) = 1, which leads to the critical 
temperature 

Tc = ^zJ = ^zngA. (3.31) 

The results (I3.30|) and (|3.31() correspond to the low cou- 
pling part (small g) of the phase diagram shown in 
Fig. fTTji . 

As we have mentioned previously, the Anderson- Abeles 
criterion fails to produce the correct phase boundary be- 
tween the insulating and superconducting states for ar- 
rays with a str onger coupling q > 1 . Accordin g to 
the experiments (jJaeeer et a^J . ll989HOrr e^"a/lll986| ) the 
superconductor-insulator transition is controlled by the 
resistance of the array only and the critical resistance 
(in the case of granular films) is close to the resistance 
quantum Ro ~ QAkft. The Coulomb interaction drops 
out from the result. 

This interesting observation can be understood using 
Eqs. (|O^ . IT?7irr^rn?|l . it is important to notice that 
the energy Eo (g) depends on the dimensionality. Using 
Eqs. (|3.23l 13.271) we see that in 3d arrays of the grains 
the function Eq (g) equals 

Eo (g) = cA/g, (3.32) 

where c is a number of the order unity, while in 2c? we 
obtain with the logarithmic accuracy 

2A 

Eo{g) = -^\n{gE,/A). (3.33) 

One can see from Eqs. (|3.32|l and H3.33|l that the effective 
charging energy is independent on the originally strong 
Coulomb interaction in 3d case and is almost indepen- 
dent on it in the 2d case. Substituting Eqs. H3.32|l and 
(I3.33|) into Eqs. (|3.29|l . we come to the conclusion that 
in the limit Ec 3> A the superconductor-insulator tran- 
sition occurs at gc ^ 1, which is in agreement the with 
experiments. 

Co nsiderations presented in the original 
work l)Chakravartv et all IT987|) lead to results qualita- 
tively similar to those presented here. In particular, the 




FIG. 17 The sketch of the phase diagram of a granular super- 
conductor in coordinates critical temperature vs. tunnelling 
conductance for two cases Ec <C Ao and Ec S> Aq. At large 
tunnelling conductances gr 1 the Coulomb interaction is 
screened due to electron intergrain tunnelling and the transi- 
tion temperature is approximately given by the single grain 
BCS value. In the opposite case g <C 1 the boundary be- 
tween the insulating and superconducting states is obtained 
comparing the Josephson Ej = gA and Coulomb Ec energies. 

phase boundary at T ^ in two dimensional case was 
found to be 

a^ = ^ln(l + ^a.^), (3.34) 

where ac is the critical value of the dimensionless con- 
ductance related to our conductance g as a = ng. 
T he phase boundary following from the calculations 
of lChakravartv et a^J |l987i) is represented in Fig.lT^l 

Summarizing our discussion we present in Fig. 1171 the 
phase diagram for a granular superconductor in the tem- 
perature vs. conductance coordinates for two cases of 
i?c < A and Ec > A. 

The strong coupling regime, g > I, can also be ana- 
lyzed using the diagrammatic technique. This approach 
is advantageous in many respects since it allows one 
to find the corrections to the critical temperature in a 
straightforward and rigorous way not relying on the mean 
field approximation used in the treatment of the effec- 
tive functional. This approach is considered in detail in 
Sec HlLDl 

We note that while the mean field theory gives quali- 
tatively correct results, it certainly cannot be viewed as 
a complete description of the model, Eq. H3.17|l because 
fluctuations near the mean field solutions can be impor- 
tant. In particular, Eq. (|3.2t)|) assumes that the system 
is not affected by the Josephson couplings so long as the 
mean value of ip is zero. However, the Josephson cou- 
pling is known to be important even in the insulating 
state. In the simplest case of the insulating state with 
the low Josephson coupling the effect of the later can be 
included via the perturbation theory. As a result, the 
charging energy is reduced, Er, En = En + SEc with 
the correction SEc given by llSachdevll20nl 

SEc = -zJ/2. (3.35) 



36 

Closer to the superconductor - insulator transition the 
charging energy gets renormalized even more strongly 
and in the vicinit y of the transitio n it is expected to have 
the scaling form, llSachdevll200l 

Eo{g)^\g:-g\\ (3.36) 

where 7 > is the critical index. One may interpret the 
energy Eq (g) as the increase of charging energy when 
adding a Cooper pair to a neutral grain. Exactly at 
the phase transition this energy vanishes according to 
Eq. (|3.36l) and this ensures the continuity of the Coulomb 
gap at the phase transition, since in the superconducting 
state this gap has to be zero. 

The Coulomb gap is an important characteristic of the 
insulating state. In particular, it controls the low tem- 
perature conductivity that has an activation form 

Experimentally measured conductivity of granular sam- 
ples is, however, usually more complicated than that 
given by Eq. H3.37|l and typically has a variable-range 
hopping type of behavior. Such a behavior can be ob- 
tained taking into account the electrostatic disorder ne- 
glected in the model, Eq. (|3.17() . 

C. Upper critical field of a granular superconductor 

As we already mentioned, the granular superconduc- 
tors have rather unusual magnetic properties. In partic- 
ular, at low temperatures, depending on the tunneling 
conductance g, the critical magnetic field is determined 
by either bulk formula H3.12|l a,t g > g* or by the critical 
magnetic field of a single grain 13.11|l a.t g < g*, where 
the characteristic conductance g* is defined by Eq. (|3.9I 
) . In this section we will discuss in detail the dependence 
of the upper critical field of granular superconduc- 
tors on the tunneling conductance g at arbitrary tem- 
peratures. We will consider granular samples that are 
relatively good metals in their normal state, such that 
g ^ 1. This will allow us to neglect the effect of the 
suppression of the critical temperature by the Coulomb 
interaction and fluctuations. 



1. Critical field of a single grain 

We begin with a discussion of the critical magnetic 
field of a sin gle grain . This problem was first considered 
long ago bv iLarkinl l)l965(l by means of averaging the 
Gor'kov equations over disorder. Another way to con- 
sider this problem is based on using a l ater technique , 
namely, a semiclassical Usadel equation, l|UsadelL ll97Cl|) . 
This way is more convenient for our purposes since it 
can easily be generalized to the case of a granular array. 
Assuming that the transition between the superconduct- 
ing and normal states is of the second order we define 
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the critic al field as a field at wh ich the self-consistency 
equation l)Abrikosov et all Il965|) 

A(r) = Ao^r^/^(r), (3.38) 

where lo = ttT (2n + 1) , acquires a nontrivial solution un- 
der decrease of the magnetic field. The function /^^ (r) is 
a quasiclassical Green function that can be obtained from 
the anomalous Gor'kov Green function F by integration 
over ^, Eqs. H3.2|) and H3.3|l . The Green function /^^{r) 
can be determined from the Usadel equation th at can be 
linear ized close to the phase transition, see e.g. ijKoDninl 

l2nnl^ 

(|cj| + D{-iW - 2eA/c)V2) /^(r) ^ A(r), (3.39) 

where D is the classical diffusion coefficient in the su- 
perconducting grain and A is the vector potential. We 
choose the vector potential corresponding to a homoge- 
neous magnetic field H in the form A (r) = [H x r]/2, 
which allows us to reduce Eq. H3.39|l to the form 

- Z?VV2 + 2i^(eA/c)2)/^(r) - A(r). (3.40) 

Equation H3.4()(l can be solved via perturbative expansion 
in A^. In the limit determined by Eq. I|3.16|l . one can 
neglect dependence of the functions /^^(r) and A(r) on 
coordinates. In the main approximation we can simply 
replace the term with A^ (r) in Eq. H3.4()|l by its average 
over the volume of the grain. Thus, we obtain the simple 
algebraic equation relating constant components of / and 
A 

/^(|u;|+a) = A. (3.41) 

where the depairing parameter a is the case of a spherical 
grain with the radius R is 

a = R^D{eH/cf /b. (3.42) 

Substituting the solution of Eq. (I3.41|l into Eg. II3.38|I we 
come to the standard equation (jAbrikosov and Gor'kovi 
EgGl^) relating the critical temperature to the depairing 
parameter a, Eq. (|3.42|) . 

ln(r,/r,o) = ?A(l/2) - V(l/2 + a/2^T,). (3.43) 

Here il}{x) is the digamma function and Tco = Tc{H — 0) 
is the single grain superconducting critical temperature 
in the absence of the magnetic field. One can see that at 
zero temperature the superconductivity is destroyed at 

aco = Ao/2, (3.44) 

where Aq is the superconducting order parameter at 
T,H — that is related in the standard way to the criti- 
cal temperature Tcq in the absence of the magnetic field 
as Aq = 7rTco/7 with \nj fa 0.577 being the Euler con- 
stant. For a single grain of a spherical shape we obtain 
Eq. H3.11|l for the critical field i/f'' at zero temperature. 




FIG. 18 Phase diagram of a superconducting grain at T = 
in coordinates Zeeman energy h vs. pairbreaking parameter 
a. The dashed line separates the gapless (GS) and gapful (S) 
superconducting states. First and second order phase tran- 
sition lines are shown by thick and thin solid lines, respec- 
tively (Bcloborodov et alll200^ . 

The above approach can be easily generalized to take 
into account the Zeeman splitting by making the shift 
|cli| \lo\ + hsignuj, where h is the Zeeman energy h — 
^ibH. Then, Eq. (j^^t^ is modified to 

ln(T,/T,o) = ^(1/2) - Re [ ^(1/2 + {a + ih)/2^T,)] . 

(3.45) 

We note that Eq. (|3.45|) is valid as long as the phase tran- 
sition is of the second order. However, in the li mit where 
the Zeeman effect dominates the orbital one ICloestonl 
(1960) the phase transition is known to be first order. 
Analysis of the first order phase transition at finite tem- 
peratures is rather involved and here we present only 
the corresponding phase diagram at zero temperature 
fFig. ll8|) that illustrates how the seco nd order phase tran- 
sition turns into the first order one ijBeloborodov et aZI . 
l2006j) . To simplify the discussion we assume that the 
Clogston limit is not reached and Eq. 13.45|l describing 
the second order phase transition is applicable. 

At zero temperature, Eq. H3.45|l reduces 
to ()Beloborodov et all l2000|) 

+ = A2/4, (3.46) 

such that Eq. H3.11|l for the critical field fff in the pres- 
ence of the Zeeman coupling can simply be generalized 
via the substitution 

Ao ^ Ao = y^A2-(2/i)2. (3.47) 

2. Critical magnetic field of a granular sample 

According to subsection IIII. G . T1 the effect of the mag- 
netic field on a single grain reduces to the renormalization 
of the Matsubara frequency 

to ^ uj = Lu + ih + a. (3.48) 

So, considering the critical temperature of the whole 
granular sample one can take into account the local single 
grain effects via the frequency renormalization, 13.48|l . 
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The critical temperature of the granular array can be 
again found from the linearized Usadel equation that can 
be written for the granular sample in the Fourier repre- 
sentation as 



(gV2)A(q-2eA/c) /(q) = A(q: 



(3.49) 



where A(q — 2eA/c) = 2 X]a(l ^ cos[(q — 2eA/c)a]) and 
q is the quasimomentum of the periodic lattice. In the 
regime under consideration, Eq. (|3.16|l , the magnetic flux 
per a unite cell of the lattice is always small. This allows 
us to find the lowest eigenvalue of the operator A(q — 
2eA/c) as the lowest Landau level of the operator for the 
corresponding continuous system a^(q— 2eA/c)^, which 
gives 2eHo? /c. 

Thus, the evaluation of the orbital effect that originates 
from flux accumulation on scales involving many grains is 
very similar to the analogous effect in the homogeneously 
disordered samples, which is a direct consequence of the 
smallness of the magnetic flux per unite cell. Now one can 
see that the orbital effect that comes from large distances 
simply gives an additional contribution gSa^eH/c to the 
depairing parameter a. Using Eq. 1)2. 3|l for the effective 
diffusive coefficient D^f / we obtain for the total depairing 
parameter including both single grain and bulk orbital 
effects of the magnetic field 



a = R^D{eH/cf/b + Def/eH/c. 



(3.50) 



We note that although Eq. (|3.50() has been obtained 
within the periodic cubic lattice array model it is in fact 
more general. Indeed, the first term in the right hand 
side of Eq. (|3.5UI) represents the single grain effect and 
is sensitive to the grain size only. The second term, be- 
ing expressed in terms of the effective diffusive coefficient 
Deff, is not sensitive to the structure of the array at all. 

Comparing the two contributions to the depairing pa- 
rameter a in Eq. (|3.5U|) we obtain immediately the value 
H* characterizing the crossover from the single grain to 
the bulk orbital effects 



H* 



5 00 Deff 



(3.51) 



Equations (|3.5U|) and H3.45|l determine the superconduct- 
ing transition line in the coordinates T vs. H plane. 

In the limit of zero temperature T ^ 0, Eq. H3.45|l re- 
duces to Eq. (|3.46|l . which allows us to find a simple ex- 
pression for the critical filed at T = resolving Eq. (|3.5U|) 
with respect to H : 



Hc{T = 0) 



5(/>o 



2I?7ri?2 V 



JJ2 , 2Aoi?2J 



D 



eff 



(3.52) 

where Ag defined by Eq. H3.47|l takes into account the 
Zeeman effect. 



D. Suppression of the superconducting critical temperature 

The BCS theory gives a very accurate description of 
superconductors in bulk well conducting metals. How- 
ever, in strongly disordered or granular metals there can 
be considerable deviations from this theory. This con- 
cerns, in particular, the superconducting transition tem- 
perature. 

In this subsection we consider suppression of the super- 
conductivity in granular metals at large tunneling con- 
ductances between the grains, g ^ 1. In this limit, all 
properties of the system should be close to those following 
from the BCS theory. The mean field transition tempera- 
ture Tc of the granular superconductors has to be close to 
that of the corresponding disordered bulk superconduc- 
tor. A difference between the granular and disordered 
bulk superconductors may appear only in corrections to 
the BCS theory. These corrections will be discussed in 
the present subsection. 

The main mechanisms of the suppression of the su- 
perconductivity are the Coulomb repulsion and fluctua- 
tions of the superconducting order parameter. For ex- 
ample, disorder shifts significantly the sup erconducting 
trans it ion te m perature in the 2d thin fi l ms llFinkelstein , 



transition tem perature m tne za tnm mr ns inKeistein 
19871 Il994t llshida and Ikeda'. Il998t iLarkinL |199C 



Maekawa et al\ | l983t iMackaw a and FukTivamaL 1198- 
Ovchinnikovlll973(l . The physical reason for the suppres- 



sion of the critical temperature is that in thin films the 
interaction amplitude in the superconducting channel de- 
creases due to peculiar disorder-induced interference ef- 
fects that enhance the effective Coulomb interaction. We 
refer to this mechanism of the superconductivity suppres- 
sion as to the fermionic mechanism. 

The superconducting transition temperature can also 
be reduced by the fluctuations of the order parameter, 
the effect being especially strong in low dimensions. The 
corresponding mechanism of the superconductivity sup- 
pression is called the bosonic mechanism. In particular, 
the bosonic mechanism can lead to the appearance of 
the insulating state at zero temperature (Efetov. 1980t 
iFisheii [T99n[ ISimanekL I1994D . 

In this subsection we discuss the suppression of the 
superconductivity in terms of the perturbation theory. 
Although this method is restrictive and cannot be used 
to study the superconductor -insulator transition, it is 
useful in a sense that both the relevant mechanisms of 
the reduction of the critical temperature can be included 
systematically within the same framework. The useful- 
ness of the perturbative calculations has been demon- 
strated in the study of properties of normal granular met- 
als in Sec. III.DI where the new important energy scale 
r, Eq. H2.2I) was introduced. In view of these findings 
one may expect that the correction to the superconduct- 
ing transition temperature can be different depending on 
whether the temperature is larger or smaller than the 
energy scale F. 

We start with the following Hamiltonian Hs describing 
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Ztl — c 



b) 



d) 



FIG. 19 Diagrams a) - c) describe the correction to the su- 
perconducting transition temperature due to Coulomb re- 
pulsion. The diagram d) describes correction to the transi- 
tion temperature due to superconducting fluctuations. All 
diagrams are shown before averaging over the impurities. 
The solid lines denote the electron propagators, the dashed 
lines describe screened Coulomb interaction and the wavy 
lines describe the propagator of superconducting fluctua- 
tions (,Bcloborodov et al... 2005a.'l . 



the granular system 



(3.53) 



where the Hamiltonian H for the normal granular metals 
is given by Eq. (|2.15|l and He-ph is an additional electron- 
phonon interaction on each grain 



i.k.k' 



(3.54) 



where i labels the grains, k e 
is the interaction constant 



(k,T), -fc= (-k,i); A> 
afi^itti^k) are the creation 
(annihilation) operators for an electron in the state k of 
the i'th grain. 

The superconducting transition temperature can be 
found by considering the anomalous Green function 
F i n the presence of an infinitesimal source of pairs 
A, llOvchinnikovill973D . Neglecting fluctuations and in- 
teraction effects th e anomalous Green fun ction F is given 
by the expression l|Abrikosov et all ll96,'J) 



F{^,io) = A/{oj^+e), 



(3.55) 



where ^ = p^/2m — /i, and ui — 27rT(n + 1/2) is the 
fermionic Matsubara frequency. The suppression of the 
transition temperature Tc is determined by the correction 
to the function F(^, w) 



AT, 
T, 



T 
A 



J d?5]<5F(e,u;), 



(3.56) 



where the function 6F{^,lu) represents the leading order 
corrections to the anomalous Green function F(^, w) due 
to fluctuations of the order parameter and Coulomb in- 
teraction. All diagrams (before impurity averaging) con- 
tributing to the suppression of the transition tempera- 
ture in Eq. H3.56|l are shown in Fig. [201 One can see 



that there exist two qualitatively different classes of dia- 
grams. First, the diagrams a) - c) describe corrections to 
the transition temperature due to Coulomb repulsion and 
correspond to the so called fermionic mechanism of the 
suppression of superconductivity. The second type, dia- 
gram d) , describes a correction to the transition temper- 
ature due to the superconducting fluctuations and repre- 
sents the bosonic mec hanism. The details of calcu lations 
are presented in Ref . ijBeloborodov et all . l2005aj) . They 
are rather straightforward and therefore we present and 
discuss below only the final results for the suppression of 
Tc. 

It is convenient to separate corrections due to the 
bosonic and fermionic mechanisms from each other and 
write the result for the suppression AT, of the supercon- 
ductor transition temperature as 



AT, 



/AT, 



/AT, 



V Tc 



(3.57) 



/ 



where the two terms in the right hand side correspond to 
the bosonic and fermionic mechanisms, respectively. The 
critical temperature T, in Eq. ()3.57fl is the BCS critical 
temperature. 

It is shown bv lBeloborodov et all l|2005al) that at high 
transition temperatures, T, > F, the fermionic correction 
to the superconducting transition temperature does not 
depend on the dimensionality of the sample and takes 
the form 



/AT, 



V Tc 



f 



2,3. 



(3.58) 



where ci — 7C(3)/27r^ — (ln2)/4 is the numerical coef- 
ficient and d is the dimensionality of the array of the 
grains. 

In contrast, in the low temperature regime, T, < F, the 
fermionic mechanism correction to the superconducting 
transition temperature depe nds on the dimensionality of 
the sample and is given by l)Beloborodov et all l2005a|) 



AT, 

Tc 



5^1n2f , 

Z-K q Ic 

1 . 1^3 r_ 



d = 3 
d = 2 



(3.59) 



where A = a'^ J d^q/{27T)^X~^ w 0.253 is the dimension- 
less constant, a is the size of the grain and Aq is given by 
Eq. (|2.3U|) . Note that in the low temperature regime, 
Tc < F, the correction to the critical temperature in 
two dimensions agrees with the one obtained for homo- 
geneously disordered superconducting films provided the 
substitution F — > is done. 

The correction to the transition temperature due to 
the bosonic mechanism in Eq. (|3.57|l is qualitatively dif- 
ferent from the fermionic one. It remains the same 
in both the regimes T < F and T > F and is given 
by l|Beloborodov et all l2005a|) 



AT, 
Tc 




d = 3 

d = 2 



(3.60) 
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where Ci^) is the zeta-function and the dimensionless 
constant A is written below Eq. H3.59|l . 

Note that the energy scale F does not appear in this 
bosonic part of the suppression of superconducting tem- 
perature in Eq. H3.60() . This stems from the fact that the 
characteristic length scale for the bosonic mechanism is 
the coherence length ^eff that is assumed to be much 
larger than the size of the single grain. The result for 
the two dimensional case in Eq. H3.60|l is written with a 
logarithmic accuracy assuming that \n{g^6/Tc) ^ 1. 

We see from Eqs. (|3.59ll and (|3.60|) that the suppres- 
sion of the superconducting transition temperature be- 
comes stronger with diminishing the coupling g between 
the grains. Of course, the calculations may be justified 
only when the correction to the transition temperature 
is much smaller than the temperature itself. 

The above expression for the correction to the tran- 
sition temperature due to the bosonic mechanism was 
obtained in the lowest order in the propagator of super- 
conducting fluctuations and holds therefore as long as the 
value for the critical temperature shift given by Eq. (|3.6()(l 
is larger than the Ginzburg region (AT)g 

f d-3 
(AT),^ (3.61) 

I g 

Comparing the correction to the transition temperature 
Tc given by Eq. (|3.60() with the width of the Ginzburg 
region, Eq. (|3.61(l . one concludes that for 3d granular 
metals the perturbative result H3.60|l holds if 

Tc < g^S. (3.62) 

In two dimensions the correction to the transition tem- 
perature in Eq. (|3.60(l is only logarithmically larger than 
(AT)g in Eq. (|XHT|) . The two dimensional result 
holds with the logarithmic accuracy in the same temper- 
ature interval (|3.62(l as for the three dimensional samples. 

Note that inside the Ginzburg region higher order fluc- 
tuation corrections become important. Moreover, the 
non perturbative contributions that appear, in particu- 
lar, due to superconducting vortices should be taken into 
account as well. These effects destroy the superconduct- 
ing long range order and lead to Berezinskii-Kosterlitz- 
Thouless transition in 2d systems. The suppression of 
the transition temperature in the limit T > g^S in 3d 
should be studied considering strong critical fluctuations 
within the Ginzburg-Landau free energy functional. 

It follows from Eq. that in the limit T > T 

the fermionic mechanism of the suppression of the su- 
perconductivity is no longer efficient. This can be seen 
rather easily in another way using the phase approach 
presented in Sec. III. After the decoupling of Coulomb 
term (|2.39|l by the integration over the auxiliary field 
V and gauge transformation (|2.42|) the phase enters the 
tunneling term, Eq. (|2.53() only. However, this term is 
not important in the limit T 3> F and we conclude that 
the long range part of the Coulomb interaction leading 
to charging of the grains is completely removed in this 



limit. Therefore, the effect of the Coulomb interaction 
on the superconducting transition temperature must be 
small and this is seen from Eq. (|3.58(l . This conclusion 
matches well the fact that the upper limit in the loga- 
rithms in Eq. (|3.59() is just F and, at temperatures ex- 
ceeding this energy, the logarithms should disappear. 

At lower temperatures T < F the phase description 
is not applicable and, as a consequence, we obtain the 
non-trivial result, Eq. H3.59|l . This result is of the pure 
quantum origin, and the interference effects are very im- 
portant for its derivation (one should consider a contri- 
bution of the diffusion modes). In the high temperature 
limit T 3> F the interference effects are suppressed and 
this is the reason why the fermionic mechanism of the 
suppression of the superconductivity is no longer efficient. 

The results obtained, Eqs. H3.58|l - (|3.62|) . suggest an 
experimental method of extracting the information about 
the morphology of the samples from the Tc data. Indeed, 
one can study the dependence of the superconducting 
transition temperature Tc of the granular metals as a 
function of the tunneling conductance g by comparing 
several granular samples with different tunneling conduc- 
tances (different oxidation coating). The experimental 
curves for the suppression of Tc should have a different 
slope at high Tc > F and low Tc < F critical temperatures 
due to the fact that the suppression of the superconduc- 
tivity is given by the two different mechanisms. The in- 
formation on the morphology of the sample, i.e. whether 
the samples are homogeneously disordered or granular 
is then obtained from the dependencies of the critical 
temperatures on the tunneling conductance Tc{g). (We 
remind the reader that the scale F exists in the granular 
samples only). 

Another consequence of Eqs. H3.58|l - (|3.62|l is the fol- 
lowing: since at low critical temperatures Tc <C F the 
suppression of the superconductivity in granular metals 
is provided by the fermionic mechanism and coincides, 
upon the substitution F r~^, with the proper result 
for h omogeneously disordered samples, ijOvchinnikovl 
I1973D . one can gen e ralize the renormalization group re- 
sult by iFinkelsteinl l)l987j) for the Tc suppression. The 
latter result obtained for homogeneously disordered films 
can be directly applied to the granular superconductors 
by making the proper substitution for the diffusion coef- 
ficient D = Fa^, where a is the size of the single grain. 



E. Magnetoresistance of granular superconductors 

In this subsection we discuss magnetoresistance of a 
granular superconductor in a strong magnetic field. We 
start our consideration with a discussion of the experi- 
ment (Gcrbcr et at, 199^ where transport properties a 
system of Al superconducting grains in a strong magnetic 
field were studied. The samples were three dimensional 
and quite homogeneous with a typical grain diameter 
120 ± 20 A. Bulk superconductivity could be destroyed 
by application the strong magnetic field leading to the 
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appearance of a finite resistivity. Magnetic field above 
17T was enough to destroy even the superconducting gap 
within each grain. 

The depe ndence of the r e sistivi ty on the magnetic field 
observed in iGerber et al\ l|l997t) was not simple. Al- 
though at extremely strong fields the resistivity was al- 
most independent of the field, it increased when de- 
creasing the magnetic field. Only at sufficiently weak 
magnetic fields the resistivity started to decrease and fi- 
nally the samples displayed superconducting properties. 
A similar behavior had been reported in a number of pub- 
licati ons ijGantmakher et aIl . ll996HParthasarathv et all 

A negative magnetoresistance due to weak lo- 
caliza tion effects is n ot uii us ual in disordered met- 
als l)Altshuler et all 119801 iLee and Ramakrishnanl 
Ho wever, the characte ristic magnetic field in 
experiment ijGerber et all Il997j) was several tesla such 
that all weak localization effects had to be strongly 
suppressed. 

We want to present an explanation for this unusual be- 
havior based on consideration of superconducting fiuctu- 
ations. At first glance this idea looks counterintuitive, be- 
cause, naively, superconducting fiuctuations are expected 
to increase the conductivity. Nevertheless, it turns out 
that the magnetoresistance of a good granulated metal 
{g ^ 1) in a strong magnetic field, H > He, and at 
low temperature, T <C TL, is negative. In our model, 
the superconducting gap in each granule is assumed to 
be suppressed by the strong magnetic field. All the in- 
teresting behavior considered below originates from the 
superconducting fiuctuations that lead to a suppression 
of the density of states (DOS) and decrease the conduc- 
tivity. 

Theory of superconducting fluctuations near the tran- 
sition int o the supercondu cting state has been developed 
long a go lAbrah ams a7l.ll970; .Aslamazov and Larkinl 
119681: iMakl Il968albl: iThompsonL ll97f]D ffor a review 
see iLarkiiTanT^/arlam l)2005|) ') Above the transition 
temperature Tc, non-equilibrium Cooper pairs are formed 
and a new channel of ch arge transfer opens (Asl a mazov - 
Larkin contribution) ijAslamazov and Larkinl Il968|) . 
This correction always gives positive contribution to the 
conductivity. Another fluctuation contribution comes 
from a coherent scattering of the electrons forming a 
Coop er pair on impurities (Maki-Tho mpson contribu- 
tion) l|Maki ll968aliol: iThompsonl Il970|) . This correction 
also results in a positive contribution to conductivity, 
though, in principle, the sign of it is not prescribed. 

Formation of the non-equilibrium Cooper pairs re- 
sults a lso in a fluctuational ga p in the one-electron spec- 
trum ijAbrahams et all Il970|) but in conventional (non 
granular) superconductors the flrst two mechanisms are 
more important near Tc and the conductivity increases 
when approaching the transition. 

The total conductivity for a bulk sample above the 
transition temperature Tc can be written in the following 



form 



5a, 



(3.63) 



where (JDrude = {e^Tn)/m is the conductivity of a nor- 
mal metal without electron-electron interaction, r is the 
elastic mean free time, m and n are the effective mass 
and the density of electrons, respectively. In Eq. I|3.63|l . 
5a is a correction to the conductivity due to the Cooper 
pair fluctuations 



5a = 5a DOS + 5a al + 5a mt, 



(3.64) 



where 5aDOS is the correction to the conductivity due to 
the reduction of the DOS and 5aAL and 5aMT stand for 
the Aslamazov-Larkin (AL) and Maki-Thompson (MT) 
contributions to the conductivity. Close to the critical 
temperature Tc the AL correction is more important than 
both the MT and DOS c orrections and its contr i bution 
can be written as follows l)Aslamazov and Larkinl Il968|) 



5aAL 

aOrude 



Tc 



T-Tc 



(3.65) 



where 7 is a small dimensionless positive parameter, 
7 <C 1, depending on dimensionality and /3 = 1/2 for the 
three-dimensional case (3d), 1 for Id and 3/2 for quasi- 
\d. Equation l|3.65|l is valid for 5a alI aomde ^ 1. 

For quite a long time the superconducting fluctuations 
have been studied near the critical temperature Tc in a 
zero or a weak magnetic field. Only recently the limit of 
very low temperature, T <C T^, and a strong magnetic 
field, H > He was considered (Beloborodov and Efetovl 
IBeloborodov et all I200C ) and some qualitatively 



neiu. 
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new effects were predicted. In this subsection we concen- 
trate on this limit. 

As in conventional bulk superconductors, we can write 
corrections to the classical conductivity ao, Eq. 1)2. l|l . as a 
sum of corrections to the DOS, Aslamazov-Larkin (AL) , 
aAL, and Maki-Thompson (MT), auT, corrections. Di- 
agrams describing these contributions are represented 
in Fig. 1201 They have the conventional form but the 
wavy lines standing for the superconducting fiuctuations 
should be written taking into account the granular struc- 
ture. It turns out that in the regime T <C Tc and H > He 
the DOS corrections plays a very important role: This 
correction reaches its maximum at H He, where He is 
the field destroying the superconducting gap in the single 
grain. At zero temperature and close to the critical field 
He such that T/Aq > h , where h = {H - He) /He, the 
m aximum value 5anns for 3rf granular superconductors 
is (IBeloborodov et aiJ . l2000|) 



5aDos 1 r / Aq 

ao 3g Ao V ^ 



(3.66) 



We see from Eq. (|3.66|) that the correction to the con- 
ductivity 5a DOS i) is negative and its absolute value de- 
creases when the magnetic field increases, ii) it is smaller 
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FIG. 20 Diagram a) describes correction to DOS, diagrams b) 
and c) describe corrections to conductivity due to supercon- 
ducting fluctuations. The wavy lines denote the propagator 
of the fluc tuations, the dashed hnes stan d for the impurity 
scattering iBeloborodov and Efetovlll999t) . 



than unity and it becomes comparable with ctq when the 
tunneling conductance g ~ 1. However, such values of 
g mean that we would be in this case not far from the 
metal-insulator transition. Then, we would have to take 
into account all localization effects, which is not done 
here. 

Even in the limit of strong magnetic fields, 
H ^ Hr, the correct i on to cto can still be notice- 



able l)Beloborodov et all . \200(j(\ 
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(3.67) 



where, for spherical grain, Eq {H) = | (^^^) D with D 
being the diffusion coefficient in the single grain. Equa- 
tion l|3.67|l shows that in the region H ^ He the correc- 
tion to the conductivity behaves essentially as Sa^os ~ 
H^^ , which is a rather slow decay. 

We emphasize that the correction to the conductivity 
coming from the DOS, Eqs. H3.66|l and (|3.t)7|) remains 
finite in the limit T 0, thus indicating the existence of 
the virtual Cooper pairs even at T = 0. 

In order to calculate the entire conductivity we must 
investigate the AL and MT contributions (Figs. l2Qh and 
I20b'l. Near Tc these contributions are most important 
leading to an increase of the conductivity. At low tem- 
peratures T ^ Tc and strong magnetic fields H > He 
the situation is completely different. It turns out that 
both the AL and MT contributions vanish in the limit 
T ^ at all > i/c and thus, the correction to the 
conductivity comes from the DOS only. 

So, at low temperatures, estimating the total correc- 
tion to the classical conductivity uo, Eq. (|2.42(l . one 
may use the formulae H3.66|) and H3.67|l . We present 
the final result for AL correction shown in Fig. 120b at 
T <C Tc and h <C F/An for Sd granular superconduc- 
tors (^Bcloboro dov et a/J . |2000|) 
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It follows from Eq. H3.68|l that at low temperatures the 
AL correction to the conductivity is proportional to the 
square of the temperature 5a al ^ and vanishes in 
the hmit T ^ 0. 

For MT corrections shown in Fig. 120b one obtains in 
the limit T <C Tr and h <C F/An for 3d granular super- 
conductors ijBeloborodov et aTll2000(l 



(3.69) 



(3.68) 



The temperature and magnetic field dependence of 
SaAL and 5<jmt is rather complicated but they are def- 
initely positive. The competition between these correc- 
tions and a DOS determines the sign of the magnetoresis- 
tance. One can see from Eqs. (|3.68|l and H3.69|l that both 
the AL and MT contributions are proportional at low 
temperatures to T^. Therefore the crjjos in this limit is 
larger and the magnetoresistance is negative for all H^- 
In contrast, a,t T Tc and close to He, the AL and MT 
corrections can become larger than aoos resulting in a 
positive magnetoresistance in this region. Far from He 
the magnetoresistance is negative again. 

In the above considerations we did not take into ac- 
count interaction between the magnetic field and spins 
of the electrons. This approximation is justified if the 
size of the grains is not very small. Then, the or- 
bital magnetic critical field H"^ destroying the super- 
conducting gap is smaller than the paramagnetic limit 
IJ.Hc ~ Aq and the orbital mechanism dominates the 
magnetic field effect on the superconductivity. However, 
the Zeeman splitting leading to the destruction of the su- 
perconducting pairs can become important if one further 
decreases the size of the grains. The results for conductiv- 
ity corr ections including Ze eman splitting can be found 
in Ref. ijBeloborodov ad 12000). 

The results presented in this subsection show that, 
even if the superconducting gap in each granule is de- 
stroyed by the magnetic field, the virtual Cooper pairs 
can persist up to extremely strong magnetic fields. How- 
ever, the contribution of the Cooper pairs to transport 
is proportional at low temperatures to T^ and vanishes 
in the limit T 0. In contrast, they reduce the one- 
particle density of states in the grains even at T = 0, 
thus diminishing the macroscopic conductivity. The con- 
ductivity can reach its classical value only in extremely 
strong magnetic fields when all the virtual Cooper pairs 
do not exist anymore. This leads to the negative magne- 
toresistance. 

Qualitatively, the results for resistivity behavior of 
granular superconductors at low temperatures T <C Tc 
and strong magnetic fields H > He are summarized in 
Fig. 1^ The resistivity R at low temperatures grows 
monotonously when decreasing the magnetic field and 
it reaches asymptotically the value of the classical resis- 
tivity Rq only at extremely strong magnetic fields. The 
transition into the superconducting state occurs at a field 
He2 that may be lower than He 
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Magnetic field 

FIG. 21 Schematic picture for resistivity behavior of the gran- 
ulated superconductors as a function of the magnetic field 
at low temperatures T <^ Tc- The resistivity at low tem- 
peratures grows monotonously when decreasing the magnetic 
field. It reaches asymptotically the value of the classical re- 
sistivity Ro only at extremely strong magnetic fields. The 
trans ition into the supercond ucting state occurs at a field 
Hc^ llBeloborodov et aZ.ll2000ll . 



IV. DISCUSSION OF THE RESULTS 

A. Quantitative comparison with experiments 

The niodels of the normal and superconducting granu- 
lar electronic systems considered in this review are rele- 
vant to many experiments on different materials and we 
have presented quite a few of them. This is not accidental 
because we selected for the presentation such theoretical 
results that could be relevant for existing experiments. 

On the one hand, the review is not as complete as it 
could be because we did not include many interesting the- 
oretical works where such physical quantities, a s tunnel- 
ing density of states llB eloborod ov et a^J■l2004a^■ thermal 
conductivity fBcloborodov et all I2005bt iBiagini et all 
l2b05b: .Loh and Tripathi. 2006), etc. were considered. 
The only reason for this omission is that at present we 
are not aware of clear measurements of these quantities 
in the granular materials. 

On the other hand, this gave us the opportunity to 
concentrate on experimentally well studied phenomena. 
Some of them, like the logarithmic temperature depen- 
dence of the conductivity of well conducting granular 
metals or the hopping like temperature dependence of 
the insulating very regular granular systems remained 
unexplained for several decades. As we have seen, these 
dependencies and many interesting features of the exper- 
iments can be explained using rather simple models of 
the metallic grains. The simplicity of the models and 
a small number of parameters characterizing the system 
suggests that there should be not only a good qualitative 
but also a quantitative agreement between the theory and 
experiment. 

So, we present now a quantitative comparison of sev- 
eral theoretical results with experimental data. Of 
course, we cannot discuss here all existing measurements 
and therefore only three different experiments are con- 
sidered. 



1. Logarithmic temperature dependence of the conductivity 

In this subsection we will compare with experiments 
the logarithmic temperature dependence of the conduc- 
tivity, Eq. (|2.72|) . This formulae should be applicable at 
not very low temperature that simplifies the experimental 
check. Specially prepared granular metals like those con- 
sidered in Refs. [A beles et at, 1975: Gerber et ai, 199'3 
ISimon et aZlll987|) should be the first objects of the ap- 
plication of the theory developed. 

As an example, we compare now Eq . I|2 . 72l) with the ex- 
perimental results of lGerber et on films made 
of Al grains embedded in an amorphous Ge matrix. At 
low temperatures, superconductivity in Al grains was de- 
stroyed by a strong magnetic field. Depending on the 
coupling between the grains (extracted from the conduc- 
tivity at room temperatu res) the samples of the experi- 
ment l|Gerber et all Il997|) were macroscopically either in 
an insulating state with the temperature dependence of 
the resistivity R ~ exp (a/T^/^) or in a "metallic" one. 
However, the resistivity of the metallic state depended 
on temperature and the authors suggested a power law, 

^ T"^ with an exponent a 1. As the exponent a 
for the "metallic" sample was small we can argue that 
Eq. (|2.72l) should not be worse for fitting the experimen- 
tal data. Then, we can estimate the exponent a without 
using fitting parameters. 

The sample of the experiment ijGerber et all Il997j) 
had the room temperature resistivity Rq = 7.3 x lO'^ficm. 
The diameter of the grains was 120± 20A, which allows, 
using the value h/e^ = 4.1 x 10~^U, and the relation be- 
tween the conductivity and conductance for 3d samples 
a = 2e^g/ha, to estimate the dimensionless tunneling 
conductance as g = 0.34. Taking d = 3 and (z = 6) in 
Eq. (|2.72|l we obtain a = 0.15, which agrees fairly well 
with the experimental value a = 0.117. Keeping in mind 
that the theoretical dependence (|2.72|) was obtained for 
a periodic cubic array while in realistic samples the co- 
ordination number can be only approximately close to 6 
one can hardly hope for a better agreement. 

A logarithmic dependence of the resistivity on temper- 
ature has been observed in other granular materials. In 
Ref. l)Simon et all Il987(l a granular cermet consisting of 
Nb grains in a boron nitride insulating matrix was stud- 
ied. Again, at small coupling between the grains the tem- 
perature dependence of the resistivity exp (a/T^/^) was 
observed in a very broad interval of temperatures. The 
resistivity of samples with a strong coupling between the 
grains was very well described by the law 

i? = i?oln(To/r), (4.1) 

which is close to Eq. 1)2.72(1 provided the temperature in- 
terval is not very large, such that the variation of the 
resistivity is small. However, Eq. (|4.1|) gave a good de- 
scription for the temperature dependence of the resistiv- 
ity in a very broad region and the changing of the resis- 
tivity was not small. The reason for the applicability of 
Eq. H4.1|) in a so broad interval of temperatures is not 
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clear because according to the results of the renormal- 
ization group analysis of Sec. Ill.Cl not the resistivity but 
the conductivity should obey Eq. H4.1|l . A more careful 
experimental study might clarify this question. 

The unusual logarithmic behavior of the type, 
Eq. (|4.1|) . has been observed not only in "standard" 
granular systems but also in high- Tc cuprates at very 
strong magnetic fields. The first observation of this de- 
pend ence was done on u nderdoped La2^xSrxCu04 crys- 
tals llAndo et aljAim^ . The superconductivity in this 
experiment was suppressed with pulsed magnetic fields 
of 61T. It was found that both the in-plane resistivity 
Pab and out-of-plane resistivity pc diverged logarithmi- 
cally with decreasing the temperature. This means that 
a 3d effect was observed in a very strong magnetic field 
and traditional explanations for a logarithmic behavior 
like a weak localization or the Kondo effect could not 
clarify the situation. 

In a subsequent publication ijBoebinger et alj Il996j) 
a metal-insulator crossover was observed in the same 
material at a Sr concentration near optimum doping 
(x ~ 0.16). In underdoped samples both pab and pc 
showed no evidence of saturation at low temperatures 
and diverged as the logarithm of the temperature. The 
authors called this state "insulator" in contrast to the 
state at high doping where the resistivity did not have a 
pronounced dependence on the tem perature. It was con- 
jectured bv lBoebinger et all lll996() that the logarithmic 
behavior they obse rved might be related to the one seen 
in the experiment ijSimon et all . Il987|) on granular Nb. 
This would demand a phase segregation throughout the 
underdoped regime of LSCO. 

We hope that our results for the model of the 
granular materials may be applicable to the exper- 
imen t s on the La2-^ . Sr^Cu 04 crystals llAndo et all 
Il995t iBoebinger et all \l99(^ . which would mean that 
the underdoped crystals have a granular structure and 
the logarithmic behavior is due to the Coulomb in- 
teractio n. The transition to the met a llic s tate of 
Refs. |Ando et all Il995t IBoebinger et all Il996j) would 
mean that at higher doping the granularity disappears. 
A quanti tative comparison of Eq. I[2.72ll with the data 
of Refs. f Ando et al, 1995; Bocbin ger et all \l99(^ was 
done in .Beloborodov et al. (,2003.) and a good agreement 
was found. 

What is interesting, a microscopic granularity was 
directly experimentally observed in the supercon- 
ductin g state of Bi2Sr2CaCu20s+s by the STM 
probe l)Lang et a/.U200^ . As this cuprate is rat her simi- 
lar in its structure to those st udied in Refs. (.Ando et all 
Il995t IBoebinger et all Il996fl the assumption that the 
materials studied there were granular does not look 
groundless. 

The logarithmic dependence of the resistivity 
on temperature has also been obse rved in many 
other experiments. For example, in iGerbed l)l990() 
this dependence was observed in granular Pb 
films. It was also observed in phase compounds of 



Nd2-xCexCu04-y l|Radhakrishnan et adll990|) . 



2. Hopping conductivity 

As we have mentioned in the previous subsection, the 
logarithmic temperature dependence is usually seen in 
samples with a good coupling between the grains. If 
the coupling is weak, one observes usually the Efros- 
Shklovskii law, Eq. (f^T^ . 

According to the theory developed in Sec. III.GI the 
characteristic temperature Tq entering Eq. H2.14|l is de- 
termined by Eqs. H2.111|l and (|2.116() for the elastic and 
inelastic co-tunneling, respectively. 

Now we want to compare our results for the VRH con- 
ductivity of granular metals obtained in Sec. III.GI with 
the m ost recent experimental data. In Ref. (jTran et all 
1200.^ the zero-bias conductivity was studied for bilayers, 
trilayers, tetralayers and thick films of Au nanoparticles 
with the particle diameters around 5.5 nm and their dis- 
persion less than 5%. 

Figure El demonstrates that the zero-bias conductivity 
follows cr(T) - exp[-(ro/r)i/2] over the range 30-90K 
for the multilayers and 30— 150i4r for the thick films. The 
fits indicated by the lines give Tq = (4.00 ±0.02) x 10^ K 
and (3.00 ± 0.01) x lO^K for the multilayers and thick 
films, respectively. 

At temperatures slightly exceeding IQQK, the conduc- 
tivity for the multilayers starts deviating from the low- 
temperature behavior (dashed lines) and crosses over to 
Arrhenius behavior a{T) ~ exp[— iJ/ZcsT] (Fig.jHJ inset) 
with the activation energies U /ks ~ 320 ± . 

Associating this high-temperature behavior with the 
nearest-neighbor tunneling between the particles we 
can use U ~ 0.2Er., i n an alogy with the monolay- 
ers jParthasarathv et all 12004,) . This gives us an esti- 
mate Ec « 1600i^, where Ec = e^jka is the Coulomb 
charging energy of an individual grain, expressed in terms 
of the grain radius a, the electron charge e, and the dielec- 
tric constant of the surrounding medium k. The charging 
energy E^ ~ 1600i^ for this system leads to k « 4. 

The energy scale Tq is related to the localization length 
^ by Eqs. H2.111ll?.ll(j|l . Using the above values for k and 
To, we find Znm < ^ < 4nm for both the multilayers and 
thick films, which corresponds to the localization within 
one grain. This is an important check because larger val- 
ues would imply strongly coupled clusters of the grains. 

The typical hopping distance r* {T) is given by Eq. ( 
I2.122|l . The number of grains N* involved in a typical 
hop is N* — r* /do, with a center-to-center distance do ~ 
8nm between neighboring grains. At T = lOK this leads 
to A^* = 4 for multilayers and iV* = 4 — 5 for the thick 
films. 

As temperature increases, N* decreases down to N* ^ 
1 and this is the point of the crossover to a standard 
activation transport, Eq. (|2.13|) . For the multilayers, this 
happens at T w 90 — 95/^ but for the thick films only 
above T « 13QK. Both the estimates are in excellent 
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agreement with the data in Fig. |H| 



3. Negative magnetoresistance 

Let us compare the theoretical predictions on the neg- 
ative magnetoresistance presented in Sec. [ II I. El w ith th e 
available experimental results of iGerber et al\ l)l997(l . 
In that work three samples were studied. We will 
concentr ate our attention on the sample 2, Fig. 4, 
of Gerber et al that has a metallic conductivity 

behavior above the critical temperature. 

We analyze the case of low temperatures T <^ Tc and 
magnetic fields H > He, where Tc « 1.6K is the criti- 
cal temperature for Al grains studied in the experiment 
and He is the critical magnetic field that suppresses the 
superconductivity in a single grain. Eq(Hc) = Aq, where 
Eo{H) was defined below Eq. H3.67|l . At temperature 
T ~ Q.3K and magnetic field H ~ 4T this sample shows 
a large negative magnetoresistance. The resistivity has 
the maximum at H = 2.5T and the value of this peak 
is more than twice as large as the resistivity in the nor- 
mal state (that is, at H ^ He, when all superconducting 
fluctuations are completely suppressed). A negative mag- 
netoresistance due to weak localization (WL) is typical 
for disordered metals and in order to describe the exper- 
imental data, its value should be estimated along with 
the effects of the superconducting fluctuations discussed 
in the previous subsection. 

The total conductivity of the granular metal under con- 
sideration including effects of WL and superconducting 
fluctuations can be written in the form: 



Now, let us estimate the corrections Sa nps and 5awL 
using the parameters of the experiment ijGerber et all 
Il997|) . Fo r the typical diamet er 120 ± 20 A of Al grains 
studied bv lGerber et al\ ljl997|) . the mean level spacing S 
is approximately S « IK. Using the critical temperature 
Tc ~ 1.6K for Al we obtain for the BCS gap in a single 
grain the following result Aq ~ l.STc ~ 3K. Substituting 
the extracted values of the parameters into Eq. (|3.66(l we 
can estimate the maximal increase of the resistivity. As 
a result, we obtain ('5p//0o)max ~ ^-4' which is somewhat 
smaller but not far from the value {Sp/p)^^^ ~ 1 observed 
experimentally. 

Although Eq. H3.66|l gives smaller values of iSp/po)^^^. 
than the experimental ones, the discrepancy cannot be 
attributed to the weak localization effects. Using the 
experimental values of the conductance g, mean level 
spacing 6 and the superconducting gap Aq we find from 
Eq. H4.3|) that SawL is 10 times smaller than Sa^os- 
The value of the correction 5<7wl/c^o, Eq. (|4.3() . near He 
equals 2.8 x 10"^ 

At the same time, we should not take this disagreement 
too seriously because all the results were obtained un- 
der the assumption of a large conductance f/ 3> 1, while 
experimentally this parameter is not large. As the ex- 
perimental value of (5ct/cto is not small, we come again 
to the conclusion that the parameters of the system are 
such that the Eq. H3.65|l is no longer applicable and one 
should develop a more sophisticated theory to describe 
this region more accurately. 



B. Outlook 



0- = (To + 5<tdos + 5a AL + 5a 



MT 



5a 



WL- 



(4.2) 



At low temperatures, T <^ Tc, the contribution 5aoos 
originating from the reduction of DOS due to the forma- 
tion of the virtual Cooper pairs is larger than the con- 
tributions 5aAL and daMT since the latter vanish in the 
limit T — !■ 0. So, let us concentrate on estimating the 
contributions 5auos and 5awL- 

It is not difficult to show that in the case under con- 
sideration the weak localization corrections originating 
from a contribution of Cooperons are strongly suppressed 
by the magnetic field. Using app roximations developed 
bv lBeloborodov and Efetovl l)l999|) one can easily obtain 
the following expression for a 3d cubic lattice of metallic 
grains at tunneling conductances g <C Eo(H)/5 



5awL 



r 



EoiH) 



(4.3) 



Equation 1)4. 3|l shows that the weak localization correc- 
tion in the strong magnetic fields H > He considered here 
is always small. Comparing Eq. H4.3() with Eq. H3.66|l for 
DOS correction one can see that in the limit g <C Aq/iJ 
the contribution from the weak localization correction 
can be neglected. 



We discussed in this review a rather simple general 
model that allowed us to understand many properties 
of granular materials. This is a model of disordered or 
chaotic grains coupled to each other by the electron tun- 
neling. A very important ingredient of the model is the 
long range part of the Coulomb interaction taken in the 
form of the charging energy of the grains. For the de- 
scription of the superconducting grains we included the 
superconducting BCS gap into the consideration. 

We considered the limit of temperatures T exceeding 
the mean level spacing 5 in the grains and neglected quan- 
tum size effects. Considering the superconducting grains 
we assumed that the superconducting gap A was larger 
than 5. These regimes are most easily achieved experi- 
mentally and this was the main reason for choosing them. 

In spite of its simplicity, a variety of interesting phe- 
nomena has been derived from the model involved. Such 
phenomena as the logarithmic dependence of the con- 
ductivity on temperature or Efros-Shklovskii law for a 
regular system of grains with an almost perfect shape re- 
mained unexplained for several decades and have been 
clarified within the model only recently. Another in- 
teresting effect concerns the negative magnetoresistance 
due to superconducting fluctuations. For these effects we 
have obtained not only qualitative but a good quantita- 
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tive description. 

The granularity is a rather general phenomenon 
and it may be encountered rather unexpectedly. 
Thin metallic films are oft en rather granula r than 
homogeneously disorde r ed. llGantmakher et all Il996t 
iGoldman and Markovid Il998j) . Another unexpected 
conclusion about the granularity of some underdoped 
high T r cuprates has be en arrived at using STM tech- 
nique, I Lang e^'a/ll2002() . 

Studying the granular systems may help to under- 
stand some properties of homogeneously disordered met- 
als and superconductors. As all scales involved in the 
model we studied are much larger than the electron wave- 
length, it might in some cases simplify explicit calcula- 
tions. This concerns, in particular, study of the regime 
of 5 ~ 1. We presented arguments that the metal- 
insulator transition is possible at least in three dimen- 
sions. The superconductor-insulator transition is also 
possible. Studying these transitions is the most chal- 
lenging extension of the present study but this may be 
simpler than investigation of strongly disordered systems 
with interaction. 

From the experimental side, an evident extension of 
the works presented here is fabricating grains made of 
ferromagnets and studying properties of such systems. 
One can couple these grains directly or put them into 
normal metals or superconductors. One can make super- 
conducting grains and embed them into normal metals 
or not very strong ferromagnets, etc. All these systems 
promise very unusual properties. 

One of the examples of an unexpected behavior in the 
system of ferromagnetic grains imbedded in a supercon- 
ductor is inducing a magnetic moment in the supercon- 
ductor over large distances of the order of the size of the 
Cooper pairs. The direction of this magnetic moment 
is opposite to that of the ferromagnet and one co mes to 
a phenomen on of the "spin screening" l|Bergeret et al\ . 
120041 120051) . Different directions of the magnetic mo- 
ments of the grains may lead to t he phenomenon of an 
"odd triplet" superconductivity, l)Bergeret et all l200ll 
120051 

The list of new theoretical and experimental possibili- 
ties that are anticipated in granular systems can be con- 
tinued. Taking into account the growing number of exist- 
ing and potential industrial applications of the granular 
materials we are confident that all this is only the begin- 
ning of an exciting development. 
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APPENDIX A: Calculation of the tunneling probability P^i 
in the elastic regime, Eq. I2.105> 

In this appendix we derive the probability of elas- 
tic electron tunneling P^i between two distant grains 
through a chain of other grains, Eq. (|2.105|l . Such a 
probability can be found comparatively easily only for 
the model with the diagonal Coulomb interaction Efj = 
EiSij. For this reason we consider first such a simplified 
model and then discuss the generalization of the results 
obtained for the case of realistic capacitance matrices E![j 
with nonzero off-diagonal elements. 

In the model with the diagonal Coulomb interaction, 
the electron-hole excitation energies E^^ that will enter 
the final results for the hopping probability are given by 



S± = El 



(Al) 



where Vi is the local potential that models the presence 
of the electrostatic disorder. These energies have to be 
strictly positive and larger than the energies of the initial 
and final states otherwise the tunneling path could have 
been cut into two or more independent parts. This allows 
us to assume that the temperature is less than all E'^. 
For this reason we will consider the probability of the 
elastic process in the limit T ^ 0. 

The tunneling probability P between the sites iq and 
iN+i is proportional to the square of the absolute value 
of the amplitude Ai^^i^^-^ of the corresponding timncling 
process 



A 



=<0|c,„^,^c+|0> 



(A2) 



Here S is the evolution operator written in the interaction 
representation that takes into account only the tunneling 
part of the Hamiltonian (|2.15|) . The conjugate amplitude 
A* can be written as the probability of the inverse process 
of the tunneling between zat+i— st and zq— th grains 



A* — A- 



=< 0| Q„ S 6+ |0 > 



(A3) 



The probability of the tunneling process P = A* A 
can be found by perturbative expansion of the ampli- 
tude in the tunneling matrix elements and further 
averaging over the realizations of tij. At the same time 
the Coulomb interaction cannot be considered perturba- 
tively and should be taken into account exactly. In order 
to construct a proper perturbative expansion we use the 
gauge transformation described in Sec. lII.B"3l that allows 
us to transfer the strong Coulomb interaction into the 
phase factors that accompany the tunneling matrix ele- 
ments (see Eqs. (ITI^IT^ . 

Fluctuations of the phases ipi (r) are governed by the 
Coulomb action (|2.56|) and we obtain the following ex- 
pression for the correlation function 11 (t) that plays a 



47 




FIG. 22 This diagram represents the tunneUing probabil- 
ity via elastic co-tunneling processes. The crossed circles 
represent the tunnelling matrix elements t'ji! j,, e"*''^' e~^^^^'^^ 
where phase factors appear from the gauge transforma- 
tion. Wavy lines represent the average of the phase factors 
^gi0{Ti) g-^<t>{r2)'^ ^j^jj rcspcct to the Coulomb action. 



very important role for a description of the insulating 
state (c.f. Eq. |T7I )) 

n(Tl-T2) EE (exp(i0,(Tl) -i(/),(T2))) (A4) 

= exp {-El\Ti - T2I - Vdn - T2)) . 

Since we are interested in the optimal tunneling path 
we will consider only trajectories with no return points 
(no loops). This simplifies the consideration substantially 
because the phases on different sites are not correlated 
for the diagonal model under consideration. 

The gauge transformation approach allows us to rep- 
resent the tunneling probability P by the diagram shown 
in Fig. 1221 In order to simplify the derivation we work 
now using the basis of the exact eigenstates of the sin- 
gle particle Green functions that automatically takes into 
account the presence of the disorder within each grain. 
An alternative procedure based on the momentum rep- 
resentation leading to the same results would require the 
averaging over disorder in each grain and inclusion of the 
diffusion propagators. The Green function lines on the 
initial zq and final iN+i grains describe the processes for 
a particle located on these sites, while the intermediate 
Green function lines represent the tunneling amplitudes 
A (upper part) and A* (lower part). The wavy lines de- 
note the correlation function HA4|I that takes into account 
the Coulomb correlations. 

It may seem surprising that the two wavy lines that 
belong to the same i— th site are consider independently, 
i.e. that four exponent correlation function is factorized 
in a certain way into the two second order ones. One can 
check however that this factorization is justified indeed 
as long as the time intervals (ri_i,Ti) and {tI_^,tI) do 
not overlap. 

This factorization can most easily be understood by 
presenting the correlation function of N phase exponents 
as the exponent of the interaction energy of N charges in- 
teracting via a linear one dimensional Coulomb potential. 
The two nearby charges with the opposite signs create an 
"electric field" only in the region between the two. Thus, 
the energy of the dipoles can be considered independently 
as long as these dipoles do not overlap geometrically. For 
a detaile d discussion of the Coulomb analogy we refer the 
reader to IBeloborodov et al\ l)2005c^) . 



It follows from this discussion that all the electron 
Green functions in the diagram Fig. 1221 are accompanied 
by the correlation function n(r). For this reason, it is 
convenient to introduce a modified Green function G(t) 
in the following way 

G(r) = G{t)U{t). (A5) 

The Green function G{t) has a very simple form 
in the time representation: For positive arguments 
the Gree n function describes the electron excita- 
tion (^Abrikosov et a/.l[T965i') 

G^.(T>0) = KC,)-l]e-«'--^^^ (A6) 

where is the bare single particle electron energy 
counted with respect to the Fermi energy, Eq. H3.2|l and 
n{^i) is the Fermi distribution function. One can see that 
the Coulomb part of the electron excitation energy 
appears naturally in addition to the single particle en- 
ergy . For negative time arguments the Green function 
G describes the hole excitation and is given by 

Ge,(r<0)=n(e.)e~«'^+^'~^ (A7) 

where < is the bare electron energy counted with 
respect to the Fermi energy and is the Coulomb part 
of the hole excitation energy. 

Calculating the diagram in Fig.|221we further note that 
the intermediate time intervals are of the order of the in- 
verse Coulomb energy: {n+i - nl, \t-^-^ - t- \ E~''- . In 
the limit of strong Coulomb interaction {Ec ^ S) consid- 
ered here, we can integrate over the intermediate times 
independently. Each intermediate block gives rise to the 
factor 

Pk^T^ dd dT,G^{Ti) Q(t2) = (A8) 

where the energy Ek is the combination of the energies 
E^ and E~ defined by Eq. H2.107|) and Tk = gkSk, with 
gk being the conductance between fc— th and k + 1— st 
grains. We note that in the insulating system under con- 
sideration the energy scale F cannot be any more inter- 
preted as the escape rate, nevertheless its introduction is 
convenient even in this regime. 

In order to determine the time dependence of the tun- 
neling process we notice that the time intervals Tq — tq 
and T'pf^i — T7V+1 coincide within the accuracy of the in- 
verse charging energy. In the limit of the strong Coulomb 
interaction one can take both the intervals to be equal to 
the "instanton" time r that electron spends out of its 
original place at i = 0. Thus, the time dependence of the 
tunneling process is simply given by 

g-(e„+i-eo)r^ (A9) 

where Sn+i and Sq are the electron and hole excitation 
energies on the cites -I- 1-st and 0— th respectively. Mak- 
ing an analytical continuation to real times t = it and 
taking the last integral over t we obtain the delta function 

2TrS{eN+i - So), (AlO) 
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that indeed shows that the process is elastic, i.e. electron 
can tunnel only to the state with exactly the same energy. 

Finally, we obtain for the tunneling probability of the 
elastic process 

N 

Pel =w5{eN+i - eo)go JJ Pk-. (AH) 
fc=i 

where the factor w = n{(^o)[l — n{S^f44^i)\ takes into ac- 
count the occupation numbers of the initial and final 
states. 

We see that the total probability P^i contains the prod- 
uct of the ratios Ti/Ei. For this reason it is convenient 
to introduce the geometrical averages of theses quanti- 
ties along the tunneling path that were introduced in 
Sec. III. G. 21 Writing the total probability Pei in terms of 
the quantities F, E we obtain the result (|2.105|l . 

The factorization of the probability Pei, Eq. ljAll|l . 
into the product of the probabilities Pk means that the 
hops from grain to grain are independent from each 
other, which is a consequence of the diagonal form of 
the Coulomb interaction. 

If the Coulomb interaction is so long ranged that it in- 
volves many granules, the situation is more complicated 
since the integrals over the time variables cannot be taken 
on each site independently. Nevertheless, one can gener- 
alize the obtained results to the case of the long range 
Coulomb interaction as follows: The most important ef- 
fect of the off diagonal part of the interaction i?*-' is the 
renormalization of the excitation energies. Since the fi- 
nal result contains the single particle excitations only, 
this effect may be included by the proper definition of 
the electron-hole excitation energies 

Ef^Eu-tM^, (A12) 

where /x^ represents the local potential formed by both 
the external potential Vi and all charges surrounding the 
grain i. This procedure, however, does not include all 
the effects of the presence of the off diagonal part of E^? . 
The problem is that the virtual process represented by 
the diagram in Fig. does not correspond to a causal 
classical process. For this reason the virtual field created 
by an electron i at a time U affects in general the same 
electron when it is present on the other site at a different 
time. 

However, such an interaction decays at least as 
with the distance. This assumes that the long range ef- 
fects of the Coulomb interaction are not important, since 
the corresponding integral along the chain converges. 
Thus, all the effects that are not taken into account 
by the substitution (|A12|) are short range and may only 
lead to the change of a constant under the logarithm of 
the effective localization length (|2.109l) . Using t he anal- 
ogy w ith classical Coulomb problem ( Bcloborodo v et all . 

one can estimate the boundaries for the factor 
under the logarithm as 0.5 < c < 1.0. 




FIG. 23 This diagram represents the tunnelling probabil- 
ity via inelastic cotunneling processes. The crossed circles 
represent the tunnelling matrix elements i'^ y e"*'*-^' g-«'/>3{T) 
where phase factors appear from the gauge transforma- 
tion. Wavy lines represent the average of the phase factors 
(e^'^'^i) e"''*'^^') with respect to the Coulomb action. 

APPENDIX B: Calculation of the t unneling probability Pi„ 
in the inelastic regime, Eq. 12.1141 

In this appendix we derive the probability of the in- 
elastic cotunneling Pj„ through a chain of grains. As in 
the case of the elastic processes considered in the previ- 
ous section, the probability of the inelastic processes can 
easily be found for the model with the diagonal Coulomb 
interaction. For this model, we follow the same steps as 
previously, making the gauge transformation described 
in Sec. lII.B'ISl and expanding the tunneling probability in 
the tunneling matrix elements. 

The diagram describing the inelastic process is shown 
in Fig. 123 We see that it consists only of short elec- 
tron loops, which means that the charge on each step is 
transferred by different electrons as it should be in the 
inelastic process. 

All intermediate Green functions with the grain index 
i = 0, ...N enter the diagram in Fig. [SHlbeing integrated 
over the internal states of the grain, such that each Green 
function is 

G(r)^^G,(r)^ (Bl) 
^ sm(7rJ r) 

The processes for the electrons located on the initial 
and final sites (iq and im) are given by the Green func- 
tions that are accompanied by the phase correlations 
n(T) exactly in the same way as in the case of the elastic 
cotunneling and, thus, can be described in terms of the 
Green functions G. 

The intervals tq — Tq and r^v — t'j^ coincide with each 
other within the accuracy of the inverse Coulomb energy 
and we just put tq — Tq ~ t^v — = r. In the same 
limit, the integrals over the intermediate times intervals 
Arfe = Tfc — Tfc-i (and over — r^_]^) can be calculated 
independently resulting in the contribution 

/oo 
d(Arfe) e-^^l^^'=l+'^'='^^'^' 2/Ek, (B2) 
-oo 

where the energy E^ is defined in Eq. H2.107|l . Collecting 
all the terms for the probability Pin one obtains 
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where Ae = en+i — £o is the difference between the en- 
ergies of the initial and final states and E and g are ge- 
ometrical averages along a the tunneling path (see Eqs. 

mm)- 

Finally, in order to find the tunneling probability Pin 
for the inelastic cotunneling one has to make the analyt- 
ical continuation in Eq. ljB3|) to the real times — ity^ 
and integrate over arriving at 



p — 



wg 



N+l 



+00 



dt 



2TTiT 



E sinh(7rrt) 



2N 



-itAe 



(B4) 



Here the singularity of the function sinh^^^(7rT't) is as- 
sumed to be in the upper half of the complex plain. At 
zero temperature one can easily calculate the integral 
over the variable t in Eq. (|B4p by shifting the contour of 
integration in the complex plain either to ~ioo for posi- 
tive Ae or to -|-ioo for negative Ae. 

In the first case we obtain Pin ~ 0. This refiects the 
fact that the real tunneling process with an increase of 
the energy of the electron is forbidden at T = 0. In the 
latter case, Ae < 0, the zero temperature probability 
is determined by the pole of the function sinh^^^(7rTt) 
that results in Eq. (|2.118|l . 

At finite temperatures the integral in Eq. ljB4|l can be 
expressed in terms of the Euler Gamma functions that 
leads to the general expression H2.114|) . 

Generalization of the result to include the long-range 
interaction consists, as in the case of the elastic co- 
tunneling, of a proper redefinition ljA12|) of the elec- 
tron/hole excitation energies. Considerations analogous 
to those presented in the previous appendix allow us to 
conclude that inclusion of the off diagonal terms in the ca- 
pacitance matrix results only in the appearance of the co- 
efficient under the logarithm (Bcloborodov et al., 2005^ 
in the expression for the localization length (|2.117() . 
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